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Abstract 



This thesis examines the applications and determinations of meson light-cone distribution 
amplitudes, which enter the theoretical description of exclusive processes at large momen- 
tum transfer. The investigation of such processes, in the context of B physics, provides 
one with a rich and extensive way of determining the Standard Model parameters of the 
CKM matrix, which are essential in describing CP violation, and searching for tell-tale 
signs of new physics beyond the Standard Model. 

We investigate the twist- 2 and twist-3 distribution amplitudes of vector mesons and fully 
examine SU(3)p-breaking effects and include leading G-parity violating terms. We use 
the conformal expansion allowing the distribution amplitudes to be described by a set of 
non-perturbative hadronic parameters which is reduced by invoking the QCD equation 
of motion to find various interrelations between the distribution amplitudes. Numerical 
values of the leading non-perturbative hadronic parameters are determined from QCD 
sum rules. 

The new distribution amplitude results find direct application in the radiative B decays 
to light vector mesons B — > V7. We examine the phenomenologically most important 
observables in this decay mode using the formalism of QCD factorisation in which the 
distribution amplitudes play a vital role. We also include long-distance photon emission 
and soft quark loop effects, which formally lie outside the QCD factorisation formalism. 
The analysis encompasses all the relevant modes, that is B u ^ — > p, ui, K* and B s — > 0, K* . 

We also calculate the B — > rft' transition form factor using QCD sum rules on the light- 
cone. The method relies on the collinear factorisation of the QCD dynamics into a pertur- 
batively calculable hard-scattering kernel and the non-perturbative universal distribution 
amplitudes. We include the singlet contribution originating from the U(1)a anomaly and 
bring the calculation consistently within the 77-77' mixing framework. 
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Introduction 



One only has to ask the question "why?" a handful of times before one reaches the answer 
"I don't know" , regardless of the topic considered and regardless of the person asked. It 
is safe to say, however, almost all questions of the structure of matter at the smallest of 
distances leads one directly to, or at least through, the field of modern particle physics. 
The beginnings of our understanding of the physical world harks back to the dawn of 
scientific reasoning in the ancient world; logic and reasoning were applied with the aim 
of describing the behaviour of physical systems in terms of simple universal axioms, a 
philosophy which still holds strong today. Through experimentation and the language 
of mathematics the scientific method has driven back the edge of ignorance to frontiers 
unimaginable to those physicists of 100 years ago, let alone the natural philosophers of 
millennia ago. The present "coal face" is known as the Standard Model [1, 2] which 
describes three of the four known forces of nature - electromagnetism, and the weak and 
strong nuclear forces - in one unifying framework. 

Frustratingly, the Standard Model does not explain many of the things which it encom- 
passes; it does not provide an origin for CP violation but only gives a parameterisation, 
nor does it explain why there are three generations of quarks and leptons, or their hier- 
archy of masses. All attempts to bring gravity into the fold have so far failed, however, 
whatever theory lies beyond must yield the Standard Model as some limiting case. 

The Standard Model has been scrutinised relentlessly since its inception. Remarkably, 
nearly without fail it has held its ground over the entire breadth of its theoretical reach 
and so the task of finding new ways to probe its structure requires ever more the creativity 
and ingenuity of both theorists and experimentalists alike. Novel experimental signatures, 
against which to pit theory, must be used to maximum potential. From a theoretical 
standpoint there are still many challenges to be met, especially in preparation for the 
next generation of collider experiments now just round the corner. Particularly, the 
control and reduction of the theoretical uncertainty of Standard Model predictions is of 
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paramount importance as only then can one hope to be in a position to discern signs of 
new physics from that of the Standard Model background. 

Some of its most challenging tests of the Standard Model fall in the field of heavy-flavour 
physics, within which B physics has proven itself to be rich and fertile. Today it is an 
area of high activity with many success stories, including the recent measurement of the 
B° s -B° s mass difference Am s at the Tevatron [3]. Moreover, two dedicated "5-factories" , 
Belle at KEK [4] and BaBar at SLAC [5], have measured a range of observables, such as 
branching fractions and CP asymmetries, of a vast number of B decay modes. Looking 
to the future, the B physics community eagerly await the forthcoming LHCb experiment, 
and beyond that so-called "superflavour factories" [6] have been championed with the 
aim of probing rare B decays to extract CP violation parameters to much higher levels 
of accuracy. It is imperative to find tests of the Standard Model which may be observed 
in these up-and-coming experiments [7] and promising modes include the rare decays 
B -> V-i and B -> Kfi+ti'. 

The strict pattern of CP violation of the Standard Model finds its origin in the Cabbibo- 
Kobayashi-Maskawa (CKM) matrix [8,9]. CP violation was discovered in B physics via 
the decay mode B® — > J/ipKg and found to be large, in contrast to K decays where the 
violation is tiny. The possible largeness of CP violation in B decays offers promising ways 
to detect new physics indirectly via CP violating observables testing the CKM paradigm. 

Theoretically, central to the description of B decays is the disentanglement of the weak 
decay process from strong interaction effects leading to a low-energy effective Hamilto- 
nian in which the physics at a scale 0(M\y) is well under control. Achieving this goal for 
the wide range of B decays of interest has only been possible through huge calculational 
effort; the availability in the literature of Wilson coefficients at next-to-leading-order, and 
in some cases next-to-next-to-leading-order, is testament to this. Furthermore, the theo- 
retical description of the matrix elements of effective B decay operators has been hugely 
improved through QCD factorisation methods. We discuss and make use of one such 
framework, namely that introduced by Beneke, Buchalla, Neubert and Sachrajda [10-12]. 
The so-called BBNS approach showed, to leading-order in a l/m^ expansion, that the 
a s corrections beyond naive-factorisation of a large class of non-leptonic B decay matrix 
elements are calculable in terms of B transition form factors and meson light-cone distri- 
bution amplitudes. Armed with the corresponding amplitudes the phenomenologist may 
construct observables, such as branching ratios, CP asymmetries and isospin symmetries, 
which may then be compared to experiment. The predictive power of the QCD factori- 
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sation framework is jeopardised by a poor understanding of both these non-perturbative 
QCD quantities and the impact of the generally unknown power-suppressed contributions 
0(1/ rub); this in part motivates the work of this thesis. 

In this thesis we investigate SU(3)p-breaking effects in vector meson distribution ampli- 
tudes which are crucial in differentiating between the particles p, K* and <fi. The leading 
non-perturbative DA parameters are determined via the method of QCD sum rules intro- 
duced by Shifman, Vainshtein and Zakharov [13-15]. The method provides a prescription 
for the systematic calculation of non-perturbative QCD parameters, albeit with an irre- 
ducible error ~ 20 — 30%, and constitutes an extremely useful theoretical tool. 

The sum rule results have a direct application in the QCD factorisation description of B 
decays to p, K* and mesons. In particular, radiative B decays to vector mesons B — * V7, 
are an excellent example of a process potentially sensitive to new physics contributions, as 
at leading order the decays are mediated at loop level in the Standard Model. We perform 
a phenomenological analysis of these decays using the QCD factorisation framework of 
Bosch and Buchalla [16, 17] including leading power-suppressed corrections for which the 
updated non-perturbative distribution amplitude parameters find use. The impact of the 
power-suppressed corrections on the key decay observables is discussed and leads to a 
better understanding of the theoretical uncertainty of the QCD factorisation predictions. 

Also, we calculate important contributions to the B — > 77 ^ transition form factors via a 
variant sum rule approach, known as light-cone sum rules, for which distribution ampli- 
tudes play a crucial role. The result of the analysis elucidates a major source of theoretical 
uncertainty of the B — > 77W form factor. The result impacts B — > K*r)('\ for example, 
where the experimental data and QCD factorisation predictions of the branching ratios 
are inconsistent. 

The thesis is structured as follows: 

• Chapter [1] introduces some of the fundamentals of the Standard Model and its 
application to B physics. We define the QCD Lagrangian and the CKM matrix, 
introduce CP violation in Standard Model B decays, and briefly discuss the structure 
of the AB = 1 weak effective Hamiltonian. 

• Chapter [2] covers the definitions of the light-cone distribution amplitudes of the light 
vector mesons p, K* and (p. We determine their structure up to twist-3 accuracy and 
using the conformal expansion and QCD equations of motion express the distribu- 
tion amplitudes in terms of a finite set of non-perturbative parameters. We extend 
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previous determinations in order to fully differentiate between the three particles 
by including all G-parity violating contributions and SU(3)p-breaking effects. 

Chapter [3] discusses the QCD sum rule method and its extension light-cone sum 
rules. The methods allow, amongst other things, the determination of the non- 
perturbative distribution amplitude parameters and transition form factors respec- 
tively, and are very widely applicable in and beyond B physics. 

In Chapter H] we apply QCD sum rules to determine the leading non-perturbative 
distribution parameters defined in Chapter [2J Consistency requires the inclusion of 
all G-parity violating contributions and SU(3)F-breaking effects to the sum rules, 
and we extend previous determinations by including higher-order strange quark mass 
effects and 0(a s ) contributions to the quark condensates. We analyse the resulting 
sum rules and provide updated numerical results for all parameters. The results 
of this section find immediate application in QCD factorisation and light-cone sum 
rule descriptions of processes involving these vector mesons. 

In Chapter [5] we calculate the gluonic flavour-singlet contribution to the semilep- 
tonic B — > t)V> transition form factor in the framework of light-cone sum rules. In 
doing so we discuss pseudoscalar meson and two-gluon distribution amplitudes. The 
new contribution is combined with the previous determination of the quark contri- 
bution, to complete the theoretical treatment of these form factors. The rft' system 
is complicated due to large mixing effects via the U(1)a anomaly. We introduce 
the phenomenological framework of 77-77' mixing and connect it to the form factor 
calculation in a consistent manner. The results of this chapter find immediate appli- 
cation in the QCD factorisation description of B — > 7]^ transitions, which in turn, in 
principle, allow a determination of the CKM matrix element \V u t,\ from B — > rj^'lv. 

Chapter [U] introduces the framework of QCD factorisation, which is an important 
application of meson distribution amplitudes and transition form factors. We briefly 
discuss the BBNS approach and then go on to discuss the leading contributions to 
QCD factorisation in the context of B — > decays. 

In Chapter [7] we investigate the impact of the relevant, power-suppressed contribu- 
tions to B — > V'-f beyond the QCD factorisation formula. We include long-distance 
photon emission from weak annihilation diagrams and soft gluon emission from 
quark loops. The non-perturbative distribution amplitude parameters determined 
in Chapter H] find use in a light-cone sum rule estimation of the latter. The key 
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observables are the branching ratios, isospin asymmetries and the indirect time- 
dependent CP asymmetry S(V~f) which, as has been know for some time, forms 
the basis of a "null test" of the Standard Model. Assuming no new physics contri- 
butions, we extract the ratio of CKM matrix parameters \V t d/V t d\ to a competitive 
degree of accuracy. 

• We summarise and conclude in Chapter [HI 

The material of Chapters [2] and H] follows Ref . [40] and the material of Chapters [5] and [7] 
follows Refs. [65] and [70] , respectively. Some of the more bulky equations, and material 
not necessary in the general flow of reading the thesis, are given in two appendices. 
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Chapter 1 



Fundamentals Of B Physics 



In this chapter we begin with the basics of the Standard Model and then go on to discuss 
two concepts which are central to the investigations of B physics, and those of this thesis: 

• CP violation in the flavour sector, which follows a strict pattern in the Standard 
Model and can readily be sensitive to new physics; 

• the AB = 1 effective weak Hamiltonian, which we briefly discuss as it is the starting 
point of many phenomenological studies in B physics. 



The Standard Model (SM) [1,2] is a model of great scope and predictive power. Despite 
its successes, however, we know it to be incomplete; for example, the recent discovery of 
neutrino oscillation and the lack of conclusive evidence for the Higgs particle providing 
two areas of intense theoretical and experimental effort. The SM describes three of the 
four known fundamental forces of nature; the strong force, the weak force and electromag- 
netism. Quantum Chromodynamics (QCD) is a Yang-Mills gauge theory based on the 
gauge group SU(3) and describes the fundamental interactions of the strong interaction 
as interactions between quarks and gluons [18-21]. The basic QCD Lagrangian is 



1.1 The Standard Model 




(l.i) 



i 



with 



Sijd^ - ig 8 (t a )ijA 



G 



•a 



d,A a v - d v A% + g s r hc A\Al , 



(1.2) 



where the sum is over all quark flavours q, i,j = {1, 2, 3} are colour indices, the t d are the 
3x3 colour matrices with d = {1, . . . , 8} and f abc are the structure constants. G^ u is the 
gluonic field strength tensor, and A^ is the gluon field. We will make use of the notation 
{G^ u )ij = G^ u (t a )ij and the relation g 2 s = 4na s (and e 2 = 47raQ ED ). The Lagrangian can 
alternatively be defined with the replacement g s — > —g s and the sign convention matters 
for the applications in Chapters H] and [71 

The non-Abelian nature of QCD leads to the possibility of gluon self-interaction and the 
celebrated asymptotic freedom property of QCD [20-23]. The coupling tends to zero, 
giving a theory of free quarks, at asymptotically high energy. On the other hand, at low 
energy, or large distances, the coupling increases. At energies for which a s > 1 pertur- 
bation theory is not applicable, and one has to resort to non-perturbative methods to 
determine the effects of QCD. Despite the simplicity of the QCD Lagrangian (11. ip an 
accurate determination of non-perturbative QCD from first principles, and hence con- 
finement, poses a major challenge. One such method, based on ideas of Wilson [24], is 
that of Lattice QCD, which aims to calculate the QCD action computationally on a grid 
of discretised spacetime points. An altogether different, and less rigourous, method is 
that of QCD sum rules, which encodes non-perturbative effects in terms of non- vanishing 
vacuum expectation values of operators with the quantum numbers of the vacuum. This 
method is central to the work in this thesis, and shall be discussed in Chapter [31 

The electroweak force is the unification of the weak nuclear force and electromagnetism 
given by the Glashow-Salam- Weinberg model. The model is based on the gauge group 
SU(2)l ® U(1)y, which is broken by spontaneous symmetry breaking to yield U(1)q - the 
gauge group corresponding to Quantum Elecrodynamics (QED). The weak interaction is 
mediated by three massive gauge bosons W and Z° and occurs between quarks and 
leptons. The quarks and leptons are arranged, within the three generations, into left- 
handed doublets and right-handed singlets under SU(2) L 



where the weak eigenstates U = {u,c,t}, D = {d,s,b} and l~ = {e~,/i - ,r - } are the 
up-type quarks, down-type quarks and charged leptons respectively. The subscript L (R) 
represents the left (right)-handed projectors gL(R) = |(1 T 7s)9 which reflect the chiral 




(1.3) 



2 



nature of the weak interaction. The neutrinos are massless in the SM, and the right 
handed neutrino does not exist. The electroweak interactions of the quarks are described 
by the following Lagrangian, which consists of a charged current (CC) and a neutral 
current (NC) 

The neutral current part of the Lagrangian is made up of the electromagnetic current J^ m 
and neutral weak current J?,: 

j; m = Qu U Llfl U L + Q D D l7a ,D l , = ^(f/ L 7^L - A^l) , (1.5) 

where Qu(D) — 2/3 (—1/3) is the electric charge of the U (D) quarks, 6w is the weak 
mixing angle and g is the electroweak coupling related to the electromagnetic coupling 
by e = gsinOw- Rotating to the basis of mass eigenstates modifies the charged current 
in the quark sector to 

Jt = U^V CKM D™, (1.6) 

where Vqkm is the Cabbibo-Kobayashi-Maskawa matrix [8,9] and the superscript m denotes 
mass eigenstates. The CKM matrix is 3 x 3 (for three quark generations), unitary, and 
its off-diagonal entries allow for transitions between the quark generations. There are no 
flavour- changing neutral-currents (FCNC) at tree-level in the SM as the neutral currents 
J° m and are invariant under the transformation to the mass eigenbasis, which is known 
as the Glashow-Iliopoulos-Maiani ( GIM) mechanism [25]. The entries of the CKM matrix 
are written as 

/ V u d V us V u b \ 



Vc 



CKM 



;i-7) 



V c d Vcs V c b 

\ v td v ts v tb ) 

and are fundamental parameters of the SM that have to be determined from experiment. 
Evidently, the matrix has n 2 = 9 parameters n(n — l)/2 = 3 of which are rotation angles 
due to its unitarity. The six quark fields in Eq. fll.3p can be re-phased, up to an overall 
phase, leaving the Lagrangian invariant and therefore 9 — 5 — 3 = 1 phase remains giving 
rise to complex entries - complex coupling constants. This is the origin of CP violation in 
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the quark sector of the weak interaction. The leptonic sector is described by an analogous 
mixing matrix which, in the absence of neutrino masses, is given by the unit matrix 
because all phases can be rotated away. 

The CKM matrix (11.71) is often parameterised to incorporate the constraints of unitarityQ 
A very useful and convenient parameterisation is the Wolfenstein parameterisation [26] 
which, along with unitarity, incorporates the experimental observations |V^ S | <C 1, \V cb \ ~ 
\V US \ 2 and \V u b\ *C \V c b\. It is an expansion in A = |T4s| ~ 0.22, and as such is only 
approximately unitary at a given order in A: 



V, 



CKM 



2 



-A 



A A\ 3 (p - i?7) \ 



A 2 



\ A\ 3 (l- p-iri) -AX 2 



AX 2 
1 



+ 0(X 4 



J 



The matrix is given in terms of the four parameters (A, A, p,r)); A and p 2 + rf are order 
unity and the hierarchy of sizes of elements can be infered from the powers of A. The 
smallness of V cb and V ub are responsible for the relatively long lifetime of B mesons (and 
baryons), which facilitates their experimental detection. The unitarity of the CKM matrix 
gives six equations that equal zero and can be represented as triangles in the complex 
plane. The most widely used of these relations in B physics is 



v ud v: b + v cd v: b + v td v; b = o 



(1.9) 



which is invariant under phase transformations and is an observable. The above relation 
is divided by V c dV* b to give a triangle in the complex plane with a base of unit length and 
upper apex at the point (p, fjj^ known as The Unitary Triangle (UT), see Figs. [Tj] and 
11.21 The sides of the UT are given by 



R b 
Rt 



\Vu d v: b \ 
\v cd v: b \ 

\Vt d V t l\ 

\v cd v; b \ 



— 9 



2 ) X 



^(l-py + n 2 . 



v, 



ub 



V, 



cb 



1.10) 

1.11) 



1 The "standard" parameterisation of the CKM matrix is in terms of the three mixing angles 9ij 

= 1, 2, 3) and the CP violating phase S [27]. 
2 The following rescaling proves convenient to the definition of the UT: p — > p = p(l — A 2 /2) and 
V V = V (1 - A 2 /2). 



1 



The angles are given by 



a 



arg (-V td V ub V t lV: d ) , = arg (-VM,V^Q) , 7 = arg (-V^fcKTO) • (1-12) 



The (over) determination of the sides and angles of the UT is a major quest in understand- 
ing the SM. To achieve this goal one must construct decay observables, which can then be 
matched to experimental results in order to extract values for the desired CKM (or equiv- 
alently UT) parameters. Such observables include branching ratios, which may appear 
simply proportional to a CKM matrix element, and CP asymmetries, which encode the 
effects of the SM predictions of CP violation, and can also be measured experimentally. 




Figure 1.1: The Unitary Triangle. The determination of the sides Rb and Rt and the angles 
a, (3 and 7 lead to stringent tests of the Standard Model. 



lE=r 



-0.5 - 




Figure 1.2: Constraints on the angles a, (5, and 7 and sides Rb and Rt of the Unitarity Triangle 
as imposed from numerous experimental sources. Complied by the CKM fitter group [31]. 
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1.2 CP Violation In B Decays 



Does the CKM matrix (11.71) account for the CP violation observed in nature? Examining 
CP violation in B decays allows one to probe the structure of the CKM matrix and is 
a very promising way to detect the effects of new physics, which many not be expressed 
through other decay observables. Consequently, the CP properties of FCNC processes, 
which are characterised by their potential sensitivity to new physics effects, have been 
under intense theoretical and experimental investigation for many years. Prime examples 
of such processes include B°-B° mixing (see for example Ref. [28]) and radiative B decays, 
see Chapter [71 

The idea that the weak interaction may violate parity was first suggested many years ago 
by Lee and Yang [29], and quickly confirmed in the (3 decay of 60 Co by Wu et al. [30]. The 
violation of the combined CP symmetry was first observed in the context of K decays in 
1964 [32] and it was not until 2001 that it was first observed outside the K system in 
B® — > J/ifiKg decays [33,34]; in both cases the CKM paradigm was upheld. Recently 
discoveries in B physics include the measurement by CDF of the mass difference Am s [3]. 
Some of the most important sources of information about the UT from B physics include: 
the determination of sin 2(3 from the "gold-plated" decay B — > J/ipKs; the extraction of 
a from non-leptonic B decays such as B — > tttt; the extraction of | Vtd\/\ Vts\ from B mixing 
and radiative B decays, such as B — > and the determination of \V u b\ from B — > tiIu. 

The Bq-B® systems, where q = {d, s}, exhibit the phenomenon of particle-antiparticle 
mixing, which, in the SM is mediated by so-called box diagrams whose amplitudes are 
~ G 2 F and therefore very small. We do not go into any detail about the theory of neutral 
state mixing and we restrict ourselves to only the formulas required in this thesis; for more 
information see Refs. [35,36]. State mixing causes, for example, an initially pure beam 
of B° mesons to evolve into a time-dependent linear combination of B° and B° mesons. 
There are four main quantities that describe the B q -B q system and its decays: the width 
difference Ar g , the mass difference Am q , the CP violating mixing phase <p q and Xf (not 
to be confused with the Wolfenstein CKM parameter A ~ 0.22). One begins by writing 
the heavy (H) or light (L) eigenstates of evolution in terms of the flavour states: 

\B U ) = p\B°) - q\B°) , \B L ) = p\B°) + q\B°) , (1.13) 



6 



with \p\ 2 + \q\ 2 = 1. The ratio q/p is given in terms of the B®-B® mixing matrix Mf 2 , by 



12J 



Mf 2 



-i(j>q 



;i.u) 



under the condition Ar g <C Am r Experimentally, there is no evidence for mixing-indiced 
CP violation in the B^-B^ systems, i.e. \q/p\ ds ~ 1 [37]. The CP violating mixing phase 
is given by <p q = arg [Mf 2 ] which in the SM and the Wolfenstein parametrisation of the 
CKM matrix can be written in terms of the UT angles as 



<> d = arg[(^ 6 ) 2 ]=2/3 ; 



axg[(V£Vfe) 2 ] 



-2A 



V, 



ub 



V, 



cb 



S1I17 . 



(1.15) 



Besides mixing-induced CP violation there also exists direct and indirect CP violation for 
B and B decays to a common CP eigenstate /. The corresponding time- dependent CP 
asymmetry is given by 



A C p(t) 



rffift) -> f) - r(Bg(t) -> /) 
r(B0(t)^f) + r(B0(t)^f) 

S(f) sm(Am q t) - C(f) cos(Am (? t) 



1.16) 



indirect 



direct 



where we have neglected the width difference AT q = 2Re [M^T^] /| Mf 2 \- The oscillation 
frequency is set by the mass difference between the heavy and light states 



Am, = m q H — m 9 L = 2 1 M\ 2 \ , 



;i.i7) 



and the current world averages are [37]: 

Am d = 0.507 ± 0.004 ps^ 1 , 
Finally, if we define the observable quantity 



stat. sys. 

Am, = 17.77 i^itHTTps- 1 



q A 
pA 



'1.181 



1.19) 



where A denotes the decay amplitude, then the two CP asymmetries can be written as 

2 MA,] 



C(f) 



l-\X f \ 2 
1 + |A,| 2 



S(f) 



1 + |A,| 2 ' 



(1.20) 
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1.3 Effective Field Theories Of Weak Decays 



A very widely used tool in the theoretical description of B decay processes is the framework 
of effective field theories [35,38]. The framework simplifies the dynamics of the weak decay 
by relying on an operator product expansion (OPE) [39] of the weak vertices to separate 
the short and long distance physics. The OPE yields a concise effective Hamiltonian 
H, e ff built from a set of local effective operators Qi multiplied by renormalisation-scale 
dependent perturbatively calculable Wilson coefficient functions C;(/i): 

{H) (H e - ff ) ~ X)CiO*) (Qi) + 0(k 2 /M 2 w ) , (1.21) 

i 

where k is the momentum flowing through the W boson propagator. The separation of 
energy scales stems naturally from the fact that the weak decay of the B meson is governed 
by physics originating at well separated scales: m t , M w 3> m b c ^> A QCD ^> m u ^ s . It is 
the interplay of weak and strong effects that complicates the treatment of these decays, 
and must be dealt with appropriately. By taking into account radiative corrections to tree- 
level and penguin diagrams, ultimately one obtains the effective Hamiltonian in terms of 
the set of all relevant local operators, which is closed under renormalisation. The full 
AB = 1 effective Hamiltonian is, for a final state containing a D quark 



neff = % £ x(c 

V 2 U=u,c 



CiQi + C2Q2 + CV7Q77 + CsgQgg + CiQi 



i=3,...,10 



1.22) 



where make use of the standard short-hand notation for the product of CKM matrix 
elements = Vj} D Vub- The form of Eq. fl 1.22ft is chosen by assuming the unitarity of 
the CKM matrix (11.91) to explicitly remove the dependence of the top quark CKM matrix 



elements which originate from penguin loops. The effective operators are 
Current — Current^ : 



(DiUjv-AiUjtyv-A , Qi = (DU) v _ A (Ub) v _ A 



QCD Penguin 



Q 3 = 
Q 5 - 

Electroweak Penguin 

Qi - 
Q9 - 

Electromagnetic Dipole : 

Q77 



(Db) v ^ A ^2(qq) v _ Aj 

g 



Q 4 = (D i b j ) v _ A ^2(q j q i ) v _ A 



V+A , 



Q 6 = (D i b j ) v _ A ^2(q j q i ] 
g 



V+A , 



{Db) v ^ A 2 e g(^)y+A > Qs = (A6j)v-a Y n e ^Qi)v+A , 
g g 
3 3 

(Db) v - A Y n e g(q<l)v-A , Qw = (Dibj) v ^ A Y ^giwdv-A , 



-^m b Da^(l + 75)^ b + Da>^(l - l5 )F^ b . 

87T Z o7T 



Chromomagnetic Dipole : 



s m.Da^il + ^G^b+^mDDa^il-j^G^b, (1.23) 



>7T 



8tt 2 



where e q is the electric charge of the quark q in units of |e| and is the photonic 
field strength tensor. The Wilson coefficients entering the effective Hamiltonian are es- 
sentially effective coupling constants of the local effective operators. One can view the 
renormalisation of the matrix elements as an equivalent renormalisation of their Wilson 
coefficients. One makes use of renormalisation-group techniques to sum the potentially 
large logarithms ~ In Myy / fi 2 that appear naturally in the evolution from weak scales 
O(Myy) to hadronic scales, such as /J ~ m&. The operators (I1.23P mix with each other 
under evolution and from the renormalisation-scale invariance of 7i e ^ one finds 



1.24) 



where 7 is the anomalous dimension matrix, which can be given as an expansion in the 
strong coupling via the renormalisation constant Z 



-1 dZkj 



7 



fOLs_ 

Uvr 



7 



(0) 



a 

Air J 



V 7 w 



0{al) . 



1.25) 
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Solving Eq. (11.24j) yields the evolution of the Wilson coefficients via the evolution matrix 



Ci(fJL) = Uij(fJL, fio) Cj(n ) , U (//, no) = exp / dg 



9 Oo) 



;i.26) 



where (3{g) is the QCD /^-function. To leading order one has 



v(0)T 



gaOfo) ^ 2,30 

a s (/j) 



y- 1 



1.27) 



J D 



where is the matrix that diagonalises and 7^°) is a vector of the eigenvalues of 
the leading order anomalous dimension matrix 7^ = Vjp V~ x . At NLO we have 



c i 0i) = c i (o) 0i) + ^cj 1 V) 



4ir 



1.28) 



and the evolution is a bit more complicated: 

a s {n) 



1 + 



Air 



J 



£/ LC Wo) 



with 



J = vsv~ l 



S. - x ..„,(°) & 



_ a a (/xo) j 

47T 



1.29) 



2/?o 2 2/30 + Tf-Tf 



G = 1/ _1 7 (1)T V . (1.30) 



To NLO the required /3-function coefficients are (3 X = fN* - fN c N f - 2C F N f and 



0o 
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Nr. 



Nf with iVy is the number of active flavours, CV = (iV^ — 1)/(2N C ) and iV c 



the number of colours. Care must be taken in evolving through "thresholds" where the 
number of active flavours Nf changes; the evolution must then be taken in stag es, as a 
change in Nf changes the /3-function coefficients and the anomalous dimension matrices. 
If there is a flavour threshold /4h between /xq and /i, which changes the number of active 
flavours from Nf to Nf + 1, then one has to make the replacement 



U(fi,fi ) —> £/(/i, /i t h) 



Nf+l 



'1.3V 



The effective Hamiltonian, combined with the renormalisation-group improvement of the 
perturbative series forms an exceptionally powerful framework. The matrix elements of 
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the local operators (Qi) are the subject of QCD factorisation theorems, such as that 
discussed in Chapter [HI which allow the calculation of B decay amplitudes. From these 
amplitudes one can construct observables such as branching fractions, CP asymmetries, 
and isospin asymmetries which can be investigated phenomenologically. 
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Chapter 2 



Vector Meson Light-Cone 
Distribution Amplitudes 



In this chapter we discuss light vector meson light-cone distribution amplitudes and via 
the (approximate) conformal symmetry of QCD present expressions for the distribution 
amplitudes up to twist-3. The method introduces a set of non-perturbative parameters 
which is reduced in size by invoking the QCD equations of motion to relate the two- 
particle twist-3 distribution amplitudes to the three-particle twist-3 and two-particle twist- 
2 distribution amplitudes. In our analysis we include all SU(3)F-breaking effects and 
G-parity violating terms thus allowing one to fully differentiate between p, K* and </> 
mesons. Moreover, a non-zero quark mass induces a mixing between twist-2 and twist-3 
parameters under a change of renormalisation scale p. To simplify notation we explicitly 
consider the K* meson, with quark composition sq where q = {u, c/}0 




There are two main applications of meson distribution amplitudes that motivate their 
study: 

• they are directly applicable to the theoretical description of exclusive decay processes 
via QCD factorisation theorems, which require the distribution amplitudes as a non- 
perturbative input, see Chapter El 

1 The notation in this thesis, K* being a (sq) bound state, is in contrast to the standard labelling, 
according to which A"* = (ds) and K*° = (sd). This is the standard notation used for light-cone 
distribution amplitudes where K* always contains an s quark, and K* an s quark. This distinction 
is relevant because of a sign change of G-odd matrix elements under (sq) <-> (qs). This notation also 
applies to calculations of form factors and other matrix elements which involve light-cone distribution 
amplitudes. 
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• they are also applicable to the determination of transition form factors from the 
light-cone sum rule approach and as such are indirectly applicable to the same 
QCD factorisation theorems for which the transition form factors are also required, 
see Chapters |3] and El 

In Chapter H] we calculate, from QCD sum rules, numerical values for the leading twist-2 
and twist-3 distribution amplitude parameters defined here. Standard notations used, 
such as the light-cone coordinates, are given in Appendix |A] The material covered in this 
chapter partially follows that of Ref . [40] . 



2.1 Introduction 

Hadronic light-cone distribution amplitudes (DAs) of light mesons were first discussed in 
the ground-breaking papers of Brodsky, Lepage, and others, see Refs. [41-48] and play an 
essential role in the QCD description of hard exclusive processes [49,50]. The amplitudes 
that describe such processes factorise in the asymptotic limit Q 2 ~ 1/x 2 — ■> oo - where 
Q 2 is the momentum transfer and x the transverse separation of the partons - and are 
dominated by contributions from near the light-cone. The factorisation is given by the 
convolution of a hard-scattering kernel, calculable in perturbation theory, and process- 
independent, universal, non-perturbative DAs. 

The study of hadronic DAs has a long history. The simplest and first to be investigated 
were the twist-2 DA of the tt [43,44,46,47]. Higher twist DAs of the n, alongside those 
of the other pseudoscalar mesons followed [51]. For vector mesons, the leading-twist DAs 
of the p were first investigated by Chernyak and Zhitnitsky in Ref. [52] and later in 
Refs. [53,54]. The formalism of higher twist-3 and twist-4 contributions, including meson 
mass corrections, was investigated by Ball et al. in Refs. [55-58]. 

The DAs of the K* (K) differ to those of the p (it) due to the non-zero strange quark 
mass which yields SU(3)F-breaking and G-parity violating corrections from a number of 
different sources,^] The study of the various contributions span many publications: 



explicit quark mass corrections to DAs and evolution equations are generated by 
the QCD equations of motion (EOM) and only affect higher twist DAs. The con- 
tributions for vector mesons were calculated in Ref. [55] up to twist-3, and those to 



2 Perfect SU(3)f symmetry is realised for equal u, d, and s quark masses. 
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the evolution equations for vector mesons in Ref. [40] and flavour-octet pseudoscalar 
mesons Ref. [57]. 



• G-parity violating contributions, which are proportional to m s — m q and hence 
vanish for equal quark masses, i.e. for p and 0, were investigated for twist-2 DAs in 
Refs. [52,55,59-62] and for twist-3 DAs in Ref. [40]. 

• SU(3)F-breaking of non-perturbative hadronic parameters entering the DAs. The 
effects for the twist-2 parameters are known from Refs. [52,55,59], twist-3 from 
Ref. [40] and twist-4 from Ref. [58]. The twist-3 vector meson parameters are 
discussed in Chapter H] where we include all these effects in a determination of 
numerical values using QCD sum rules. 

The objects which define the DAs are vacuum-to-meson matrix elements of non-local 
operators at strictly light-like separations z 2 = [52]. Two examples we shall encounter 
are 

(0\q(z)r[z,-z}s(-z)\K*(p,X)), (0\q(z)[z,vz}g s G liV (vz)r[vz,-z}s(-z)\K*{p,\)), 

(2.1) 

where T is a general Dirac matrix, A = {||, _L} is the polarisation of the K* meson and the 
quark fields are taken at symmetric separation for simplicity]^ The first (second) matrix 
element above corresponds to a two- (three-) particle Fock state. To render the matrix 
element gauge invariant the path-ordered gauge factor is included 



[x, y] = P exp 



i 



(2.2) 



ig s / dt(x-y) ll A fi (tx + (l-t)y) 
Jo 

For convenience we work in the fixed-point gauge 

(x-x o )"Afa) = 0, (2.3) 

and by choosing x = we have [x, — x] = 1. The gauge factor will be implied unless 
otherwise stated. The DAs are dimensionless functions of the collinear momentum frac- 
tions of a fixed number of constituents within a meson, at zero transverse separation. For 



3 The Dirac matrices T = {<j^ Vi 275, 1} give rise to so-called chiral-odd distributions because they are 
chirality- violating. Likewise, distributions generated from T = {7^, 7^75} are chiral-even. 
4 also known as the Fock-Schwinger gauge. 
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two-particle DAs the constituent strange quark and antiquark (q) share u and u = 1 — u 
of the meson momentum p respectively. For three-particle DAs we have a = (a\, a2, C13) 
corresponding to the momentum fractions carried by the strange quark, antiquark (q) and 
gluon, respectively. For a minimum number of constituents, the DAs are related to the 
Bethe-Salpeter wavefunction (p BS by integration over the transverse momenta 



where /x is the renormalisation scale. The price to pay for integrating out k± below 
/j is a renormalisation-scale dependence of the DAs governed by renormalisation-group 
equations. The DAs have to be evaluated at the scale /i 2 ~ x~ 2 i.e. of the order of the 
deviation from the light-cone [63]. 

Non-local operators that appear at finite Q 2 or mass scales are expanded near the light- 
cone x 2 7^ as an OPE in terms of the renormalised non-local operators on the light-cone 
- the light-cone expansion [63] If] After taking matrix elements the resulting Lorentz- 
invariant amplitudes are matched to the definitions of the DAs with the coefficient func- 
tions of the expansion taken at tree-level, to leading logarithmic accuracy. 

The structure of vector meson DAs follows the same pattern as the nucleon structure 
functions and can be classified in the same way [64]. They are described by separate DAs 
for each polarisation and thus there are more vector meson DAs than pseudoscalar DAs. 

Lastly, we briefly mention some other DAs. Flavour-singlet pseudoscalar meson DAs are 
complicated by the U(1)a anomaly of QCD and are discussed in Chapter in the context 
of the B — > 77 transition form factor [65]. Much work has been done concerning the 
DAs of heavy mesons, such as the B meson [66,67]; indeed, the DAs of B mesons enter 
the QCD factorisation framework of radiative and non-leptonic B decays, as discussed 
in Chapter El and a variant light-cone sum rule method devised in Ref. [68]. There also 
exist DAs of the photon which describe its "soft" hadronic components, along with the 
usual "hard" electromagnetic components [69]. The photonic DAs can be important in, 
for example, B — > V'-f decays [70] as investigated in Chapter [71 and B — > ^ev [71,72]. 
Finally, the field of baryon DAs is also active and many of the tools and concepts we cover 
in this thesis find application there, see for example Ref. [73] for a review. 

5 The expansion is facilitated by using light-cone coordinates which are given in Appendix \Xi 




(2.4) 
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2.2 The Conformal Expansion 



The standard determination of meson DAs proceeds by making use of the conformal 
symmetry of massless QCD at tree-level. The conformal expansion is analogous to the 
partial wave expansion of wave functions in quantum mechanics in spherical harmonics 
ij}(r,0,(f)) — > R^Ylnui^L^^)- The expansion uncovers a simple multiplicative renor- 
malisation at leading-order, and as such different partial waves, with different conformal 
spin, do not mix under a change of renormalisation scale. At next-to-leading-order this 
is not the case, because strictly speaking the conformal symmetry of a quantum theory 
requires its (3 function to vanish. Proximity to the conformal limit in QCD is therefore 
governed by the value of the strong coupling constant, becoming true as a s — > and we 
pass to the free theory^] Using the QCD equations of motion we can elucidate this mixing 
order-by-order in the conformal expansion. 

The application of conformal symmetry to exclusive processes has recieved a lot of atten- 
tion in the literature, see Refs. [74-78]. The main benefit of the conformal expansion is the 
systematic separation of the longitudinal and transverse degrees of freedom in meson DAs. 
The former correspond to the longitudinal momentum fractions and is given by irreducible 
representations of the relevant symmetry group, SL(2,R). The latter are integrated out to 
yield a renormalisation-scale dependence of the DAs, described by renormalisation-group 
equations. Here we focus on the most important points, see Ref. [79] for a detailed review. 

2.2.1 Conformal Group 

The conformal group is defined as all transformations that change only the scale of the 
metric and as such preserve angles and leave the light-cone invariant g'^ u (x') = u(x)g^ u (x); 
the spacetime interval ds 2 = g^(x) dx fl dx u is conserved up to scaling. These transfor- 
mations form a generalisation of the Poincare group. The full conformal algebra in 4 
dimensions includes fifteen generators 

P M — > 4 Translations, 

M M!/ — > 6 Lorentz rotations, 

D — > 1 Dilatation, 

K M — > 4 Special conformal translations. (2.5) 

6 It must be noted that mass terms break the conformal expansion immediately at the classical level. 
This does not upset the conformal expansion, however. See Ref. [79] for details. 
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Our hadronic picture is of partons moving collinearly in, say the p M direction, existing 
near the light-cone. We therefore restrict the fundamental fields of the conformal group 
to the light-cone <3>(x) — > $(az), where a is a real number, and we assume fields to be 
eigenstates of the spin operator 

= % -<j^ , (2.6) 

so as to have a fixed Lorentz-spin projection s in the ("plus") direction S + _$(«2;) = 
s$(az). For leading-twist operators this is automatically satisfied and for higher-twist 
operators projections are used to separate different spin states, as we shall discuss shortly. 
The full conformal symmetry (12.51) is now modified and it turns out that the resulting 
group of transformations form the special linear group SL(2,R), or so-called collinear 
conformal group, given by just four generators. They are written in standard form by 
constructing the following linear combinations 

L + = Li + iL 2 = -^P+ , L_ = Li — iL 2 = -K_ , 

L = ^(D + M + _), E = ^(D-M+_). (2.7) 

which leads to the familiar relations 

[L , L T ] = =fL t , [L_, L + ] = -2L . (2.8) 

The operators act on the fundamental fields as 

\L+,$(az)] = -<9 a $(cm), (2.9) 

[L_,<$>(az)) = (a 2 d a + 2ja)$(an) , (2.10) 

[Lo, = {ad a + j)$(an) , (2.11) 

[E,$M] = i(Z-s)$(an), (2.12) 

where t = I — s is the £u>zsijll I is the canonical mass dimension^ s the Lorentz-spin 
projection, and j = \{l + s) the conformal spin of the field $. The conformal spin 
specifies the representation of the collinear conformal group. The operator E commutes 
with all Lj and therefore twist is a good quantum number for each conformal field. The 



strictly it is the collinear twist which is defined as "dimension minus spin projection on the positive 
direction" . There also exists geometric twist which is defined as "dimension minus spin" . 
8 For example, / = 3/2 for quarks and I = 2 for gluons. 
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Casimir operator commutes with all Lj and is given by 

[Li, *(<*z)]] = j(j - l)$(az) = L 2 <l>(az) . (2.13) 

i=0,l,2 

At the origin of the light-cone a = and the field $(0) is killed by the lowering operator 
L_ and as such has the minimum spin projection j min of states of conformal spin j. One 
can define a conformal operator O n = 5>(0) by requiring that it transforms just as the 
fundamental field, Eqs. (12.101 - I2.12p . and is killed by the lowering operator L_. The 
raising operator L + can be repeatedly applied to $(0) to give 

O n , n+k = [L + ,..., [L + , [L+ , $(0)]]] = (-id + ) k ® nj (2.14) 

V v ' 

A' 

where O n>n = O n and the subscript n defines the conformal tower of states, of conformal 
spin jmin < j m i n + k < oo, generated by the collinear conformal algebra. This is an infinite 
dimensional representation of the collinear conformal group. 



2.2.2 States of Definite Spin 

Now the main language of the collinear conformal group is defined, if one can relate the 
fundamental fields $ to the operators of hard processes in QCD one can export all the 
machinery above and immediately reap the benefits. To this end, consider the non-local 
two-particle operator at light-like separation (12. ip and expand at small distances 

q(z)Ts(-z) = ^ 1 g(0)(B -^r S (0) , (2.15) 

k 

where D^—D^ — D hl - The question is, how does one express these local operators in 
terms of conformal operators and thus separate all the different twist contributions? To 
proceed one decomposes the quark fields into definite Lorentz-spin components using the 
projection operators 

n + = i^, n_ = i^, n + + n_ = i. (2.16) 

2p ■ z 2p ■ z 
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which project onto the "plus" and "minus" components of the spinor respectively. Using 
the generator of the spin rotations of a spinor field (12. 6p one can show 





(2.17) 



The composite operators ( 12. 15ft have conformal spin projection j = j q + j s + k where the 
subscripts correspond to the separate quark fields. The composite operators are ordered 
by increasing twist 

q+Ts + , g + rs„+g„rs+ , . (2.18) 

t=2 



t=3 



t=4 



It can be shown that the corresponding local conformal operators are 



r 2 ( 

nt=4 



X) 



X) 

4 (x) 



(id + r [fah+Cft 2 [D+ ld + ) s(x) 
(id + ) n [q(x) 1+1±1 .P^ (p+ /d+) s(x) 
(id + ) n g» 7 -Cy 2 (d+ /d+) s(x) 



(2.19) 
(2.20) 
(2.21) 



where d+ =D+ + D+, C™(x) are G eg enbauer polynomials and Pn' s \x) are Jacobi polyno- 
mials. There is another twist-3 operator corresponding to Eq. ( I2.20j) with the replacement 
P^' 0) (x)^P^ 1 \x). One can now connect firmly to QCD with the specific example of 
the leading-twist operator; consider again Eq. (12.11) and specify the twist- 2 Dirac matrix 
r — > (projected onto z^) and define the DA as 

(0\q(z) lz s(-z)\K*(p,\)) = flm K ,(e^ ■ z) [ due^ az+ <-^4 ;K *(u, fx) . (2.22) 

Jo 

Then, using Eq. ( 12.151) and comparing the result to Eq. ( 12.19ft . one finds 

duCT{0<p{ K ^^) = m t u 2 )), (2-23) 

where we introduce the shorthand C, = u — u = 2u — 1, and ((Q^ =2 )) are the reduced 
matrix elements of the operator Q^ =2 . The Gegenbauer polynomials form a complete set 
of orthogonal functions over the weight function 6m (1 — u) in the interval < u < 1 



1 du uuCl'%) CT{0 = <W ( " + 1)( " + 2) 
o 4{2n + 3) 



(2.24) 
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and so one can invert (12.231) to find 



J 2;K 



. (u, /i) = 6«(1 - u) \l + C 3 n /2 (0 \ , (2.25) 



n=l 



where the Gegenbauer coefficients a„,(/i) are related to the reduced matrix elements, 
Eq. (TZ23D, as 

'^>-wfw) mT2)) ' oJ = 1 ' (2 ' 26) 

The result is that the conformal symmetry has separated the longitudinal degrees of 
freedom - as contained in the orthogonal Gegenbauer polynomials which are function of 
the momentum fraction u for the twist-2 distribution - from the transverse degrees of 
freedom, which now show up as the renormalisation-scale dependence of the Gegenbauer 
coefficients at. The Gegenbauer coefficients contain the non-perturbative information of 
the DA and, for the leading- twist DA, are of conformal spin n + 2. The higher- twist 
two-particle DAs are expanded analogously in Pn and Pn' 1 ^ for twist-3 and ClJ 2 for 
twist-4, see Eqs. (12.201) and (12.211) respectively. The explicit expression for the DA of an 
m-particle state with the lowest possible conformal spin j — ji + ■ ■ ■ + j m , the so-called 
asymptotic distribution amplitude, is given by 

if \ r(2j! + ■ • • + 2j m ) 27'i-l 2j' 2 -l 2?™-l fr, r>rr\ 

<f) as {a 1 ,a 2 ,- ■■ ,a m ) = - - - a/ 1 a/ 2 . . . a^ m , (2.27) 

where J^^Lx ctk = 1 [77] . For the twist-2 two-particle DA considered j = 1 and we recover 
the weight function Qu(l — u). Analogously, for a multi-particle DA, states higher in con- 
formal spin are multiplied by polynomials orthogonal over the weight function Eq. (I2.27p . 
The matrix element Eq. (I2.22p with the chiral-odd Dirac matrix V — > o~ zp starts at twist-2 
also, and gives rise to the second two-particle twist-2 DA 

0^K / i) = 6 M (l- M )|l + fja^( / i)C^ 2 (o| • (2.28) 

The expansion of DAs in terms of an infinite sum of partial waves, as in Eq. (12.281) . is 
very general, and is valid at the level of operators. In practice the concept of G-parity 
allows one to classify which Gegenbauer coefficients contribute for specific matrix elements 
of those operators. The G-parity operator Q is defined as Q = Ce~ mT2 where T 2 is the 
isospin generator of the 2 axis and C is the charge conjugation operator. G-parity is the 
generalisation of charge conjugation to particle multiplets, for example, Q\tt ± ' q ) = — 1^'°) 
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and is conserved in QCD. It effectively swaps quarks for anti-quarks and therefore for equal 
mass quarks u <-» u and consequently, for ir, p, to and 0, the odd Gegenbauer coefficients 
a2„+i vanish. 



2.3 Two-Particle Twist-2 Distribution Amplitudes 

As mentioned in the last section, there are two two-particle matrix elements that begin 
at twist-2 [55] @ 

(0\q(x)^s(-x)\K*(P,\)) = 




All other DAs in the above relations are of twist-3 or -4 and all terms in the light- 
cone expansion of twist-5 and higher are neglected. The twist-4 DAs are shown for 
completeness. The normalisation of all DAs is given by 

f du<p(u) = l. (2.31) 
Jo 

The conformal expansions of the leading-twist DAs 4>\. K » and are gi ven by Eqs. (12.251) 
and (12.281) respectively. In the local limit — > the matrix elements (I2.29P and (12.301) 

9 The vacuum-vector meson matrix elements vanish for T = ij^ because it is impossible to construct a 
pseudoscalar quantity from the three available 4- vectors p M , z M and e^ A \ 
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reduce to the longitudinal fy, and transverse decay constants; 

(0\q(0)^s{0)\K*{P,X)) = e^m K *fl, (2.32) 
(0|g(0)^s(0)|iT(P,A)) = i(eMp v -ePPjfa{ii). (2.33) 

Note that f^* (aO is scale-dependent because the tensor current is not conserved. Numer- 
ical values of the decay constants are discussed in Chapter HI The above DAs are related 
to those defined in Refs. [55, 56] by 

02;K* = 0||(JL) > = ^ 7 ^\k* = 93 , 

= A (T) , <f>i, K * = , = h 3 . (2.34) 



2.4 Two-Particle Twist-3 Distribution Amplitudes 

The two-particle twist-3 DAs ^'k* have already been defined in Eqs. (I2.29H and (I2.30p . 
There are two more two-particle DAs, V^il* > defined as ^ 

(0|^) 7m 75s(-^)|^*(P, A)) = i fj c .m K *e i r ft ti x) PaZp Jjue^^u) , (2.35) 
(0|?(z) S (-z)|iT(P,A)> = -ifkM X) -z)™ 2 K * [ due^ z 4. K *(u). (2.36) 
The normalisation is given by 

/ du ^&i( u ) = 1 -7i7ii)— i > ( 2 - 37 ) 

JO /if* x * 

which differs from Ref. [55], where all DAs were normalised to 1; here we keep the full 
dependence on the quark masses. 



10 In the notations of Ref. [55], ip^ K , = {1 - (fk*/.fl*)( m s + m q )/m K *}g { l\ ^| ;if , = {1 - 
(/i//ir.)K + m g )/mW/4 s) . 
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2.5 Three-Particle Twist-3 Distribution Amplitudes 



There are also three three-particle DAs of twist-3: 

(0\q(z)g s Gp z (vz) lzl5 s(-z)\K*(P,\)} = f\ .m K .(p ■ z) 2 e%$\. K {v,p ■ z) + . . . , (2.38) 

(0\q(z)g a G Px (vz)i>y z s(-z)\K*(P,X)) = A*™K*(p ■ zfe^\. K (v,p ■ z) + . . . , 
(0\q(z)g s G z p(vz)a z ps(-z)\K*(P,\)) = /^m|.( e W • z){p ■ z)^. K .(v,p ■ z) + . . . , 

where the dots denote terms of higher twist and we use the short-hand notation 

F(v,p-z) = J Vae- ip - z{a2 - ai+vai) F{a) (2.39) 

with J 7 (a) being a three-particle DA and the integration measure Da is defined as 

J Va = J daida 2 da 3 6(1 — ai) . (2.40) 

The twist-3 three-particle DAs correspond to the light-cone projection r y z G z ± and a± z G± z , 
respectively, which picks up the s — | component of the quark fields and the s = 1 com- 
ponent of the gluonic field strength tensor. According to Eq. (I2.27p . the (normalised) 
asymptotic DA is then given by 360 aic^ct^. To NLO in the conformal expansion, each 
three-particle twist-3 DA involves three hadronic parameters, which we label in the fol- 
lowing way: (, k are LO and u, A NLO parameters. ( and u are G-parity conserving, 
whereas k and A violate G-parity and hence vanish for mesons with quarks of equal mass, 
i.e. p, to and (p. We then have 



360aia 2 «3 | «3ftr» + w\ K . («i - a 2 ) + A| K . ^ (7a 3 - 3) | 



®\. K *(°L) = 360aia 2 «3 jcL» + A L*( a i - «a) + ^\ K *\ ( 7 "s - 3) 
$£jc.(qO = 360a 1 a 2 a3{4jf.+W3j i :.(ai-a 3 )+A^.^(7a3-3)| . (2.41) 

The parameters defined above are related to those of Ref. [55] by: 

C 3 4 = Cl, C 3 y = ^/i4, 

C 3 T = Loi/U, Cl <o = 4 • (2-42) 
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G-parity breaking terms were not considered in Ref. [55]. For equal mass quarks, $ 3 .#* 
are antisymmetric under oi\ <-» «2, whereas $3.^* is symmetric. All these parameters can 
be defined in terms of matrix elements of local twist-3 operators. For chiral-odd operators, 
for instance, one has 

(0\qa^g s G^s\K*(P,X)) = f^m 2 K , ( e W ■ z)(p ■ z)k£ k . , 
(0|g<7*[i£>„ g a G*} 8 - | z^g^C^s] 1T(P, A)) = /^.m^. (e< A > • z)(p • zf A A ^ , 

(0\qiD z a z( g s G z ^s - qa z ^g s G z ^D z s\K*(P, A)) = /^»m^»(e (A) • z)(p • z) 2 — ^3K* I 

the formulas for chiral-even operators are analogous. In Chapter H] we calculate numerical 
values for all the parameters in Eq. (12.411) from QCD sum rules. 



2.6 Relations Between Distribution Amplitudes 

The QCD EOM are a crucial ingredient in simplifying the kinematic contributions of 
different operators, a task which is facilitated by the fact that they are preserved to 
all orders in the conformal expansion. The EOM relate via integral equations the two- 
particle twist-3 DAs, defined in Sections 12.31 and I2.4[ to the two-particle twist-2 DAs, 
Eqs. (12.251) and (12.281) . and three-particle twist-3 DAs, defined in Section 12.51 We do 
not quote the EOM themselves for which we refer the reader to the literature. The 
framework for the procedure was developed in Ref. [77] based on deriving the EOM for 
non-local light-ray operators [63]. The operator relations are then sandwiched between 
the vacuum and meson states and the definitions of the DAs used to convert them into 
integral equations, making use of partial integration to remove explicit dependence on 
co-ordinates and momentum 4-vectors. The resulting expressions are then solved order- 
by-order in the conformal expansion, see Ref. [79] for an overview. The EOM contain 
mass dependent contributions oc m s ±m q that were calculated in Ref. [55]. In the present 
analysis, G-parity breaking terms of the three-particle twist-3 DAs are included, which 
then, via the EOM, impact on the two-particle DAs [40]. The resulting integral equations 
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arc: 
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T(«) 



and 



yIk*( u ) — u / ^rTW+« / du - T(u) 



A K *(u) 



5< 
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with 



<t>i K *( U ) 



it / dv - Q(v) + u I dv — Q(v) , 
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J 



+ 2 
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da>i / d«2 — a± — h a 2 
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Oi3 



d«i 
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dai 



d 
da 2 
d 

d«2 



(2.46) 



Using Eq. (I2.4ip . and the corresponding relations for twist- 2 DAs, one obtains expressions 
for the twist-3 two-particle DAs, which are valid to NLO in the conformal expansion. As 
discussed in Ref. [55], the structure of this expansion is complicated by the fact that 
these DAs do not correspond to a fixed Lorentz-spin projection s of the quark fields. The 
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resulting expansion is in C^ 2 (£) for ip^-K* an< ^ ^n 2 (0 for ^t'tc*'- 



15 



, 4a'* - 4 Kk* ) m 2 - 3) + - c^. (3 - 30^ + 35^ 



3 m s + m q f r + g ^ a j + 3 ( 7 _ 30mi )al + £ In u(l + 3aJ + 6aij) 
-£]nu(l - 3a) + 6a|)| 



2 m K * / 



+ lnw(l - 3a| + 6a5)| , (2.47) 



+ 3 ms + mg \uu{\ + 2£a] + 3(7 - 5uu)al) + uln u(l + 3aj + 6o|) 
+u lnw(l — 3a| + 6a|)| 

-3 ms ~ mq |nn(9a! + 10fa|) + ulnu(l + 3aJ + 6a|) 

™>K* J K * L 



-«lnu(l -3aj + 6a!J)} , (2.48) 

+ go! + ^cL + ^4.* - |sL.) ^ 2 3/2 (0 + QaL* - ^aL.) ^ /2 (o} 

^2 —I— 777, 

+ 6 — q - ^f- {uu{2 + 3£af + 2(11 - 10tm)a£) + u In u(l + 3af + 6a±) 



+u In u(l — 3a± + 6a^") } 



{uu(9a^ + 10£<4) + u lnw(l + 3a] 1 + Qa^ 



-«lnu(l -3a^ + 64)} , (2.49) 
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15 \ll 

)K 3K* A 3K* 



~Mu) = ^(l + a + ^M + {^4 + 5CL}(3e 2 -l) + {5. 

15 ~ii 1 f 9 ii 15 ii 15 ii 1 

+y xIk* j £(5£ 2 - 3) + | — 4 + - o4. - - 4 K * | (35£ 4 - 30£ 2 + 3) 



+ - ms + m « 1*1 {2 + 9^ + 2(11 - 30uu) 
2 m K * fl , 



a 2 



+ (1 - 3a] 1 + 6aa ) lnw + (1 + 3a] 1 + 6a^) lnw} 

q - J -^f- |2£ + 9(1 - 2uu)a± + 2£ (11 - 20uu)a^ 

2 m K * /|* 

+ (1 + 3af + 6a^) lnu - (1 - 3a^ + 6a^) lnw} . (2.50) 

These expressions supersede those given in Ref. [55] where G-parity violating terms in 
«3 and A3 were not included. The DAs given above now contain a minimum number 
of parameters which can be determined from one's favourite method, such as QCD sum 
rules or Lattice QCD. We briefly mention that the matrix elements of QCD operator 
identities can also be used to relate twist- 2 and twist-4 DA parameters to each other; such 
investigations were performed for the G-partiy violating twist-2 K* parameters ^[^(K*) 
in Refs. [60,62]. 



2.7 Evolution Equations 

The scale dependence of the leading- twist DAs of Eqs. (12.251) and (12.281) can be investigated 
using perturbation theory. The resulting renormalisation-group equation is the Efremov- 
Radyushkin-Brodsky-Lepage (ER-BL) evolution equation [43-45,47] 

d f 1 

^ 2 ^2 0(u,^) = y dvV(u,v;a s (fi))<i)(v,fi) , (2.51) 
which completely specifies (j)(u, /S) given <fi(u,[i ). The kernel is given by an expansion in 



V(u, v; o.(jjl)) = (u, v) + 2 V M(u, «) + .... (2.52) 

The evolution equation (12.511) can be solved readily at leading-order using the conformal 
expansion [74,76,78,80]. This amounts to finding its eigenfunctions, which we already 
know to be Gegenbauer polynomials, and using this fact it can be shown that the leading- 
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order kernel can be written as 

VW(u,v) = -6uuJ2 (w + n ( t + 2) 1 - C - 2{2U ~ 1)CI ' 2{2V - l) ' (2 ' 53) 

71=0 

giving the LO anomalous dimensions of the Gegenbauer coefficients j n . The renormalisa- 
tion is hence multiplicative at leading-order 

a L n °(^)=aM)L^ /m K (2.54) 

where L = a s (fj, 2 )/a s (nl). The one-loop anomalous dimensions of the twist-2 Gegenbauer 
coefficients are [23,81] 

- ^(gM- 2(n+1 ; (n+2) ). ("5) 

■C = 8C '(§H)- (2 - 56) 

At next-to-leading-order the scale dependence is more complicated [82,83] 

n-2 

^ LC V) = «.(^)^ NLO + ^ £ , (2.57) 



47T 

fc=0 



where 



7n are the diagonal two-loop anomalous dimensions which are available for the vector 
current [84], and the tensor current [85]. The mixing coefficients d^l are given in closed 
form in Refs. [82,83] where the formulas are valid for arbitrary currents by substitution of 
the corresponding one-loop anomalous dimension. For the lowest moments n = {0, 1, 2} 
one has 

lim lid) 23110 512 Ar iim 34072 830 XT 

7! (1, = 0, 7 }« = _-_ W/ , 7 |« = _-_ W/ , (2.59) 

± m 724 104 ±m , AT 38044 904 Ar 

7o ( } = -g" - 27 ^ ' 71 = 124 " 8Nf ' 72 = "W ~JT Nf > (2 ' 60) 
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and 




35 20 - 3/3, 
¥ 50 - 9/3, 



(1 _ £50/(9/3o)-l 



) , d 



,±(i) 
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28 16 - 3/3, 
¥ 40 - 9/3, 



1 _ ^40/(9/3 )-l 



) • (2-61) 



It is evident from Eqs. (12.551) and (12.561) that the anomalous dimensions of the Gegenbauer 
coefficients increase logarithmically with conformal spin. This implies that as \i — > oo the 
coefficients higher in conformal spin are damped and the DA approaches its asymptotic 
form (123711 



This limit offers a great simplification in that if it can be verified, experimentally of other- 
wise, that a given process is well described by the asymptotic form of the DA at hadronic 
scales, there are no non-perturbative parameters to be determined. The convergence of 
the conformal expansion in general has to be verified case by case and there is no a priori 
reason why it should do so. In practice one has to truncate the expansion at some order in 
conformal spin, usually n = 2, and as this constitutes an approximation it thus introduces 
a model dependent assumption. In Ref. [86] for example, these issues are discussed and 
an alternative method suggested. 

The scale dependence of the three-particle twist-3 DAs can in principle be deduced from 
evolution equations for the non-local operators in (12.391) using the techniques of Ref. [63]. 
The evolution of the parameters in ( 12.411) could then be found by projecting out the 
desired conformal spin, as for the leading-twist DA. Another approach is to consult the 
literature of results for the corresponding nucleon structure functions [87]. The three- 
particle twist-3 parameters Qk*> ^k* an d ^k* renormalise multiplicatively in the 
chiral limit, and the others mix with each other. For non-zero strange quark mass, there is 
additional mixing with twist-2 parameters with the mass corrections featuring as m s ± m q 
depending on the G-parity of the parameter. Here, we write down explicitly only the 
renormalisation-group improved relations for the above 5 parameters. The relations can 
be written in compact form as 




(2.62) 



[l—fOC 



3 



P,(/i 2 ) = L^P° PM) + °ij (£ (7Q) - //3 ° - L^P°) QM) , (2.63) 
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where the parameters are given by: 



P 

Ql(2) 

7p 
(To) 1,2 



{fl<*C?,K*i /i^* K 3A'*' fk* K ZK*^ fK* u 3K*y flC*^3K*} 

— — {m s ± m q , (m s =F m q )a^, (m s ± mja^ } , 

/"-* II II 

— — {m s - m 9 , (m s + m q )a\, (m s - m q )a\} , 

m K * 

{m s + m q , (m s - m q )a\, (m s + m 9 )a|} , 
rriK* 

77 77 55 73 104 



fl 



9 ' 9 ' 9 ' 9 ' 9 



16 



68 86 



(tq)3,4,5 - { 4, — , , , 



(2.64) 
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(2.65) 



We refrain from delving into a full discussion here of the mixing including the remaining 
parameters A|, cj|, and guide the reader to Appendix A of Ref. [40] for details. 
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Chapter 3 
QCD Sum Rules 



The original QCD sum rule approach was introduced by the revolutionary work of Shif- 
man, Vainshtein and Zakharov in Refs. [13-15], and has proven itself to be one of the most 
effective tools for determining non-perturbative parameters of low-lying hadronic states. 
It does so in terms of a finite number of universal non-perturbative input parameters, 
and as such has great predictive power. The approach has been massively successful in 
ascertaining a wide range of phenomena of non-perturbative origin. QCD sum rules are 
particularly advantageous for B physics because the presence of an intrinsic heavy-quark 
mass scale mj provides the necessary conditions required for the application of the short- 
distance OPE or light-cone expansion of relevant correlation functions, from which the 
relevant quantities can be extracted. The heavy-quark limit mj — > oo is not necessary 
and sum rules can be derived in full QCD for finite mj. Despite its successes, the method 
is limited by an inherent irreducible systematic uncertainty of 20 — 30%. However, such 
is the relative ease of the QCD sum rule method, as compared to, for example, Lattice 
QCD, that its place in the tool-box of the QCD practitioner is ensured. 

Firstly we discuss step-by-step the methodology of the original QCD sum rule approach. 
Secondly we discuss its modification to accommodate non-local correlation functions which 
aides the extraction of DA parameters of beyond leading-order in conformal spin. Thirdly, 
we outline an extension of the original approach; light-cone sum rules. All three methods 
find application in this thesis: 

• in Chapter H] we make use of the non-local formalism to extract numerical values 
for the leading twist- 2 and twist-3 DA parameters defined in the Chapter [21 
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• In Chapter [5] we calculate important contributions to the semileptonic B — > 77 ^ 
transition form factors in light-cone sum rules. 

• An example of the original method is presented in the last section of this chapter, 
section I3.3[ where we calculate the a s corrections to the gluon condensate contribu- 
tion from a local correlation function. 

We focus solely on the points required for future chapters. For more information on sum 
rules, see for example Refs. [88-91]. 

3.1 SVZ Sum Rules 

The original sum rule method, which we refer to as SVZ sum rules, parameterises un- 
known non-perturbative QCD vacuum effects in terms of the so-called universal vacuum 
condensates. These quantities are vacuum expectation values of local operators Oi that 
vanish in perturbation theory by definition and are ordered by their dimension D. 

The calculation of a QCD sum rule starts from the calculation in QCD of a suitable corre- 
lation function in which the mesons are represented by interpolating currents possessing 
the correct quantum numbers. The method proceeds by equating two different represen- 
tations of the correlation function. The first is obtained by performing a short-distance 
OPE, the result of which is matched to a second representation, in terms of a dispersion 
relation over physical hadronic states, leading to a sum rule from which various properties 
of the hadronic states can be extracted. 

The SVZ sum rules find an important application in determining the universal hadronic 
parameters that appear in meson DAs. Indeed, some of the first SVZ sum rule calculations 
were performed to extract decay constants f n>p [15], and Gegenbauer coefficients a n [94-96] 
of light meson DAs. The method can only be applied to parameters of the lowest few 
orders in conformal spin; parameters higher in conformal spin must be determined from 
other methods because the sum rules become unreliable. One such method is that of 
non-local condensates, see for example Ref. [97]. 
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3.1.1 Correlator 



The following two-point correlation function describes the propagation of a quark-antiquark 
pair 



where q is the incoming momentum and possible Lorentz indices are omitted for simplicity. 
The local currents Jj are chosen to have the correct quantum numbers and particle content 
corresponding to the particular hadronic parameters under investigation. The physical 
picture of a hadron is of quarks, antiquarks and gluons confined within a typical hadronic 
size R which is large when compared to the scale associated with perturbative effects. 
If one can show, however, that the correlation function is dominated by small spacial 
distances and time intervals 



for a certain momentum configuration, then one has ensured the small size of the strong 
coupling a s and hence the use of perturbation theory in our calculation. One begins by 
noting that after contracting any Lorentz indices which may appear in the currents J, 
in Eq. (13.1 1) , the correlation function can only depend on the interval x 2 = x\ — x 2 . By 
taking the Fourier transform, completing the square, and shifting the variable x one finds 



The integral is dominated by the region where the arguments of the exponential vary 
slowly. This condition requires k ~ 1/x 2 and k ~ Q 2 which are both fulfilled for x 2 ~ 
1/Q 2 ; for large momentum transfer the quarks propagate near the light-cone 3 To find 
the true short-distance dominance one needs to dig a little further and by choosing the 
Lorentz frame qo = one finds x 2 ~ l/Q 2 as required (13. 2h . In the case of light quarks 
one needs the momentum transfer to the quarks to be large, Q 2 = —q 2 ^> Aq CD . In 
the case of heavy quarks a large energy scale is introduced through the quark mass, 
for example mf,, which then serves to set the characteristic distances for the correlation 
function |x| ~ Xq ~ l/(2mf,); one is thus automatically in the perturbative regime. 

1 An expansion round x 2 — > is the basis of QCD sum rules on the light-cone - see Section f3T2l 






(3.2) 




(3.3) 
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3.1.2 Short-Distance OPE 



The first of the two representations of the correlation function is obtained by performing 
the QCD calculation, valid for Q 2 = — q 2 ^> Aq CD , using the short-distance OPE 



where the non-perturbative long distance effects of QCD are encoded in the condensates 
(Oj) and the short-distance effects are included in the Wilson coefficient functions Cj which 
are calculable in perturbation theory. Both the condensates and their coefficients are in 
general renormalisation scale dependent. Perturbative corrections to the condensates are 
calculated when necessary. The perturbation theory contribution to Eq. (13 .4p has D = 
and corresponds to the unit operator (Opt) = 1. The condensates play the role of power- 
corrections and are suppressed by inverse powers of the hard scale as (Q 2 )~ D / 2 . In the 
asymptotic limit Q 2 — > oo only the unit operator survives, corresponding to asymptotic 
freedom. 

3.1.3 Condensates 

The condensates represent the effects of non-perturbative QCD and they cannot be de- 
termined from first principles due to the unknown nature of the QCD vacuum. The 
determination of the condensates is an industry in itself. The light quark condensate 
(0|gg|0) has been known for a long time [98] and it drives the breakdown of the chiral 
symmetry of the light quarks q = {u, d} and its value can be extracted from experiment: 



where we use the notation (Oj) = (0|Oj|0). To define other condensates, one notes that 
the only vacuum expectation values of operators that can survive are those which are 
Lorentz invariant, spin zero, colour and flavour-singlets i.e. possess the quantum numbers 
of the vacuum. The complete set of condensates (Oj) that contribute with D < 6 are 




(3.4) 




(m u + m d ) (qq) , 



(3.5) 



m q (qag s Gq) , 



D=6 



D=0 



(1) , 
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where q = {u,d, s} is a light quark spinor and all indices are contracted]^ We assume 
isospin symmetry for q = {u, d} and one must differentiate q = s when SU(3)p-breaking 
effects are taken into account. Higher dimensional condensates D > 6 are not very 
well determined and generally unknown. If required, however, they can be estimated 
by employing the vacuum saturation hypothesis whereby the operator fields are simply 
split to form products of known condensates; for example, the quark- ant iquark D = 6 
operator can be simplified to the product of two qq operators [14,99]. In practice, the 
OPE is truncated to a given order, and is usually justified by the stability of the resulting 
sum rule. The series, Eq. (I3.4j) . is then given in terms of a limited number of condensates 
allowing sum rules to be written in terms of a small set of parameters incorporating the 
general features of non-perturbative QCD, while retaining its predictive power. 

The procedure works in reverse, of course, where the values of condensates are deduced 
from sum rules for which the hadronic parameters are known from other methods; two- 
point correlation functions featuring irf^b or c , ~f fl c currents correspond to the T and J/^J 
resonances respectively, of which the decay constants and masses are known. Values for the 
condensates are given, along with other input parameters, in Appendix [B] Uncertainties in 
the values of the condensates and other input parameters constitute part of the reducible 
theoretical uncertainty of the sum rule approach. 

3.1.4 Dispersion Relation 

To proceed we need to relate the result of the OPE to a second representation of the 
correlation function which is obtained in terms of the spectrum of hadronic states in the 
physical region q 2 > 0. This is done via a dispersion relation, which is derived from the 
analytic properties of the correlation function as follows. The function II(g 2 ) is analytic 
in all q 2 except on the real axis starting at a pole corresponding to the ground state 
particle. At higher energy higher mass excited states and a continuum of many-particle 
states also feature. The higher mass resonances give poles above the ground state, the 
details of which depend on the physical spectrum of particles which possess the correct 
quantum numbers to couple to II. The continuum of many-particle states, correspond to 
a continuous cut, see Fig. 13.11 

2 The heavy quarks c, b and t do not form condensates because they are too massive to interact 
non-perturbatively with the QCD vacuum. 
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Figure 3.1: The general features of a spectral density function p (s) in the complex plane. 
The blob represents the pole due to the ground-state, the cross possible poles due to higher mass 
resonances, and the thick line the cut due to the continuum of multi-particle states. The dotted 
line is the integration contour. 



Using Cauchy's formula we can write 

rw 2x 1 / , n(z) 1 f R , n(z + ie)-II(*-ie) , N 

|*|=A 

where the region of integration is split into the parts just above and below the positive 
real axis and the circle of radius R. Provided that the correlation function vanishes at 
least as quickly as q~ 2 as \q 2 \ ~ R — > oo then the integral over the circle at radius R 
goes to zerolfl The remaining integral can be simplified using the fact that below the first 
pole at q 2 = s m i n , Il(g 2 ) is real and above this point, according to the Schwarz reflection 
principle, 11(2: + ie) — Ii{z — ie) = 22lmll(2; + ie). Hence 

I% 2 )= r ds P( 2 S) , , (3.7) 
/, . s — q l — ie 

where the function p(s) = ilmll(s) is the spectral density and describes the physical 
particle spectrum as a function of energy s. 



3 If II does not vanish quickly enough we subtract the first few terms in its Taylor expansion as required. 
We shall see that this does not matter in the end, due to the Borel transformation. 
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3.1.5 Unitarity Relation 



As we have seen, for large negative q 2 our correlation function is dominated by short- 
distance physics. As q 2 becomes more positive the separation of the quarks increases. For 
large enough positive values of q 2 long-distance QCD interactions become more important 
and the correlation function then describes the creation of hadrons, which is the basis of 
its second representation. As discussed in the last section, II uncovers a very complicated 
spectrum of states for q 2 > 0. We describe this situation by using the unitarity relation, 
which allows in insertion of a complete set of states into the correlation function 

1 = J2 I dn n \n(p))(n(p)\, (3.8) 

where dQ n includes all phase-space factors and momentum conservation and the sum runs 
over all possible particles and polarisations, starting from the ground state M of mass 
ttim- Inserting (13 .Sp between the currents of our original correlation function (13.1 1) yields 
an expression which we can relate to the hadronic spectral density 

n ha V) = j ( ^ 4 m /_ p2 j d*xe^(0\J 1 (x)\M(p))(M(p)\J 2 (0)\0} + ... , (3.9) 

where the dots denote higher mass states which contribute to the continuum. We are 
usually interested in the ground state, and can insert the expressions for the matrix 
elements on the right hand side. The local matrix elements considered here can be used 
to extract vacuum-meson decay constants, for example. Using the unitarity relation (13.91) 
one can single out the ground state M by comparing it to ( 13. 7\) and writing the hadronic 
spectral density as: 

P had (s) = fu S(s - m 2 M ) + p cont (s), (3.10) 

where fu is directly related to the matrix elements of the currents J\ and J 2 in Eq. (13. 9p . 
For example, one could choose J\ = j\ = cpj z s to extract (f^*) 2 c.f. Eq. (12.321) . The exact 
form of the spectral density beyond the ground state is mostly unknown and the higher 
mass states and continuum contributions are usually lumped together in one function 
p cont (s). If the next highest particle above the ground state occurs at an energy not very 
much higher than tum then it is possible to explicitly include this particle as another delta- 
function term, analogously to the ground state. This procedure was used, for example, 
while investigating the leading-twist K* and p DA parameters for which the relevant 
correlators couple to the K\ and b\ resonances respectively [54,61]. 
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3.1.6 Quark-Hadron Duality 

It is possible to write the result of the OPE as a dispersion relation, with spectral density 
p OPE (s). As p cont (s) is mostly unknown we replace it by p OPE (s) above a certain energy 
so 

p cont (s) ^ p OPE (s)e(s _ So)> (311) 

This assumption relies on the validity of the hadronic representation being approximated 
by the partonic representation at higher energies. Inserting Eqs. (13.101) and (13.111) into 
Eq. (13.71) one finds 

n had (g 2 ) = -J»— + r ds pOP } S \ . (3.12) 
™ M ~ 1 J so s-q 2 -xe 

Now the assumption is not so strict because we only require a duality between the inte- 
grated spectral densities, not the spectral densities themselves. This is called semi-global 
quark-hadron duality. The parameter sq is called the continuum threshold and its value is 
specific for each particle spectrum being roughly equal to the energy of the next highest 
resonance above the ground state: s ~ (ttlm + A) 2 where A ~ C(Aqcd)- Ultimately 
it must be determined from the sum rule itself by requiring the numerical value of the 
determined quantity to be largely insensitive to its variation and this introduces the first 
source of systematic uncertainty to the sum rule method. We are now in a position to 
equate both representations 

n had = n OPE , (3.13) 

to derive our sought after sum rule, however, before we do so, there is one last procedure 
to discuss, which greatly improves the behaviour of the sum rule. 



3.1.7 Borel Transformation And The Sum Rule 

The sum rule can be improved by suppressing the continuum contribution, which we have 
assumed to be well described by p OPE (s > sq) and the possible detrimental impact of this 
assumption is thus reduced. We do this by performing a Borel transformation to both 
sides of the sum rule. The transformation is obtained by applying the operator 

B= lim LAI , (3.14) 



-qi/n=M 2 
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which takes a function of q 2 and produces a new function of the Borel parameter M 2 . 
One frequently encountered example is 



B 



1 




(3.15) 



(m 2 — q 2 ) k 



{k-l)\ (M 2 ) k 



providing an exponential suppression of the unknown continuum contributions, and a 
suppression of the power-corrections by factorials thus reducing the impact of neglected 



to Eq. (13. 7p . which can only appear as polynomials in q 2 , are killed off. The Borel 
transformation improves the stability and accuracy of the sum rule. 

The Borel parameter M 2 is the second and last sum rule specific parameter to be intro- 
duced; along with sq it is required to impact very little, when varied, on the numerical 
value of the quantity being determined. The variation of M 2 changes the relative impact 
of the power-corrections and perturbation theory contributions. In evaluating sum rules 
one looks for a Borel window which is usually in the range 1 GeV 2 ^ M 2 ^ 2 GeV 2 for 
a typical mesonic DA parameter. The sum rule should be reliable if a weak dependence 
(a plateau) is found, the contribution from the continuum is small, and there are no 
unnatural numerical cancellations. 

We now equate Eqs. (13.71) and (13.41) to reach the sum rule 



where the hadronic quantity Jm is given as a function of the universal non-perturbative 
condensates, the perturbative short- distance coefficients as calculated from QCD, and 
the sum rule parameters Sq and M 2 . The sum rule is saturated by the ground state and 
higher mass states are suppressed. As the correlation function ( 13. ip does not depend 
on the renormalisation scale, the /i dependence of the condensates, when multiplied by 
their coefficient functions, must cancel in the sum of (13. 4p . The sum is always truncated, 
however, and the residual /x dependence will be a source of theoretical uncertainty. 



higher dimensional condensates. Also, as B(q 2 ) k = 0, any subtraction terms introduced 




(3.16) 
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3.1.8 Non-local Formalism 



One way to gain access to parameters higher in conformal spin is to calculate sum rules 
involving operators which are related to moments of DAs 

(0|g(0)(S -zfTsmV) ~ /* du (2u - l)fy(u) = (£ fc ) . (3.17) 

Jo 

For the K* for example the first few moments of both the leading-twist DAs are = 1, 
(?) = fai(#*), (O = ^(7 + 12a 2 (K*)) and (£ 3 > = ^(Tla^K*) + 20a z (K*)). A 
more elegant method, enabling the DA parameters to be extracted individually, relies on 
calculating a correlator of two currents, one of which is non-local, with fields at light-like 
separations [z 2 = 0) [59]. Consider the following 

IL(q-z)=i J d 4 xe iq -' x (0\TJ(x)s(0)~f z q{z)\0) , (3.18) 

where J(x) is local, and the non-local current yields the leading- twist DA ( 12. 25ft . The 
sum rule ( I3.16P then reads 

fjfl* e~<> /M2 f du e- taq - z 4- K * = r ds e~ s/M2 [ du e^ z p OPE {s, u) . (3.19) 
Jo Jo Jo 

The integration over u on the right hand side naturally arises via the Feynman parame- 
terisation used in the calculation. At this point one can exploit the orthogonality of the 
Gegenbauer polynomials by replacing the exponential weight function e~ %( - q ' z Cn^iO 
on both sides to project out a)n(K*) via Eqs. (12.231) and (I2.26p . In Fig. 13.21 we show the 




Figure 3.2: A generic non-local diagram. The dotted line denotes the path ordered gauge 
factor [z, —z] between the two quark fields. The momentum q is injected at point y - the vertex 
on the right hand side. 

leading diagram of the non-local correlation function (13.181) . The dotted line denotes the 
path ordered gauge factor [z, —z] between the two quark fields. The non-local formalism 
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allows, in principle, an extraction of parameters of arbitrary order n. In practice, however, 
only the parameters of the lowest few orders n are accessible due to instability of the re- 
sulting sum rules. One finds that the power-corrections in p OPE grow with positive powers 
of n compared to the perturbative contribution. For high enough n this behaviour upsets 
the hierarchy of contributions to the OPE, where non-perturbative terms are expected to 
be moderately sized corrections to the leading term. Hence the method is justified for 
low-order coefficients n < 2 where the non-perturbative coefficients describe the general 
features of the DA. It breaks down for higher-order coefficients n > 3 because the local 
vacuum condensates appear with 5-functions which cannot accommodate the information 
needed to describe the more detailed shape of the DA, see Refs. [59, 102]. 

3.2 QCD Sum Rules On The Light-Cone 

A modification of the QCD sum rule method known as QCD sum rules on the light- 
cone, or light-cone sum rules (LCSRs) [103-105], was developed to overcome difficulties 
encountered when calculating transition and electromagnetic form factors in the SVZ 
method^ The problems are related to the asymptotic scaling behaviour of the form 
factors in the heavy-quark limit m;, — > oo. LCSRs rely on the use of DAs as their 
universal non-perturbative hadronic input and lead to the correct asymptotic behaviour 
in the heavy-quark limit. The DAs represent a partial re-summation of the operators 
appearing in the condensates and appear ordered in contributions of increasing twist [102]. 
We can view LCSRs as a marriage of the SVZ technique and the theory of hard exclusive 
processes [41,43,44,46]. In the case of the "heavy-to-light" B — > M transition form factors, 
LCSRs have been applied successfully to pseudoscalar transition form factors [107-110] 
and vector transition form factors [111,112]. 

For LCSRs to become competitive with the SVZ sum rules, a good knowledge of higher- 
twist DAs is required. This motivates the determination of the non-perturbative DA 
parameters via SVZ sum rules and via LCSR, the DAs themselves to determine other 
non-perturbative parameters, such as transition form factors. As with SVZ sum rules, 
the starting point of LCSRs is with a suitable correlation function. For the extraction of 
B — > M transition form factors we require a two-point correlator, this time sandwiched 
between the vacuum and the meson state M, which is the example considered in this 
section. One employs much of the same methodology as in the last section, although now 

4 The term "light-cone sum rules" first appears in Ref. [106]. 
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one requires the correlation function to be expanded in an OPE on the light-cone. In 
doing so one finds that the correlation function factorises and can be written in terms of 
a convolution of hard scattering kernels and the universal non-perturbative DAs of the 
light-meson. To that end, consider a correlation function of two quark currents taken 
between the vacuum and an on-shell meson M 

n(g,PB) = i j d'xe^iMip^TJ^jUoM , (3.20) 

where j B = m^qi^b is the interpolating current of the B meson which defines the B 
meson decay constant 

(B(p B )\j B \0) = m 2 B f B . (3.21) 

The current J\{x) is chosen to project out the form factor of interest. The momentum 
q is injected into the weak vertex, p B is the momentum of the B meson and p is the 
momentum of M with q + p = p B . The correlation function is dominated by light-like 
distances for virtualities 

m 2 b - p% > 0(A QCD m b ) , m\ - q 2 > 0(A QCD m 6 ) , (3.22) 

which ensures the slow variation of the exponential in Eq. (I3.20p and its suitability for 
an expansion around the light-cone. The light-cone expansion results in the transverse 
and "minus" degrees of freedom being integrated out, leaving the longitudinal momenta 
of the partons as the relevant degrees of freedom. As a result a cutoff /i is introduced 
below which the transverse momenta are included in the resulting light-mesons DAs. The 
contributions from momenta above this cutoff are calculable in perturbation theory. The 
procedure yields the collinear factorisation of the correlation function 

W,Pl) = Yl f duT^^q^pl^)^^^), (3-23) 

n Jo 

where u (1—u) denotes the momentum fraction of the outgoing quark (antiquark) and the 
sum is over all twist and possible polarisation contributions. The scale dependence of the 
hard scattering kernels must cancel that of the DAs 4> n -M- The factorisation formula 
has to be verified by direct calculation and a proof to all orders in a s does not exist. 
The verification relies on the cancelation of divergences, of which there are two types: 
the IR and UV singluarities arising from loop calculations and so-called soft singularities 
which appear when the convolution over u does not converge at the end-point regions 
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[u ~ or 1) i.e. when one of the quarks is soft. In terms of kinematics there are two 
main contributing processes: the hard-scattering mechanism and the soft contribution or 
Feynman mechanism. Both mechanisms are included in the LCSR approach for which 
there are no soft divergences and the IR/UV divergences can be treated in dimensional 
regularisation. 

One can write the result of the light-cone expansion (13.231) as a dispersion relation in p 2 B 

n LC (p|,g 2 )= Hds^p-. (3.24) 

Jrn\ S Pb 

Taking the imaginary part, to obtain p LC , is straight forward after integration over the 
momentum fraction u is performed. The correlation function has a cut in p 2 B starting at 
over the physical region. One now matches this calculation to the hadronic representation 
of the correlation function, which can also be written as a dispersion relation 

n had (p|,g 2 )= rds ^fK (3.25) 

Jm\ s Pb 

where the physical spectral density is given by the ground state B meson plus higher mass 
states forming a continuum as 

p had (s, q 2 ) = F M 6(s - m 2 B ) + p LC (s, q 2 ) Q(s - s ) . (3.26) 

The quantity Fm will contain the form factor we require. We perform the Borel transfor- 
mation to arrive at the LCSR 

F M e - m B' M2 = ds e- s / A V C (s, q 2 ) . (3.27) 

Jrnl 

To extract the form factor we need to find a sets of parameters M 2 and sq such that 
the form factor is largely insensitive to their variation. As with SVZ sum rules, there 
is no rigourous way to do this and so the procedure introduces the irreducible source of 
uncertainty to the method. 
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3.3 Example Calculation - The Gluon Condensate 



Here we present an example calculation within the SVZ sum rule framework. The result 
of the calculation is used in the sum rule for the G-even K meson three-particle twist-3 
DA parameter f 3K , see Ref. [57]. We calculate the a s correction to the gluon condensate 
(c^G 2 ) proceeds from the following local correlation function 

n (G2) = i J dSe^(0\Tm^ z 9sG,Ms(0)KyK z 9sG u M<l(y)\0) , (3.28) 

for which the leading-order contribution vanishes. A convienient way of extracting the 
gluon condensate is to make use of the back- ground field technique in which the fixed- 
point gauge allows the Taylor expansion of quark and gluon fields to be written in a 
gauge-covariant form, see Ref. [113] for details. The gluon field in the QCD Lagrangian 
(11. ip is split into "quantum" and "classical" (background) fields 

- a- + A* , (3.29) 

where the background field A a „ is taken in the fixed-point gauge at xq = 0. The quantum 
field a° is taken to be in the Feynman gauge, thus requiring the gauge fixing term (£ = 1) 

£ fix = -~(<9X + 9sf abc A b ^) 2 , (3.30) 

to be added to the QCD Lagrangian. The quantum field propagates perturbatively and 
we may use the standard expression 

afaTaM = z5 ab J ^D, u (De- lHx ~ y) , D^l) = (3.31) 

The background field does not propagate perturbatively, and is the field that goes to 
form the condensate; it represents the low-energy, long distance modes of the gluon field 
that probe the non-perturbative structure of the QCD vacuum. The fixed-point gauge 
condition allows A a Ax) to be expressed in terms of the gluonic field strength tensor as 

oo 1 

A » = E^TT-^^ 1 -^' 1 [^(°)> [^(°)> [- •••]]] , (3-32) 

n=0 V ' 
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and translating to momentum space one finds 



AW) = -"^(O) [W^W] ^ + ■ ■ ■ . (3-33) 

where we only require the first term; higher order terms give rise to higher dimensional 
condensates which we do not consider. As we have to introduce two condensate gluons 
to construct (G 2 ) we introduce two auxiliary vacuum momenta k and k' for which the 
fixed-point Xq = is a sink. After integration over coordinates these momenta appear 
in the quark and gluon propagators. The two corresponding derivatives are then taken, 
and the vacuum momenta set to zero. The following expression proves very useful in 
managing derivatives of quark propagators 

■£-S®(p) = -S^(p)rS^{p) , S®(p) = 4^ , (3-34) 

OPfj, V m q 

for arbitrary quark flavour q. The gluon condensate is finally extracted using 

G£ M (0)G^(0) = J ab _ - (wv - g w g^) (G 2 ) , (3.35) 

where D is the spacetime dimension. Due to Eq. f l3.29f) the expansion of £qcd yields 
"interaction" terms in which background fields are radiated from the propagating gluons 
at single or double vertices, both of which contribute to the 0(a s ) correction to the gluon 
condensate. These vertices are shown in Fig. 13.31 and the corresponding terms are 

££f = -\g s f abc [(d»A av - dTA a »)ay v 

+ ^a au - d u a a »)(Ay v + a\Al) + 2(d IM a a ")A bl/ a c u ] , 
C^ aa = -\g 2 J abc f ade [A\A d »a%a cv + A\a^A%a cv + A\a^A d v a w \ , (3.36) 

where terms which vanish eventually via Eq. (I3.35P due to f abc S bc = are omitted. 
Contributions also stem directly from the gluonic field strength tensors in Eq. (13. 28ft which 
give rise to gluon emission of either one or two fields from the vertices at co-ordinates 
and y. Due to the gauge condition there is no "left-right" symmetry and all diagrams 
with two gluons, of which at least one is a condensate gluon, emerging from the vertex at 
x = vanish due to ^4^(0) = 0. Diagrams with two condensate gluons at point y, which 
originate from the non-linear part of the gluonic field strength tensor, also give zero due 
to f abc S bc = 0. There is an "up-down" symmetry where diagrams related by a reflection 
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Figure 3.3: The interactions of the background field (denoted by a cross) with the quantum 
field corresponding to Cf^ t a and Cf^ aa respectively. 



in the central horizontal axis are equal. Overall we find there to be 10 distinct non-zero 
diagrams to be calculated which are shown in Fig. 13.41 

Some of the diagrams are divergent, however, all divergences cancel in the sum of all 
diagramsjj For an explicit example consider the last diagram in the second line of Fig. 13.41 
It is evident that we require C^ t a to be contracted in all possible ways with quantum fields 
originating from the linear part of the gluonic field strength tensors at points and y. 
This, multiplied by the condensate field originating from the quark loop yields the gluonic 
part of the calculation 

~ Aj{v) {d,a a M - 0,a»(O)) C tTH (d u a b z (y) - ^(^)) | aUcont _ tioils , (3-37) 

which is eventually given in momentum space by (omitting Lorentz indices) 

d d f(l,k>) 



" dkdk'P(l-kf^ ' (3 - 38) 

where the condensate gluon within C^ t a (w) is expressed by Eq. (13.331) with momentum 
k' and f(l, k') is a function of the loop momentum / and the vacuum momentum k' . The 
quark loop yields a usual trace 

tr 1^ + 4-^^ + ^^*] 

(p + q - i) 2 {p + k) 2 p 2 ' 1 J 

and after multiplying together Eqs. ( 13.381) and ( 13.391) . performing the derivatives in k and 
k' and integrating over the momenta p and I we find 

n cxamplc - ^ - { - G ) — ^ . (3.40) 
In this way we can include all the other diagrams shown in Fig. 13.41 to obtain the contri- 



3 We use dimensional rcgularisation and the M S rcnormalisation scheme throughout this thesis. 
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Figure 3.4: The diagrams contributing to the gluon condensate at G(a s ) for the SVZ sum rule 
of the K twist-3 DA parameter fox ~ see Ref- [57]. For each diagram the fixed-point xq = is 
at the left most vertex and the right most is at y. 



bution to the sum rule 

5184 n \ n I q 2 ' K J 

which differs from the result obtained in Ref. [96]; the logarithmic term is not reproduced: 

-logf^/^CV (3.42) 



/i 2 J \ TV 
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Chapter 4 



The Determination Of Vector Meson 
Twist-2 And Twist-3 Parameters 

In this chapter we determine the leading twist-2 and twist-3 two- and three-particle vector 
meson DA parameters using the non-local modification of SVZ sum rules. The parameters 
are defined in Chapter [5] and the sum rule method is outlined in Chapter [31 We express the 
relevant correlation functions, via the OPE, in terms of the perturbative and condensate 
contributions. Key to the analysis is the inclusion of all G-parity and SU(3)F-breaking 
effects which, as discussed in Chapter [21 come from a variety of sources, and allow a 
consistent determination of the parameters for the p, K*, and <fi. Motivation for the 
present analysis comes from various sources, including: 

• values for the decay constants and leading-twist DA Gegenbauer moments are 
required as input for QCD factorisation frameworks which provide a systematic 
method for the calculation of B decay matrix elements. We discuss one such frame- 
work in Chapter [HI 

• Twist-2 and twist-3 DAs provide the leading non-perturbative input within the 
method of LCSR, as discussed in Chapter El and as such are applied to many 
problems in heavy-flavour physics, such as the calculation of B transition form 
factors and the estimation of B decay matrix elements including power-suppressed 
contributions to QCD factorisation frameworks, see Chapter [71 

• A full determination of the twist-3 DA parameters, including SU(3)p-breaking and 
G-parity violating effects, and the inclusion of 0(a s ) and 0{m 2 s ) corrections to the 



18 



quark condensate contributions to the twist-2 DA parameter sum rules are new to 
the present analysis, allowing a| ,_L (0) to be determined, for the first time, to the 
same accuracy as a| ,_L (p, K*). 

All input parameters for the sum rules, and useful formulas, such as those required to 
take the imaginary parts of intermediate results, and various relevant integrals, are given 
in Appendix [B] In performing the calculations we find Refs. [99, 100] very useful. The 
material covered in this chapter partially follows that of Ref . [40] . 



4.1 Twist-2 

In this section we focus on the determination of the twist-2 DA Gegenbauer coefficients 
a\' ± defined by Eqs. (12. 25ft and ( 12.28ft . The sum rules for f^t, including SU(3)F-breaking 
corrections, were calculated in Refs. [61,62,114]. Those for the G-parity violating a^^K*) 
in Refs. [61,62] and those for a!j (K*) in [59] apart from perturbative terms in m 2 and 
the 0(a s ) and 0{m 2 s ) corrections to the quark condensate, which are new to the present 
analysis. Motivation for including these corrections is found by examining the individual 
contributions to the sum rules for ajj' 1 " (K*) given in Ref. [59]. They are found to be 
dominated by (ss) as we can see from the following explicit break down of contributions: 

PT (lf G2 ) (s 9s <TGs) (qq) 2 (ss) (ss) 2 (qq)(ss) 

al(K*) = ao5+ ao8 + an + o^-OA6 + om- cTo5 

= 0.06+ 0.10 + 0.25 +0.03 - 0.27 + 0.02 -0, (4.1) 

for the reference point s = 1.2 GeV 2 , M 2 = 1 GeV 2 and /j, — 1 GeV. Moreover, for the <p 
the impact of a finite strange quark mass may be even more pronounced with respect to 
perturbation theory and the gluon condensate. 

Firstly, we give an overview of the calculation of the 0(a s ) and 0(m 2 ) corrections to 
the quark condensate (ss); the calculations for (qq) are analogous. We only need extract 
terms proportional to m s as the contributions proportional to m q are identical; we can 
find the contributions for by simply replacing (qq) — > (ss) and doubling the terms 
in m s (ss), m s (qq) and m s (sg s Gs). Contributions for p are found by setting m s — > 0. 
Secondly, we go on to analyse the sum rules for a\' ± (</>). We end this section by presenting 
the results. 
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4.1.1 Calculation 



For both polarisations we begin from the diagonal correlation function 



U 2 -K*(q ■z)=i 



d 4 ye-^(0\Tq(y)Ts(y)s(0)nO,z]q(z)\0) 



(4.2) 



where r" = 7 2 and T 1 - = a^ z . For the longitudinal parameters the sum rule is exactly 
that given by Eq. (13.191) with fj — > f^* and for the transverse parameters the sum rule 
is analogous. Both polarisations have the same projections onto the DA parameters 
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where Im s denotes taking the imaginary part with respect to s. The fact that we are 
dealing with non-local correlation functions means that we do not integrate over the co- 
ordinate z. The resulting residual exponential function remains throughout the calculation 
and can contribute to the momentum integrals yielding powers of ic(q ■ z), where c is a 
constant. Ultimately the exponential functions can be cast into the "canonical form" set 
by the exponential appearing in front of the leading- twist DA i.e. e~ mq ' z - see Eq. (13.1 9ft . 

Quark Condensate 




Figure 4.1: The leading-order diagram contributing to the quark condensate {ss). 



The tree- level diagram is shown in Fig. 14.11 To extract the quark condensates to 0(m 2 s ) 
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we use the following expansion of the quark fields (for general quark flavour q) 
<0|:^(*iK(* 2 ):|0> = ^ 



- ™ g -^(7A)/3aA ; 



777 

2(2 + D) q 



(4.4) 



where A M = (x2 — £i) M and i,j are colour and a,/? spinor indices. One can deal with the 
co-ordinate A M by trading it, via partial integration (PI), for a derivative of the trace that 
arises from the quark loop. A convenient way to do so is via an auxiliary momentum Q 







(4.5) 



Diagrams for the 0(a s ) corrections to the strange quark condensate are shown in Fig. 14.21 









Figure 4.2: Diagrams contributing to the quark condensate (ss) at 0{a s ). The crossed circle <8> 
depicts the emission of a gluon from the non-local gauge factor - see Eq. (|4.6p . The corresponding 
diagrams for (qq) are identical but reflected top to bottom. 

Importantly there are contributions from the gauge-factor which need to be included 

[0, z\ = P exp | -ig s dt z^A^tz) j = 1 - ig s dt z^A^tz) + ... . (4.6) 

Calculating 0(a s ) corrections leads to divergent diagrams and the dependence of the 
condensate on the spacetime dimension D leads to 0(e) contributions at tree level, that 
then cause finite counter-terms upon renormalisation. Also, the derivative with respect 
to Q K in Eq. (14. 5p yields r y K in the trace via Eq. (13.341) which can also give a finite counter- 
term. This happens for the vertex correction diagrams. 
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4.1.2 Evaluation of The Sum Rules 



The new quark condensate contributions are added to the results presented in the litera- 
ture, see Refs. [40,61]. For f^t the sum rules read 
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and for a\' {K*) 
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To obtain the sum rules for /J' and a 2 (0), one has to substitute (gg) — > (ss) and 
to double the terms in m s (ss), m s (qq) and m s (sagGs) , and replace the perturbative 
contribution by 



for (f^) 2 : 
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(4.11; 



We have derived sum rules for the decay constants /]>' , however, numerical values can be 
extracted from experiment for the longitudinal decay constants. The perpendicular decay 
constants, on the other hand, must be determined from non-perturbative methods; results 
are available from Lattice QCD calculations and previous QCD sum rule determinations. 
A detailed discussion of the latest numerical values of the decay constants can be found 
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in Ref. [70] from which we just quote the following 

fl = (215 ± 5) MeV, # = (186 ± 9) MeV, (4.12) 

where /jj is an experimental result, and f£ is from Lattice QCD [115]. We can compare 
these results to the sum rules of Eqs. (14. 7p and (14. 8 p which are plotted in the upper row 
of Fig. 14.31 The sum rule determinations of a| ,_L (0) are plotted in the lower row. 




i ■ ■ ■ t ■ ."".'T 



1.2 1.4 1 .6 1 .8 2. 1.2 1 .4 1.6 1.8 2. 

M 2 M 2 

Figure 4.3: The decay constants f\ (upper left) and (upper right) and the Gegenbauer 
coefficients a|(0) (lower left) and a^i^) (lower right) plotted as a function of M 2 . The continuum 
thresholds are sjj = 1-85 ± 0.05 GeV 2 and = 1.40 ± 0.05 GeV 2 - see text. Solid line: central 
input parameters of Tab. IB. 1L Dashed lines: variation due to the uncertainties of m s and the 
gluon condensate. All quantities are evaluated at fi = 1 GeV. 

In all the plots the dashed line and shaded region represent the central value and uncer- 
tainty of the parameter in question. To evaluate the sum rules we use the input parameters 
of Tab. IB. 11 For the continuum threshold we note that for the sum rule determination 
of in Ref. [61] it is taken to be Sq(K*) = 1.7 GeV 2 , and we expect for it to be 
slightly larger. Indeed, by taking sjj(^) = 1-85 ± 0.05 GeV 2 we find a stable plateau and 
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excellent agreement with the experimental result for f) (upper left plot). Likewise, guided 
by sfr(K*) = 1.3 GeV 2 [61] we find s^(<p) = 1.40 ± 0.05 GeV 2 yields a result consistent 
with that from Lattice QCD (upper right plot). We use these thresholds in evaluating the 
sum rules for a| ,_L ((f)) and also replace the decay constants by their sum rules, which helps 
reduce dependence on the Borel parameters. The results are plotted for a\(4>) (lower left 
plot) and a^"(0) (lower right plot). It is found that the longitudinal parameters exhibit 
a stronger continuum threshold dependence, which is reflected in the larger uncertainty 
of the determined value of a| (</>). The sum rule determinations of the other particle DA 
parameters follow analogously and all the numerical results are given in Tab. 14.11 
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Li = 1 GeV 
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At = 2 GeV 










0.03(2) 


0.02(2) 
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0.04(3) 


0.03(3) 
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0.15(7) 


0.10(5) 


0.11(9) 


0.08(6) 


0.18(8) 


0.13(6) 




0.14(6) 


0.11(5) 


0.10(8) 


0.08(6) 


0.14(7) 


0.11(5) 



Table 4.1: Results for the twist- 2 hadronic DA parameters at the scale fi = 1 GeV and scaled 
up to = 2 GeV using the evolution equations (|2.57p . Note that a|'" L (if*) refers to a (sq) bound 
state; for a (qs) state it changes sign. 



4.2 Twist-3 

In this section we determine the twist-3 three-particle parameters of the DAs , &\-k* 
and $\.k* as defined by Eq. (12.411) . Previous determinations of these parameters are rather 
few and far between, thus motivating the present analysis. The chiral-even p parameters 
dp) u \ P i an d oj\ p were obtained in Ref. [96], and uo^ p was obtained in Ref. [55]. We make 
a comparison with these results in Section 14.2.21 

Firstly, we outline the calculation of the three functions it^k* which all proceed in a 
similar manner, and secondly we explicitly discuss the sum rules for §\. K » and present 
the results. In the diagrams that follow, q is the upper line and s is the lower line. 
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4.2.1 Calculation 



Each DA is accessed via a correlation function featuring its defining current. The chiral- 
even twist-3 parameters (l K *, w\ K *, \\ K * can be determined from 



^U-iw*) = (q ■ zf{2- D) J Ae-^(0|Tg(z)^G-M 7z 75 5 (0)%) 7 ^(2/)|0), 

(4.13) 

where the definition of g^ v is given in Appendix [All 1 ! The parameters k\ k * , w\ K * and A^* 
can be obtained from the correlation function nl,^, obtained from nj.^* by making the 
replacement 

g s G a zizi5 -> g s G az i-f z . (4.14) 
Lastly for the chiral-odd operator 

Ili. K .(v,q-z) = — ^ / d%e- iq - y (0\Tq(z)a z ,g s G z ,(vz)s(0)s(y)a qz q(y)\0}. (4.15) 



(<?• z ) 3 j 

All three correlation functions II can be written as 

n 3; x*M ■ z) = J Vae-^ z ^ +va ^7r 3 , K 4a) , (4.16) 

where the exponential function is due to the fact that we keep the correlation functions 
non-local. The calculation proceeds for each correlation function analogously. Considering 
Eq. (I4.13P for instance, firstly we express it in terms of hadronic contributions 

fiU(M^) = ( 4* )2m T / V(a)e- i ^ + ^^ K ,(a) + ... , (4.17) 

where the dots denote contributions from higher-mass states. To derive the sum rule we 
tread down a well worn path; express Eq. (14.161) as a dispersion relation and equate to 
Eq. (I4.17p . subtract the continuum contribution for s > s , perform the Borel transforma- 
tion and project out the desired DA parameter by substitution of the relevant polynomial. 



1 We also make use of the relation 7^75 = ^e M Aiy7r7 A 7 iy 7 ir denned in D dimensions. 
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The three hadronic parameters Qk*> uj \k*i *3K* are projected out like so: 
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The formulas for the other parameters are analogous. In calculating the functions tt^-k* 
we keep explicit mass corrections 0(rn^,m^,m s m q ) and all operators up to D = 6 except 
the triple gluon condensate {g 3 fG 3 ) which is expected to yield a negligible contribution. 
By retaining all mass terms the resulting formulas for tt^-k* can be used to derive sum 



rules for all the DA parameters for K*, p and <fi by setting m q = 0, m q 



m.. 



and 



m s respectively. For p and expressions for the three-particle twist-3 DAs are 



analogous to Eq. (I2.4ip . except that the G-parity violating parameters k and A vanish. 



Perturbation Theory 

The perturbation theory calculation is given by the two diagrams shown in Fig. 14.41 As an 




Figure 4.4: Diagrams contributing to perturbation theory. 

example, consider the first diagram, which up to an overall factor can be written generally 

as 

■[l„D vP (l)- l v D^(l)]e iz < w+p K (4.19) 

where the Dirac matrices depend on the correlation function. In performing the 
two successive integrations over / and p, Feynman parameterisation leads to shifting the 
variables I — > I — px and p — > p — qy respectively. Each time the exponential in (14.191) 
is also shifted. In expanding the part of the exponential that contributes to the integral, 
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for example, for I we have e ll ' zv = 1 + i(l • z)v + . . . , only the first two terms contribute; 
higher order terms are killed off either via z 2 = or because integrals with odd numbers 
of open indices, for example / Ml / M2 / M3 , in the numerator vanish due to symmetry. After 
the integrations any terms (T) including factors of i(q- z)v are dealt with by trading them 
for derivatives of T by using partial integration of the final exponential 

l ( q . z )v = - e -^y^-™) < q . z )vT^ -— T, (4.20) 

y ox yox 

where surface terms do not contribute as they vanish for x = {1,0}. The exponential can 
be matched to the "canonical form" by writing 

dx [ dye~ iq -^- SfS) = f dx [ dy f Va 6(^^)6 (a 2 -xy)6(a 3 -xy) e^*^-™^ . 
Jo Jo Jo J 

(4.21) 

Performing the x and y integration of the whole expression gives the desired result 

J Vae- iq - z(a2+va3) n {PTl \a) . (4.22) 

The second diagram follows analogously. Both diagrams are divergent and need to be 
renormalised separately. We find finite counter terms which are proportional to the quark 
masses. 



Gluon Condensate 

The leading order contribution to the gluon condensate (^G 2 ) is found using the back- 
ground field method as outlined in Section [331 There are only two diagrams contributing 
as depicted in Fig. 14.51 One vacuum momentum k, from the gluon attached to the quark 
line, is introduced and hence one derivative is taken. As the gluon emerging from the 
non-local vertex G(vz) carries no momentum these diagrams are proportional to 5(as) 
and the remaining momentum fractions are related by 1 — a\ = c^; the identification of 
the momentum fractions with the Feynman parameters is therefore straightforward. The 
calculation requires the integration over one momentum p and the result can simply be 
written unexpanded in the quark masses. 
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Figure 4.5: Diagrams contributing to the gluon condensate (^-G 2 ). 



Mixed Condensate 



The mixed condensates (qug s Gq) and (sag s Gs) originate from the diagrams shown in 
Fig. 14.61 To extract the mixed condensates one uses the first non-local term in the 





Figure 4.6: Diagrams contributing to the mixed condensates (qag s Gq) and {sag s Gs) 



expansion (D = 4) [99] 



(0|:gi(x 1 )^(G /1I ,)y(y)gj(x 2 ):|0> = 

' (qg s aGq) 



(4.23) 
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3q 



The first er^ does not contribute, but the term ~ m q does. The A M s can be expressed 
as derivatives of the trace via partial integration which is dealt with simply by using 
Eq. (13.341) . Along with the condensate gluon, the quark condensate lines carry no mo- 
mentum. There is therefore no loop integration to perform and the results are proportional 
to 5(a 3 )5(ai >2 ). 



Quark Condensates 

The diagrams of Fig. 14.71 generate the condensates m q ^ s (qq) and m q ^ s (ss). We do not 
consider 0(m 2 qs ) corrections, which are however of dimension six, as they are very well 
suppressed with respect to the other contributions. To extract all 0{m qs ) mass corrections 
the first non-local term in the expansion of the quark fields, given by Eq. (I4.4p . is needed. 
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There is one loop momentum to integrate over and one finds contributions from the 
exponential which can be dealt with via partial integration in the same way as with the 
perturbation theory calculation, see Eq. (14.20p . The results are proportional to S(ai^)- 
The diagrams in Fig. 14.81 generate the condensate (qq) (ss) which is already of dimension 
six, so we do not require mass corrections. The two diagrams are of equal magnitude 




Figure 4.8: Diagrams contributing to the quark condensate {qq) (ss). 

and cancel, however only for n\. K , they add. There is no loop integral to perform and 
the result is proportional to 5 (0.1)6(0.2). The four quark condensate is simplified via the 
vacuum saturation hypothesis (VSH) [14, 99] 

(0|: g l(xi)^(x 2 )4(x 3 )4(^):|0) ^3 (OhgKxOgjCx^ilOXOh^xsMC^):^). (4.24) 

The diagrams in Fig. 14.91 generate the condensates (qq) 2 and (ss) 2 . They stem from 
the operator (6) appearing in the expansion of the mixed condensate, Eq. (14.241) . which 
simplifies as 

(6) = (qiXq 97^3) ^ ~ <<?<?> 2 ! (4.25) 

u,d,s 

thus at higher order the mixed condensate also contributes to the quark condensates. The 
light-like co-ordinate of the gluonic field strength tensor vz^ simplifies the resulting trace 
via z 2 = from Eq. (I4.24p and the other co-ordinates are dealt with as before. 
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Figure 4.9: Diagrams contributing to the quark condensates (qq) 2 and (ss) 2 from the expansion 
of the mixed condensate - see Eq. (|4.24|) . 



Results 



For the functions 7r 3; ^*, given by Eq. (14.191) . we find (dropping all terms that vanish upon 
taking the imaginary part): 
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4.2.2 Evaluation of The Sum Rules 



In the following we consider tv%. k * ; the sum rules for the other DA parameters and particles 
p and follow similarly. The values of the input parameters and the continuum thresholds 
used for all sum rules are given in Appendix [B] 

One subtlety must be noted: upon integration over and subsequent expansion in powers 
of the quark masses, the gluon condensate contribution yields terms in m 2 s ln(m 2 s / {—q 2 )), 
which are long-distance effects and must not appear in the short- distance OPE of the 
correlation functions of Eqs (14. 13)) and (14.151) . The appearance of these logarithmic terms 
is due to the fact that the expressions of Eqs. (I4.26tl4.28l) are obtained using Wick's 
theorem which implies that the condensates are normal-ordered: (O) = (0|:O: |0) [116]. 
Rewriting the OPE in terms of non-normal-ordered operators, all infrared sensitive terms 
can be absorbed into the corresponding condensates. Indeed, using, 

o 

(0\sg s Gs\0) = (0| : sg s Gs : |0> + ^ In ^ (0| : — G 2 :|0) , (4.29) 

2 jji 71 

and the corresponding formula for q quarks, all terms in In m 2 s can be absorbed into the 
mixed quark-quark-gluon condensate and the resulting short-distance coefficients can be 
expanded in powers of m 2 . 

In Fig. 14.101 we plot the sum rules for \l K ,, w\ K * and CIk*i gi ven by Eqs. (I4.19p . which 
are evaluated for the central input parameters of Tab. IB.ll and at a scale fi = 1 GeV. The 
parameters unfortunately exhibit very strong M 2 dependence, which leads to increased 
uncertainty of their values; we do not find a stable plateau in the region 1 GeV 2 ^ M 2 ^ 
2.5 GeV 2 . On the other hand, there is only a very small sq dependence ~ 1% over the 
range 4(K*) = (1.3 ± 0.3) GeV 2 . The curves flatten at high M 2 which is expected, 
as the power corrections become negligible compared to the perturbative contribution^ 
The sum rules for the other parameters and particles show the same general behaviour 
which is fairly typical of non-diagonal correlation functions. If one were to use diagonal 
correlation functions then it is possible that the sum rules would be better behaved and 
thus the uncertainties would be reduced somewhat. The calculation of diagonal correlation 
functions of three-particle operators, as we saw with the gluon condensate in Chapter [31 
is rather more involved, especially when calculating radiative corrections, which may very 
well be necessary in this case. 

2 The quark condensates survive as M 2 — * oo as B [q~ 2 ] = — 1 but perturbation theory ~ M i - see 
Appendix [Bj 
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All the numerical results, including the uncertainties from the variation of M 2 , so, and 
input parameters, are given in Tab. 14.21 The results are presented at the scale \i — 1 GeV 
and scaled up to /i = 2 GeV, using the evolution equations, Eq. (12.571) . The only previous 
determination for comparison is for the chiral-even p parameters, Cj P (l GeV) = 0.033 ± 
0.003, Jj p (lGeV) = 0.2, and £ 3 ' p (lGeV) = -0.1 [96] and 0^(1 GeV) = 0.3 ± 0.3 [55]. 
These results agree with ours, although we consider the uncertainty of Q to be optimistic. 
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-0.025(20) 


-0.015(10) 









Table 4.2: Results for the leading three-particle twist-3 hadronic parameters of the DAs of 
Eq. (|2.4ip . The results are presented at the scale ll = 1 GeV and scaled up to fi = 2 GeV using 
the evolution equations (|2.57p . The sign of the parameters corresponds to the sign convention 
for the strong coupling defined by the covariant derivative = — ig a A ( ^ J t a ; they change sign 
if g s is fixed by D^ = d^+ ig s Ap a . 

In Fig.mjjwe plot the two-particle twist-3 DAs as defined by Eqs. ( !2~471 - [23U|) . G-parity 
violating effects cause the small asymmetry of the K* curves. The effects of SU(3)f- 
breaking are larger and cause the pronounced difference between and 03 for the p 
and 0. We notice in particular the end-point behaviour of the DAs is greatly modified. 
The fact that both and (p%!~ K * diverge as u — > 1 and (jr^ for u — > is in itself not a 
problem. It is only the leading-twist DA that can be considered a probability distribution 
and likewise there is no cause for concern that <p\. p takes negative values. Moreover, 
in practical calculations we are only interested in convolutions of the DAs with hard 
scattering kernels, which are generally finite. If not, this signals a problem with the hard 
scattering kernel, rather than the DA, as happens with end-point divergences within the 
QCD factorisation framework for non-leptonic B decays, see Chapter [HI 
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Figure 4.10: Hadronic parameters of the twist-3 distribution amplitude ^3.^-. as functions 

of M 2 . Upper: \\ K *, middle: £> 3 #*, and lower: Q K *- The solid curve is for central input 
values for fi = 1 GeV and outer curves take into consideration their uncertainties - see Tab, IB. ll 
Horizontal dashed line is the extracted DA parameter value and shaded region its uncertainty - 
see Tab. H2J 
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Figure 4.11: Left: (f)\ as a function of u for the central values of hadronic parameters, for 
fj, = 1 GeV. Red line: (j)\. p , green: 4>\. K ,, blue: 4^.^- Right: same for ^J. 
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Figure 4.12: Left: <^g as a function of u for the central values of hadronic parameters, for 
fj, = 1 GeV. Red line: (j>s ;p , green: blue: ^(b- Right: same for ^3 . 
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Chapter 5 

B — > rjy) Form Factors in QCD 



In this chapter we discuss the semileptonic B — > 77W form factors f+^ ( ) in the LCSR 
approach. The previous LCSR determination of the B — > 77 ^ form factors presented 
in Ref. [110] is completed by calculating the gluonic contribution, the mechanism for 
which involves the annihilation of the B meson to two gluons. The 77 W particles undergo 
pronounced mixing with each other due to the U(1)a anomaly of QCD and the 77-77' system, 
after many years of investigation, has succumbed to the phenomenologically motivated 
mixing scheme proposed by Feldmann, Kroll and Stech [117,118]. The consideration of 
this mixing scheme is central to the correct description of the B — > 77 ^ form factors. 

Motivation to complete the calculation of /+ ^ T?( ' comes from a variety of sources, with 
probably the most prominent being: 

• the flavour- singlet contributions to the QCD factorisation framework to be discussed 
in Chapter [6] were added by Beneke and Neubert in Ref. [120]. It is found that the 
branching ratios of B — > r]'(V, P) are very sensitive to f+^ v{ ' as the leading-order 
annihilation diagrams can be interpreted as a gluon contribution to the B — > 77 W 
form factors [121]. Therefore a consistent estimation of the annihilation diagrams 
necessitates the inclusion of the gluonic contributions to the form factor. 

• There exists a "tension" in the determinations of \V U b\ from inclusive semileptonic 
decays B — > X u lv and their exclusive counterparts, namely from B — > irlv. The 
former have led to larger values than the latter, and the reason for the discrepancy 
is unclear. B — > 77 W transitions are at leading order a b — > u transition and so 
sensitive to | which can, in principle, be extracted from B — >• if-'^lv. An improved 
calculation of f+^ r ' ( ' would reduce the theoretical uncertainty of the result. 
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• Finally, the observation that exclusive B —* rj'K and inclusive B — > rj'X decays 
have shown unexpectedly large branching ratios with respect to B — > 7r transitions, 
for example, is an unresolved issue which an improved calculation of ' may 

help clarify. 

We begin by introducing the 77 ^ system and define two closely related 77-77' mixing schemes. 
We then discuss the calculation of the flavour- singlet contribution to the form factor before 
lastly we discuss the results of the LCSR analysis, the framework for which was covered 
in Chapter [31 The material presented in this chapter follows that of Ref . [65] . 



5.1 The 77-77' System 

The approximate chiral symmetry of light quarks u, d and s in QCD seems to be broken 
by Nature to reveal the pseudoscalar mesons (ir°, 7r + , 7r~, K + , K~ , K°, K°, 77) as the cor- 
responding octet of Goldstone bosons (all massless in the chiral limit m u ^ s — > 0) of the 
broken SU(3) <S> SU(3) symmetry. There is another symmetry of the QCD Lagrangian 
(II. ID : a global U(1)a symmetry which exists at the classical level in the chiral limit. Due 
to non-vanishing quark masses, the broken U(1)a symmetry creates a Goldstone boson, 
but such a light particle does not appear in the physical spectrum and this embodies 
the U(1)a problem. At the quantum level, however, the U(1)a symmetry in the massless 
limit is broken due to the QCD anomaly and so was not present in the first place; thus a 
ninth state, the 77', exists as a singlet and only becomes massless in the chiral limit and as 
N c — > 00, causing the effects of anomaly to vanish. The situation is complicated by instan- 
ton effects, but was ultimately resolved by 't Hooft with the same conclusion [122, 123]. 
It has been known for a while that the U(1)a anomaly plays a decisive role in the rj^ 
system with the 77' consisting of a large gluonic component [124, 125]. The large mass of 
the 77' is mostly generated by the anomaly and SU(3)F-breaking effects^] 

The 77-77' system has been of considerable interest for a number of years [126-128]. Vast 
simplifications can be made in studying the low-energy particle spectrum of QCD by 
employing Chiral Perturbation Theory (ChPT) which is an effective theory in which the 
heavy quarks are integrated out and the dynamically relevant light quarks remain at a 
scale fi ~ Aqcd after an expansion in powers of energies, momenta and quark masses. 

lr rhe particles rj^ have masses m v = 547.51 ± 0.18 MeV and m v > = 957.78 ± 0.14 MeV and quantum 
numbers J PC = 0~+ [27]. 
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Alongside the 1/N C expansion, ChPT is the method of choice for analysing the light 
pseudoscalar mesons^] We do not discuss ChPT in any detail although we do quote a few 
of its constraints; for more details see for example [130-132] 

Concerning 77-77' mixing, ChPT requires a description in terms of two mixing angles beyond 
leading-order [133,134]. How this is implemented in practice has caused some confusion 
in the past but a consistent picture has emerged [117,118]. Key to the phenomenological 
picture of the 77-77' system is the understanding that the main contributions to the mixing 
are due to the U(1)a anomaly of QCD, and so-called OZI-rule violating processes. Named 
after Okubo, Zweig and Iizuka the OZI-rule states that strong interaction processes that 
must proceed via the annihilation of all initial state quarks to gluons are suppressed 
[135-137]. In Fig. 15.11 we show the unsuppressed process <fi —* K + K~ (left) alongside 
the suppressed process — > 7r + 7r~7r° (right) for which the rule was originally formulated. 
Such processes are shown to be 0(1/N C ) in a 1/N C expansion and phenomeno logically they 
are found to be small ~ 10%; they can be safely neglected, leaving the U(1)a anomaly 
as the only mixing mechanism. For the mixing schemes we discuss in the next section, 
this assumption has been confronted with experimental data and holds to the expected 
accuracy. 

(PS! 

Figure 5.1: Examples of strong interaction decays. Left: <fi — > K + K~, right: (ft — > 7r + 7r~7r . 
The former occurs preferentially over the latter due to the fact that the annihilation of the <fi 
requires all gluons to be hard, yielding a suppression via a small a s which need not be the case 
for the first decay. This forms the basis of the OZI-rule. 

A schematic picture of the U(1)a anomaly at work for B — > 77W i s shown in Fig. 15.21 . 
where the flavour-singlet contribution is defined as the amplitude for producing either a 
quark- ant iquark pair in a singlet state which does not contain the 5's spectator quark, 
or two gluons, which then hadronise into an 77 W. 

What about mixing between other pseudoscalar mesons? In 77 - 77' - ir° mixing the gluonic 
component present in the ir° is found to be at the level of a few percent and so can be 

2 Another interesting approach to understanding the 77 M system was given in Ref. [129]. 
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Figure 5.2: B — ► 77W via the U(1)a anomaly. The b — > u transition allows for an annihilation 
of the .B meson's quarks to two gluons, thus probing the gluonic content of 77^. 

neglected [118,138,139]. There also exists a cc component to 77 ^ (?7 C ) which is considered 
in Ref. [117] and found to be small with the conclusion that it is not the solution to the 
abnormally large B — * Krf branching ratio. Sometimes other particles are included as 
possible glueball candidates produced via OZI-rule suppressed processes in J /ip decay, see 
for example Refs. [125,140]. Although it is unclear whether pseudoscalar mesons contain 
pure glueball properties, Ref. [141] concludes that it is unlikely. Thus the 77-77' system 
stands out on its own. 

Phenomenologically, the semileptonic decay B — > rf'Hui can be used to determine the size 
of the CKM matrix element \V u b\ from the spectrum 

j-p r<2 it/ 12 

where P = {77,77'} and Xp(x) = (m? B + m 2 p — x) 2 — 4m|mp. Alternatively, as we shall 
see, the ratio of branching ratios B{B — > rj'£u)/B(B — > f]£u) can be used to constrain the 
gluonic Gegenbauer moment B%. 

5.2 State Mixing 

The first step in describing 77-77' mixing is to decompose the two physical states |t7^) into 
other, more convenient orthogonal states. As proposed in Refs. [117,118] one can proceed 
in two ways; either by employing the singlet-octet scheme (SO) or the quark-flavour 
scheme (QF). The SO axial- vector currents are respectively 

J% = ~j= (^7/,75M + dj^d + sj^ 5 s) , Jj 5 = -j= (uj^u + dry^d - 2s-f^ 5 s) , 

(5.2) 
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and their couplings are given by 



(0\JU\ p (p)) = ifpP* (< = 0,8), (5.3) 

where J^ 5 denotes the SU(3)p-octet and J° 5 the SU(3)p-singlet axial-vector current. The 
four quantities are related to the decay constants of a pure singlet or octet state \r]i) by 
two mixing angles 9i 



fv $ \ _ / cos ^ -sin^o W /s 



# /J y V sin 8 cos^o 7 V o /, 



(5.4) 



o 



Evidently SU(3)p-breaking effects cause 9t ^ and / 8 7^ /„■, and as such the SO scheme 
is very natural. In fact, at leading-order in ChPT an expansion in quark masses and 1/N C 
gives [133] 

MOo - 9s) = 2 ^ { / K ~P +■■■, (5-5) 

where the dots denote neglected higher-order terms which are required to match phe- 
nomenology [142]. The impact of the U(1)a anomaly is plainly localised in /o via the 
divergence of the singlet current J° 5 which can be written 



P J; 5 = 2 q [t a mi l5 ] q + S°° ^GG , (5.6) 

where a = {0, 1, . . . , 8}, Tr[t a t b ] = ^5 ab , t° = I/a/3 and the mass matrix rh = diag[m u , m^, m s 
The SO scheme diagonalises the renormalisation-scale dependence of parameters; /§ an d 
9i are scale- independent, whereas fo renormalises multiplicatively 



"|-^(v) 2a+0 ^- < 57 > 



In the QF mixing scheme, on the other hand, the basic axial-vector currents are 



and the corresponding couplings are 

(0\j; 5 \P(p))=if P p, (r = q,s). (5.9) 
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The mixing is analogous to (15. 4p with 



Both quark flavour states \f] q ^) have vanishing vacuum-particle matrix elements with the 
opposite currents 

(0|J^> = <0|^ 5 k) = 0, (5.11) 

which is an assumption that has been tested. It is in part motivated by the observation 
of near ideal mixing in vector and tensor mesons. It implies that the mixing of states is 
the same as that of the decay constants and moreover leads to the diagonalisation of the 
mass matrix, which we come back to shortly. This hypothesis does not hold for the SO 
basis. It is found by Refs. [134,138] that the difference between the two mixing angles of 
the QF scheme <p q — 4> s is generated by OZI-rule suppressed processes and is not caused 
by SU(3)F-breaking effects, as for the SO scheme (15.51) . While the numerical values of 
9i differ largely, with typical values 8% ~ —20° and 6 m —5°, one finds <p s — <f) q < 5°, 
with <p q « <p s « 40° [117, 118, 134]. This observation led the authors of Refs. [117, 118] to 
suggest the QF scheme as an approximation to describe 77-77' mixing, based on neglecting 
the difference (ft q — <p s (and all other OZI-breaking effects): 

4> = 9>s , <j) q -(j) s = 0. (5.12) 

The state mixing is then given by 

\v) \ _ ( cos ^ -sin</> \ / \r} q ) 
IV) J \ sin0 cos0 J y \rj s ) 

The renormalisation-scale dependence of fo given by Eq. ( 15. 7ft is not reproduced as it 
is induced precisely by neglected OZI-breaking terms [138]. Numerically, this is not a 
problem as the scale-dependence of fo is a two-loop effect. In the case of non-local 
matrix elements, the DAs, this lack of scale dependence of the QF scheme is somewhat 
problematic. We come back to this point in the next section. 

Returning to the diagonalisation of the mass matrix; from Eq. (15. 3p one finds the quadratic 
diagonal mass matrix, for example 

(0\d»J^\ V (p)) = M 2 v f°, (5.14) 



■ sin ( 
cos< 




fs 



(5.10) 
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which, via Eq. (15.61) . gives the mass matrix in QF basis 



* f(0\%GG\ V .) ml + f s (0\^GG\ Vs ) 



with the short-hand notation 

v^2 2 

m \ q = -^{0\m u uij 5 u + m d dij 5 d\r} q ) , m 2 ss = — (0\m s si'j 5 s\r] s ) . (5.16) 

Jq Is 



From Eq. (15.151) the crucial impact of the anomaly, as the only term in the off-diagonal 
elements, is evident. To first order in SU(3)p-breaking, the decay constants and quantities 
m qqss are fixed giving the theoretical estimate 



fq — fn i fs — y2f K f% , 

m 2 qq = M 2 n , m 2 s = 2M 2 K -M 2 , (5.17) 

which also leads to a fixed value of 0; there is no free parameter left and thus the QF 
scheme is totally determined [117]. We do not work in this limit, however, and take 
numerical values of the decay constants and mixing angle from phenomenology. Given 
enough data to fix all independent parameters, there is no reason to prefer the QF over 
the SO scheme. The QF scheme is beneficial when considering DAs as the SO scheme 
leads to a proliferation of unknown parameters. For this reason we decide to use the QF 
scheme for the analysis. Its basic parameters have been determined as [117,118] 

/,= (1.07 ±0.02)/,, f s = (1.34 ±0.06)/,, = 39.3° ± 1.0° . (5.18) 

This can be translated into values for the SO parameters as 



l/^ + |tf = (1.26 ±0.04)./; 



fo = \j\fq+\fs= (1-17 ±0.03)/., 

# 8 = <f> - arctan[V2/ s // ? ] = (-21.2 ± 1.6)° , 

9 = - arctanfv 7 ^//,] = (-9.2 ± 1.7)°, (5.19) 

Note that in the QF scheme / 9)S are scale-independent parameters, and so is fo as obtained 
from the above relations. The SO decay constants are related to those of the QF scheme 
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by a change of basis 




(5.20) 



The last matrix originates from the ideal mixing angle #idcai = arctan \/2 which rotates 
from the QF basis to the SO basis. 



5.3 Pseudoscalar Meson Distribution Amplitudes 

As discussed in Chapter [3j the method of LCSRs relies on the non-perturbative universal 
light-cone DAs; specifically here we require pseudoscalar meson DAs including the two- 
gluon DA. At leading-twist both these DAs contribute and indeed mix with each other un- 
der renormalisation. The quark-antiquark DAs are extensions of the matrix elements given 
by Eqs. (15.31) and (15.91) to those of non-local operators on the light-cone. Pseudoscalar 
mesons' quark-antiquark DAs have been investigated previously in Refs. [51,57,77]. The 
two-gluon DAs of leading and higher twist have been investigated in Ref. [143]. In this 
analysis we only include the effects of the leading-twist two-gluon DA, which is justified 
as its effects turn out to be fairly small and higher-twist DAs are estimated to have even 
smaller impact. Following Ref. [144], the twist-2 two-quark DAs of rf'> are defined as 

(0\^(z)C lzl5 [z,-z]^(-z)\P(p)) =i(p-z)f i P [ due^% P (u). (5.21) 

Jo 

4>2-p( u ) 1S the twist-2 DA of the meson P with respect to the current whose flavour 
content is given by C l , with \I> = (u, d, s) the triplet of light-quark fields in flavour space. 
For the SO currents, one has C° = l/y/3 and C 8 = \/2t 8 , while for the QF currents 
C q = (V2C° + C 8 )/V3 and C s = (C° - V2C 8 )/V3- Due to the positive G-parity of 77 and 
t]', the two-quark DAs are symmetric under u «-> 1 — u, and hence all odd Gegenbauer 
moments vanish: 

<t>iA u ) = 4p(! - u ) > ( 5 - 22 ) 
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and the DAs are expanded in terms of Gegenbauer polynomials in exactly the same way 
as for the vector mesons 

cj>% P (u) = 6u(l -u)ll+ a n l ^) C n 2 (0 J (i = 1, 8, q, s) , (5.23) 

\ n=2,4,... / 

where the quark Gegenbauer moments. The gluonic twist- 2 DA is defined aq_l 

(0\G, z (z)[z,-z}G^(-z)\P(p)) ^\{p-zf^f P j\ue^% P {u). (5.24) 

In order to perform the calculation of the correlation function defined in the next section, 
we also need the matrix element of the meson P over two gluon fields. Dropping the 
gauge factor [z, —z] one has 



1 z p jf C F n S ab r .* W%. P (u) 
The two-gluon asymptotic DA is n 2j ~ 1 (l — n) 2 -?" 1 with j = 3/2 the lowest conformal spin 

5 /2 

of the operator G^ z and the expansion goes in terms of Gegenbauer polynomials C n , see 
Eq. (jl.2ip . One can show that ip^-p * s antisymmetric: 

V4pH = - «) , (5.26) 

and in particular duip2-p( u ) = an d the local twist-2 matrix element {Q\G ^ Z G^ Z \P) 
vanishes. The non-vanishing coupling (0\G a gG a ^\P) induced by the U(1)a anomaly is a 
twist-4 effect. The corresponding matrix elements are discussed in Refs. [117,118] and 
are given, in the QF scheme, by: 

(Ol—^GGI^g) = f s (m 2 v - m*,) sin0cos0, 

(0\-^GG\t] s ) = U ml - ml) /V2sin0 cos <f> . (5.27) 

In taking the ratios of both sides of the above relations one can see that SU(3)p-breaking 
in the decay constants f q / f s is driven by the anomaly. There are no twist-3 two-gluon DAs 
and the remaining twist-4 DAs also have vanishing normalisation [143]. The conformal 



3 This definition refers to the "a- rescaled" DA <j) a g in Ref. [144] with a = \fl/C F - It agrees with that 



used in Refs. [143,149], which means that we can use their results for the two-gluon Gegenbauer moment 

2 



_B| without rescaling 
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expansion of the twist-2 two-gluon DA reads 



2;P 



u 2 (i 



u 



£ ^(A*)CK(0 

ra=2,4,... 



(5.28) 



with the gluonic Gegenbauer moments B^' 9 - In this analysis, we truncate both (f>\. p and 



ljj 9 . p a ^ n 

B 



2. An estimate of the effect of higher Gegenbauer moments in </> 2:7r on the 
7r form factor has been given in Ref. [86], based on a certain class of models 
for the full DA beyond conformal expansion. The effect of neglecting a^ >4 was found to 
be very small ~ 2% hence we expect the truncation error from neglecing B^ >A to be of 
similar size. 

<p2-p and ip2-p m i x upon a change of scale \i and as discussed in Refs. [144, 145] this 
amounts to a mixing of a^'° and B^' 9 , resulting in the renormalisation-group equation to 
LO accuracy 

d 



dp, 



( <$ N 






1 = 4vr 



100 
IT 

-36 



10 \ 

81 
22 



B 9 



(5.29) 



where for simplicity we have dropped the superscript P. The solution for a° reads 



_ ] ^7 2 + /(2/3o) 



+ 



1 49 

2 2^2761^ 



9 V / 276T 




^/(S/So) 



(5.30) 



with the anomalous dimensions 7^ = (149 ± V2761)/9. The octet Gegenbauer moment 
does not have another DA with which it can mix and so its evolution is simpler 



(5.31; 



The mixing amongst the DAs complicates matters; as the scale dependence of the decay 
constants is lost in the QF scheme, one expects to have to lose scale dependence in the 
DAs too, and we must be careful to be consistent. The verification of the anomalous 
dimensions in Eq. (15.291) from the singlet and octet parts of the form factor calculations is 
a crucial test of the LCSR analysis. For this reason, we discuss the implications of mixing 
on the twist-2 DA parameters, and only briefly cover higher-twist quark DAs which are 
included in the octet part; for a detailed discussion one is referred to Ref. [65]. Following 
Ref. [144], for the DAs introduced by Eq. (15.231) we have, in terms of the quark valence 



76 



Fock states \qq) and \ss) 

|%) ~ <j>l(u)\qq) + 2 OZI («)|^) , |77 S ) ~ |?5> + , (5.32) 

where qq is shorthand for (u« + dd) / y/2 and 

# = i($ + 2$), # = i(2$ + $), 0^ = ^(0°-^). (5.33) 

In the QF scheme, the "wrong-flavour" DA 0° ZI , which is generated by OZI-violating 
interactions, is set to 0. Once this is done at a certain scale, however, the different 
evolution of a° n and a\ will generate a non-zero <^ ZI already to LO accuracy. A consistent 
implementation of the QF scheme hence requires one to either set a^ 8 = and also 
B g n = 0, or to set a\ = a° n and neglect the different scale-dependence of these parameters. 
The induced non-zero DA <^ ZI is numerically very small for the scales relevant for our 
calculation, /x = 1 GeV and 2.4GeV@ The left panel of Fig. 15.31 shows a plot of A = 
100 K^O-O — a 2(A i ))/ a 2(/ i )l as a function of scale /i, according to Eqs. (15.301) and (15 .3 1 f) . 
for a|(l GeV) = a°(l GeV) and B\ = 0. We see that A is less than 0.25 % over the range 
lGeV < n < 2.4 GeV. Choosing a|(l GeV) = 0.25 ± 0.15, guided by our knowledge of 
twist-2 DAs of the tt; we have af (2.4 GeV) = 0.171 from Eq. (OIL and a°(2.4GeV) = 
0.171 for B 9 2 = 0, from Eq. (15.301) . Evidently, the impact of the different anomalous 
dimensions of a 2 an d d\ is negligible. Also, the evolution of a 2 is not hugely different to 




Figure 5.3: Left: A = I00\(a 2 (fi) — I ' a 2^)\ as a function of scale /i, according to 
Eqs. (POD and (JEST]) with B 9 2 = 0. Right: dependence of a§(2.4GeV) on Sf(lGeV) for 
a° 2 {\ GeV) = 0.25 according to Eq. (OOP 

4 2.4 GeV is a typical scale in the calculation of form factors from LCSRs: [i = \/rn? B — ml is chosen 
as an intermediate scale between nib and the typical hadronic scale 1 GeV. 
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that of a 2 , for a wide range of values of B 2 The right panel of Fig. l5.3l shows the evolution 
of the singlet Gegenbauer moment a 2 from — \ GeV - 2.4 GeV, from Eq. (15.301) . for the 
range of gluon Gegenbauer moments |5f(lGeV)| < 20, which is a very conservative 
estimated range, as discussed below. The mixing of B 2 into a 2 is up to 20% for B 2 = 20 
and 40% for B% = -20. 

From the conclusions of the above discussion we are justified in implementing the QF 
scheme for DAs as follows: we set <p 2 = 4> 2 at the scale \x = 1 GeV, which, by virtue 
of Eq. (15.331) . implies <fi 2 = <p 2 at the same scale. We then evolve 02 according to the 
scaling-law for the octet Gegenbauer moment (15.311) PI We also set ip 2 . v = again any 
SU(3)F-breaking of this relation is expected to have only very small impact on /+^ T? . 
The twist-2 parameters used in our calculation are then reduced to two: 02 and B 2 ■ 

Concerning numerical values, we assume that the bulk of SU(3)p-breaking effects is de- 
scribed by the decay constants via f q ^ f n , and that SU(3)F-breaking in Gegenbauer 
moments is sub-leading [57]. Sum rules for a 2 and a\ would essentially be the same, 
with /,r 7^ f q driving the SU(3)F-breaking and any small differences in so and M 2 being 
negligible. This motivates setting a 2 = a 2 , with a 2 (1 GeV) = 0.25 ± 0.15 as an average 
over a large number of calculations and fits to experimental data [57]. 

For B 2 , however, no direct calculation is available. Results from fits to data have been 
obtained from the 7/7 transition form factor, yielding B 2 (l GeV) = 9 ± 12 [144], and 
the combined analysis of this form factor and the inclusive decay T(IS') — > rj'X yielding 
-Bf(1.4GeV) = 4.6 ± 2.5 [143]. Caution must be taken when considering these results as 
they are highly correlated with the simultaneous determination of a 2 and a 2 from the same 
data, yielding a°(l GeV) = -0.08 ± 0.04, a|(l GeV) = -0.04 ± 0.04 and <z°(1.4GeV) = 
02(1.4 GeV) = —0.054 ± 0.029, respectively. The same analysis, applied to the 717 form 
factor, returns a^lGeV) = —0.06 ± 0.03 [147]. These results are not really compatible 
with those from the direct calculation of a 2 from Lattice QCD and QCD sum rules; in 
particular the sign of a 2 is unambiguously fixed as being positive. A possible reason 
for this discrepancy is the neglection of higher-order terms in the light-cone expansion 
and that, in addition, as one of the photons in the process is nearly real with virtuality 
q 2 ~ 0, one also has to take into account long-distance photon interactions, of order 
lyfq 2 , as discussed in Ref. [146]. For this reason, we assume the very conservative range 
(2.4 GeV) = ± 20 in the analysis. 

This is equivalent to imposing the QF-schcme relation the scale [i = 2.4 GeV and defining 

B 9 2 as B 9 2 (2.4 GeV). 
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As far as higher-twist quark DAs are concerned, we only need those involving currents 
with flavour content qq = (uu + dd) / y/2. In line with the implementation of the QF 
scheme for twist- 2 DAs, we include SU(3)F-breaking only via the decay constants. The 
precise definitions of all twist-3 and 4 DAs, as well as up-to-date numerical values of the 
7r's hadronic parameters can be found in Ref. [57]. A discussion of the correct treatment 
of these DAs within LCSR, as modified to describe rf'\ can be found in Ref. [65]. 



5.4 Calculation 

We define the B — > P form factors analogously to those of other pseudoscalar mesons 
as [110] 

(P( M B(P + ,)> = {(2 P + ,)„ - =^ ,„} ^ + „ . 

(5.34) 

where the factor of l/y/2 on the right-hand side is to ensure that in the SU(3)p symmetry 
limit, without 77-77' mixing, /+ = /£. For semileptonic decays B — > rj^lvi the form factor 
fo appears proportional to q 2 ~ m 2 which is negligible for light leptons / = {e, u} for 
which only f+ is required. Using the LCSR method outlined in Chapter [3] we extract the 
semileptonic form factor ff from the following correlation function 

KM = iJ*xe^(P(p)\T[&yM*)jU0)\V) (5-35) 
= n^(g 2 ,p|)(2p + g) M + ... , 

where js = rribui'j^b is the interpolating current for the B meson and p 2 B = (p + q) 2 its 
virtuality. In calculating the correlation function, we use Eq. (15. 13ft which relates the 
physical states and the QF basis states |r/ 9iS ) so that 

= -L (n* cos - n* sin 0) , nj = -L ( n * sin + n* cos 0) . (5.36) 

The interpolating current uj^b only probes the uu quark component of the 77 W so II* 
vanishes to leading order in a s and at 0(a s ) is due only to gluonic Fock states of the 
meson. IK, on the other hand, receives contributions from both quark and gluon states. 
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The final LCSR for then reads 

"P( n 2\ rso 



e -^ M2 m y B i±^l = f so ds e -s/M* I Ims n p (Si g2) (5 37) 

with the usual sum rule specific parameters M 2 , the Borel parameter, and Sq, the contin- 
uum threshold. 



Quark Contribution 

The quark contributions follow from the studies already undertaken for the tt, for more 
details see Ref. [110]. We briefly cover the general features of the calculation to put the 
singlet contribution in context. The leading quark contributions to IT^ originate from the 
diagrams of Fig. 15.41 where first order 0(a s ) corrections are shown. The external quarks 
have momentum fractions up and (1 — u)p and are on-shell; p 2 = m 2 P . The two-particle 





Figure 5.4: The quark-antiquark contributions to the semileptonic B — > 77W form factors 
/+ {q 2 ) from light-cone sum rules. The top left diagram is the leading one, the others are 
0(a s ). The double line corresponds to the b quark and the dashed lines the injection of the 
weak vertex momentum q, and the momentum of the B meson ps- 



DAs are projected out by using the general spinor decomposition of quark fields 

M' b = \{i)ba{qq') - \{^)ba{qiim') + ^(i^baiqiW) - \{i^)ba{q^im') 

- o («V75)b {qio^lhq') ■ (5.3£ 
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The vacuum-meson matrix elements of each term above either vanish or yield a DA 
depending on the quantum numbers of the meson in question. For pseudoscalar mesons 
the leading- twist contribution comes from 7^75, whereas 275 and 207^75 give two-particle 
twist-3 contributions, and although the two-particle twist-3 contributions appear in the 
sum rules as formally l/m^, they are chirally enhanced by numerically large factors [51] 
and so are included in typical LCSR analyses [110]. Three-particle twist-3 and two- and 
three-particle twist-4 DAs are also included; all twist-2 and -3 contributions include 0(a s ) 
corrections twist-4 contributions are to tree level accuracy. The corresponding expressions 
yield 11+ , with the replacement f n —> f q . 

Gluonic Contribution 

In order to obtain the gluonic contribution to 11+ , one needs to calculate the diagrams 
shown in Fig. 15.51 The last diagram is divergent and the other two are finite. The gluon 




Figure 5.5: The leading diagrams for the flavour-singlet contribution to the semileptonic 
B — > 77W form factors from light-cone sum rules. The double line corresponds to the b quark. 
The dashed lines the injection of weak vertex momentum q, and momentum of the B meson 
interpolating current ps- 

fields are introduced in the standard way 

K Igiuon = 1 J w » V\ ^(P(p)|T[n 7 ,6](x)Mz 7 5H(0) SCf(w)Cf(y)\0) , 

with the usual interaction Lagrangian C^(x) = ig s [qi'-f a A'^t a qi](x) with q^ = {u, b} and 
the statistical factor S takes values 1 if q± ^ q2 and 1/2 if qi = q^. The integral is over 
each co-ordinate separately. To extract the gluon contribution we need the projection 
onto the twist-2 two-gluon DA, which can be read off Eq. (15.251) . which amounts to the 
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following replacement of the gluon fields (up to the numerical factor) 



A a a (w)A%y) ™ S ab e af3pa ^ f du ^ 2;P ^ e ip < uw+Uy) , (5.39) 



1 

p ■ z j uu 

where the separation z is light-like i.e. z 2 = (w — y) 2 = 0. Via partial integration we can 
simplify the resulting expression for n^| gluon ; the co-ordinate z is traded for a derivative 
of the hard scattering kernel with respect to the momentum of one of the emitted gluons; 
and the dot product l/{p ■ z) can be traded for an integral with respect to the DA 
momentum fraction. As the boundary terms vanish due to the leading-twist gluon DA 
being antisymmetric, the calculation takes a rather simple form: 



l 

= / du 

I gluon / 
'0 



dTPjup) 
d(up) p 



u n l,9 



vv 



(5.40) 



where T?(up) is the hard scattering kernel. Both the gluonic and quark contributions are 
renormalisation scale dependent. The relevant term concerning the quark Gegenbauer 
moment 02 is 

Ul ~ 18/ 9 a 2 (l + ^ y In ^) F(p% q 2 ) , (5.41) 

where F(p 2 B ,q 2 ) is a function of p\ and q 2 . The logarithmic terms in the convolution of 
the gluonic diagrams of Fig. 15.51 with ip^.p are 



fP ]^^ B9f ln l^^^ „2> 

9v^4vr 2JP m 2 



I? ~ ~ S B lfp ln £2 F{vl, q 2 ) ■ (5-42) 



By expressing f q via Eq. (15.201) in terms of /° and fZ, respectively, and inserting Eq. (15.411) 
into Eq. (I5.36p . one verifies that the renormalisation-group equation, Eq. (15.291) . is fulfilled. 
The twist-2 two-gluon contribution to the correlation functions 11+ , Eq. (15.361) . is given 
in terms of a spectral density as 

n JL-=r« fa ^# (5.43) 
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with the result being 



5 ml 



±m^s(164g 2 + 13s) - m 4 (82g 2 + 95s)} 
+ fif a s f P m b -A= K 7 92)(g ~ m ' } {5m 4 ± g 4 ± 3g 2 s + s 2 - 5m 2 (g 2 + s)} 



6^ (s - <? 

x ( 2 in^-pi-in4!>- 



5.5 Discussion 



For the evaluation of the LCSR, Eq. f)5.37p . as with any sum rule, optimum values of 
M 2 and so need to be found. The standard procedure [110] is to replace fs by its sum 
rule, derived via SVZ sum rules, thus reducing the dependence of the LCSR on m& for 
which we use the one-loop pole mass m b = 4.80 ± 0.05 GeV [88]. From the fs sum rule 
the optimum threshold parameter sq = 34.2 ± 0.7 GeV 2 is found, and this value is taken 
over to the LCSR. As mentioned before \x = 2.4 GeV is chosen as an intermediate scale 
between m b and 1 GeV. The Borel parameter is taken to be M 2 > 6 GeV 2 and is varied 
in the range 6 GeV 2 < M 2 < 14 GeV 2 to reflect the corresponding uncertainty. In Fig. 15.61 
we plot f+^ r '(0) and f+~^ n (0) respectively as functions of M 2 , making explicit the result 
of varying so by ±0.7 GeV 2 , a-z by ±0.15 and £>f by ±10. As expected, /+(0) is not very 




Figure 5.6: /+(0) (left) and /? (0) (right) as a function of the Borel parameter M 2 and 
various choices of input parameters. Solid curves: central values of input parameters and so = 
34.2 GeV 2 . Long-dashed curves: sq varied by ±0.7 GeV 2 . Short-dashed curves: 02(1 GeV) varied 
by ±0.15. Dash-dotted curves: -Bf varied by ±10. 
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sensitive to the gluonic twist-2 DA parameter B 9 2 (dashed-dotted curves), but is quite 
sensitive to the Gegenbauer moment a2 (short-dashed curves). For /+(0), on the other 
hand, the dependence on B 9 2 is more pronounced than that of 02- Varying all relevant 
parameters within their respective ranges, i.e. Amj = ±0.05 GeV, Aa2(lGeV) = ±0.15 
and A5| = ±20, as well as all twist-3 and twist-4 parameters within the ranges given in 
Ref. [57], we find 

M 2 s "2 B 2 M T3 T4 f B ,m b 

fl(0) = 0.229 ±frS)5 N ±0~bl^±0~o1l)±0^^ 

= 0.229 ±p_ : 024±p_ : 011, (5.45) 

param. syst. 

M 2 so 4 B a 2 f q ,<f> T3 T4 f B ,m b 

/£(o) = o.i88±o"^±o~S^±ol^±ol)^ 

= 0.188 ±0 : 043±0 : 019 ±(1009 . (5.46) 

B? 2 param. syst. 

The entry labelled "T4" also contains an estimate of the possible impact of the local 
twist-4 two-gluon matrix elements (15.271) . For this estimate, we exploit the fact that 
the asymptotic DA of the non-local generalisation of Eq. ( 15. 2 7ft is the same as for the 
twist-2 two-quark DA: 6w(l — w)E We then assume that the corresponding correlation 
function is the same as that for the leading conformal wave in the two-quark twist-2 
contribution, i.e. the coefficient in the Gegenbauer moment ao = 1, and replace ao by 
{0\a s GG / (47r)\ri q)S ) / (f q)S ml) . The factor 1/m^ comes from the fact that this is a twist-4 
effect and hence suppressed by two powers of mj with respect to the twist-2 contribution. 
This is only a rough estimate, of course, as the true spectral density will be different. The 
results ( I5.46P show that for small 5f ~ 2 both twist-2 and -4 two-gluon effects can indeed 
be of similar size. In this case, however, the total flavour- singlet contribution to f_{_ will 
also be small, ~ 0.008. In the third lines, we have added all uncertainties from the input 
parameters (param.) in quadrature and the sum- rule specific uncertainties from M 2 and 
s (syst.) linearly. For /+ (0), we have displayed the dependence on B% separately. The 
new result for /+(0) is, within errors, in agreement with the previous one from LCSR, 
/+(0) = 0.275±0.036, obtained in Ref. [110]. That for (0) is new to the present analysis. 
The results agree well with those obtained in Ref. [149], from pQCD, /+(0) = 0.208 and 
/+ (0) = 0.171, including a rescaling by a factor \pl to bring their definition of the form 
factors into agreement with Eq. (15.341) . We confirm the finding of Ref. [149] that the 



3 This follows from Eq. For G±x, one has I = 2 and s = 0. 
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range of the singlet contribution to the form factor estimated in Ref. [120] is likely to be 
too large, unless B\ assumes extreme values ~ 40. 
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Figure 5.7: (0)//+(0) as a function of the Borel parameter M 2 and various choices of input 
parameters. Solid line: central values of input parameters, which corresponds to f\ (0)//+(0) = 
tan<^> = 0.814. Dot-dashed curves: -Bf varied by ±10. Dotted curves: al' 11 (lGeV) varied 



independently: al = 0.1, 



0.4 and al 



0.4, a\ = 0.1. 



In Fig. 15.71 we plot the ratio f\ (0)//+(0) as a function of the Borel parameter. In the 
ratio, many uncertainties cancel, in particular that on As we have chosen B% = 
as central value, /+(0)//+(0) = tan</> = 0.814 exactly, see Eq. (15.361) . The figure also 
illustrates the change of the result upon inclusion of a non-zero B\ (dot-dashed curves). 
The ratio is actually rather sensitive to that parameter. While the dependence on 02 
largely cancels when a\ and al are set equal, there is a considerable residual dependence 
on al — al 7^ (dotted curves). While \al — al | = 0.3 as illustrated by these curves 
is rather unlikely, and would signal very large OZI-breaking contributions (recall that 
al 7^ al or, equivalently, a\ 7^ a\ signals the presence of "wrong-flavour" contributions 
to the r] q)S DAs and is set to in the QF mixing scheme), one should nonetheless keep 
in mind that moderate corrections of this type are not excluded and compete with the 
OZI-allowed corrections in B\. 

Finally, in Fig. 15. 81 we show the dependence of the ratio of branching ratios R m i = B(B — > 
rfev)/B(B — > rjeu) on JB|. The advantage of this observable is that all hadronic effects 
are encoded in the form factors and that \V u b\ cancels. The solid curve corresponds 
to the branching ratios obtained from the central values of input parameters; the long- 
dashed curves illustrate the dependence on parameters originating from the model used 
to extrapolate the q 2 dependence of the form factor from beyond the limit of the LCSR 
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Figure 5.8: The ratio of branching ratios R m ' = B{B — > rj ' ev) /B{B — > rjev) as a function 
of the singlet-parameter £?f. Solid curve: central values of input parameters; Long-dashed 
curves: including variation of extrapolation model parameters, See Ref. [65]. Short-dashed 
curves: theoretical uncertainty of R vv i for £>f = 0, for a™ (1 GeV) varied independently, as in 
Fig. 1^71 



approach, in this case q 2 = 16 GeV 2 , to the maximum possible value q 2 ^^ = (m B — m^,)) 2 . 
It may be noted that the dependence on these parameters is very small. We do not go 
into detail about the extrapolation procedure and refer the reader to Ref. [65]. On the 
other hand, R m i also depends on al ^ a\ ■ This dependence is shown by the short-dashed 
curves. The conclusion is that large values of Sf, |J3f| > 5, can be distinguished from 
the OZI-breaking parameter \a^ — al |, once an accurate experimental value of R vr) i is 
available, but that for smallish B 9 2 and unknown \a\ — a\ \ only mutual constraints on 
these parameters can be extracted from the data. In this case also twist-4 gluonic DAs 
can become important. 
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Chapter 6 

QCD Factorisation 



In this chapter we discuss the framework of QCD factorisation which was introduced in 
the context of exclusive two-body non-leptonic B decays by Beneke, Buchalla, Neubert 
and Sachrajda in Refs. [10,11]. We shall refer to the the original implementation of the 
framework as the BBNS approach. We also focus on its application to the radiative B 
decays B — > V 1 -/, as presented by Bosch and Buchalla in Refs. [16, 17]. 

QCD factorisation allows a rigourous determination of the B decay matrix elements of 
the weak effective Hamiltonian (11.221) to leading order in the heavy-quark limit of QCD 
mi, ^> Aqcd, and yields a neat factorisation formula. It relies on the factorisation of 
hadronic matrix elements into universal non-perturbative hadronic parameters, given by 
transition form factors and meson light-cone DAs, and process dependent hard-scattering 
kernels, calculable in perturbation theory. The validity of the QCD factorisation formula, 
to all orders in a s , and the impact of generally unknown power corrections, formally 
suppressed by powers of l/m b , must be addressed case by case. The introduction of 
the QCD factorisation framework has made more discerning phenomenological studies of 
exclusive B decays possible whereby key observables, such as branching ratios, CP and 
isospin asymmetries, can be calculated and confronted with experimental data. 

The dependence of the factorisation formula on meson DAs, either directly or via LCSR 
calculations of the transition form factors, greatly motivates their study, with their better 
determination reducing the theoretical uncertainty of the QCD factorisation predictions, 
and aiding the quest to discover new physics effects from decay observables. 

We begin with a short introduction, in the context of B — »• M\M% decays, of the general 
features of QCD factorisation, and in particular, discuss the appearance of meson DAs. 
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We then discuss the framework as applied to the radiative B decays B — > V7. We 
postpone all discussions of phenomenology to Chapter [7| in which we perform an analysis 
of the decays B u ^ — > (p, u, K*)^ and B s — > (K*, 0)7 using QCD factorisation, augmented 
by the inclusion of the dominant power-suppressed corrections. 

6.1 Introduction 

QCD factorisation (QCDF) [10,11] was introduced in the context of the "heavy-to-light" 
decays B — > nn where the factorisation of the relevant QCD matrix elements was shown 
to apply, to leading order in a l/m& expansion, to a large class of non-leptonic B decays. 
Consequently, QCDF has opened up the rich and varied landscape of B decays to a more 
complete quantitative analysis. The existence of factorisation in non-leptonic decays is 
non-trivial and complicated by the possible gluonic interactions amongst the initial and 
final states. Conversly, leptonic and semi-leptonic decays factorise much more easily into 
the product of a quark current and a leptonic current, which cannot interact via gluon 
exchange. 

Phenomenologically, QCDF has been remarkably successful, especially given the range of 
processes for which the method holds. After its introduction, it was swiftly generalised 
to encompass irK final states [12], pseudoscalar- vector final states [121] and vector- vector 
meson final states [150]. The gluonic flavour-singlet contributions to B — > K^'rf'' decays 
were added by Ref. [120]. To date, the framework has been extended to many other 
processes, including for example, (double) radiative B decays B — > 7(7, V) [17, 151] and 
B — > jlv [71]. Also, other factorisation frameworks have since been developed and applied 
to the same problems: 

• Soft Collinear Effective Theory (SCET) [152-155] makes a careful distinction be- 
tween a hierarchy of "hard" (mj), "hard-collinear" (a/Aqcd?^) and "collinear" 
(Aqcd) scales via contributions of internal quark and gluon lines. Details of the 
differences between the SCET and BBNS approaches to QCD factorisation can be 
found in Refs. [156-158]. 

• The Perturbative QCD (pQCD) approach [159], which yields a factorisation formula 
that depends on the mesons' transverse momenta. 

• The method of LCSRs, although having existed before the advent of QCDF, was 
applied to B — > titt, both to the matrix elements which exhibit factorisation and 
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also a class of power corrections, providing some useful complementary insights, see 
Refs. [160,161]. 

We now go on to discuss the general features of QCDF. 

6.2 General Structure 

Consider the case of non-leptonic decays where the B meson decays into two mesons. The 
simplest way of dealing with the resulting matrix elements is to employ naive factorisation 
[162,163]. Simply put, naive factorisation splits each local operator Qi of the effective 
Hamiltonian into two colour-singlet currents, whose matrix elements are proportional to 
a decay constant and a transition form factor respectively. For example, consider the 
four-quark operator Q% = (DU)v-a(Uo)v-a then 

(M l M 2 \(DU) Y _ A (Ub) Y _ A \B) ^ (M 2 \(DU)v. A \0) (MtHUtyv-AlB) . (6.1) 

V v ' N v ' 

The motivation for factorising in this way comes from the colour transparency argument 
[164]. It follows that a major shortcoming of naive factorisation is that it assumes the 
exchange of gluons of virtualites \i < to be negligible and hence rescattering between 
the decay products is not considered; there is then no mechanism for the generation of 
strong phase effects between different amplitudes. Also, the matrix elements (16.1 p do not 
display the correct renormalisation-scale dependence. 

The framework of QCDF allows the calculation of 0(a s ) corrections to naive factorisation, 
which occur at scales fi < m^. It is constructed by observing the cancelation of infrared 
(IR) and collinear divergences, via consistent power-counting arguments, allowing the 
use of perturbation theory to describe the hard-gluon exchanges. The resulting intuitive 
factorisation formula thus presents a massive simplification of the long-distance QCD 
effects, with QCDF recovering naive factorisation in the limit — * oo. In terms of two- 
body non-leptonic B decays to light pseudoscalar mesons B — > MiM 2 the factorisation 



s9 



formula, as presented in Ref. [10], reads schematically as 

(M 1 M 2 \Qi\B) = 

+ 
+ 

where F B ~ iM ^ is the relevant form factor, T { ' are the hard-scattering kernels, 4>b is 
one of the leading-twist DAs of the B meson and <p2;P the leading-twist DA of the final 
state meson P, and the Qi are the operators of the effective Hamiltonian. The matrix 
elements are given as the convolution of the universal DAs and the process dependent hard- 
scattering kernels, with respect to the meson momentum fractions. Since the transition 
form factor and the DAs are real functions, all strong phases are generated by the hard- 
scattering kernels and are suppressed by powers of a s . Factorisation has be proven to 
one-loop for "light-light" final states and two-loop for "heavy-light" final states [11]. It 
has be proven to all orders in a s for B — > Dn using SCET [153]. 

The ability of QCDF to accurately describe B decay processes is limited by two main 
considerations; firstly, by the nature of the factorisation formula itself, which is valid up 
to power corrections 0(l/rrib) and to a given order in a s ; and secondly by uncertainties 
of the necessary input parameters, such as the DAs, the transition form factors, the 
strange quark mass, the B meson decay constant fg etc. Whether a discrepancy between 
experiment and QCDF predictions can be put down to new physics, or not, requires an 
estimation of neglected power corrections; certainly the b quark mass is not asymptotically 
large ~ 5 GeV and power corrections are therefore expected to feature at the level 
of 0(AQCT>/ m b) ~ 10%. The size and nature of power corrections can be probed via 
phenomenology, however, the task is not straight forward; even the initial focus of the 
approach, the decays B — > n(K, n), which stands as a crucial test, has not been resolved 
satisfactorily, see for example Ref. [165] and Refs. [166,167]. Better determined input 
parameters will nevertheless shed light, case by case, on whether power corrections are 
important, and the QCDF predictions must be used to determine or constrain CKM 
matrix elements (UT angles), or detect signs of new physics, with that in mind. 



F B ^ Nh [ duT^u) <j> 2 ; M M + (Mi " M 2) 
Jo 



d£ du dvT t H (£, u, v) 0b(£) fo;Ah( v ) 02 ; m 2 O) 
0(A QCD /m b ) (6.2) 
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6.3 Light-Cone Distribution Amplitudes 



To leading-order in the heavy-quark limit the leading-twist final state meson DAs con- 
tribute to the factorisation formula and can be safely truncated after the second Gegen- 
bauer moment a^- For pseudoscalar meson final states the two-particle twist-3 DAs come 
with large normalisation factors r£ and are said to be chirally enhanced, and are therefore 
included even though they are formally 1 / m;, suppressed. The vector mesons do not have 
the same large normalisation factors but their two-particle twist-3 DAs are included in 
the BBNS approach for consistency. For a pseudoscalar or vector meson, with valence 
quark content qq', the normalisation factors are respectively 



p, \ 2m p A QCD v ( v 2m y fy-fa) 

x m b (n)(m g + m q ,)(ij) m b x rribW /j> 



Three-particle twist-3 DAs are neglected because they do not come with large normalisa- 
tions. The inclusion of the chirally enhanced DAs leads to end-point divergences from the 
convolutions of the two-particle twist-3 pseudoscalar DAs with the corresponding hard- 
scattering kernels originating from both the hard-spectator scattering and annihilation 
contributions. The resulting divergent integrals signal the breakdown of factorisation and 
are parameterised by two universal unknown parameters Xh,a-, introducing a source of 
theoretical uncertainty to the BBNS approach [10]. 

At leading-twist the B meson is described by two DAs, only one of which is required as 
input for Eq. ( 16.21) and appears in the hard-spectator diagrams contributing to Tf 1 . The 
DAs of the B mesons are complicated by the fact that the momentum of the meson is 
shared in a highly antisymmetric way: the b quark has most of it. The B meson DAs are 
given, at leading-order in l/m&, by 

(0\q a (p)bp(z)\B(p B )) = i^f W B + m b )^] M j\ie-^+ z - [^(O+MiafOU > 

(6.4) 

with the decay constant fg given by Eq. f)3.2ip . With a careful choice of n_ = (1, 0, 0, —1) 
only the following normalisation conditions are required 



/ d£<M£) = l, f ^$ B2 (O=0 
Jo Jo 



(6.5) 
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along with the first inverse moment of $£i which is parameterised as 

f\*ifL s ^L, (6.6) 

and the numerical value of As is a source of uncertainty in the QCDF framework for both 
B — > MiM 2 and B — > Vj. We now discuss the radiative decays B — > V'-f within QCDF. 



6.4 Radiative B decays to Vector Mesons 

We consider the leading contributions to the B — > V7 QCDF factorisation formula as of 
Refs. [16,17,148,168] in which a model independent framework is presented. Contributions 
that are power-suppressed by one power of l/m& or more and are 0(a s ) are not considered. 
At the quark level the decays are b — > Dj transitions, where D = {s, d}. If otherwise 
not stated, in the following we refer to B — > V7 decays where B (V) denotes a bq (Dq) 
bound state. For B — > V7 decays the matrix element of each relevant local operator in 
the effective Hamiltonian factorises as 

(V-y\Qi\B) = e* ■ 

where e M is the photon polarisation vector and T^ v (0) is the relevant form factor. 
the leading-twist DA of the perpendicularly polarised final state vector meson (12.281) : 
contributions from <p\. v are power-suppressed in the heavy-quark limit. Problems of end- 
point divergences are not encountered in B — > decays and the twist-3 vector meson 
DA does not feature - the B meson DAs (16.61) do however. The factorisation formula is 
accurate up to corrections suppressed by powers of l/m,b, as shown, and was proven to 
hold to all orders in a s in SCET [169]. The form factor T^ v (0) has been calculated, for 
example, from LCSR in Ref. [112]. 

The B —* V'j decay produces either left- or right-handed photons, which therefore consti- 
tute, in principle, two separate observable processes. In practise the direct measurement 
of the photon's helicity is very difficult; indirectly, however, it can be accessed by mea- 
surement of the time- dependent CP asymmetry in B° — > V 0, ~f, which vanishes if one of 
them is absent, see Chapter [71 We define the two amplitudes as 

A m = A(B -> V lL{R) ) , A L{R) = A{B -> V lL[R) ) . (6.8) 
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If- v (0)l/+ f dtduTffo 
Jo 



u)(t>B{i)(t>iv{u) 



+ 0(l/m b ) , (6.7) 



For (B) B decays the production of the (left-) right-handed photon is suppressed by l/m b 
with respect to the opposite helicity. The decays are dominated by the electromagnetic 
dipole operator Q 7l , and as such are penguin mediated and so loop-suppressed. The 
operators Q 7 ^ are given by 

Qrj R) = A m b Da, v (1 ± 75) bF^ , (6.9) 

and generate left- (right-) handed photons. Their matrix elements can be parameterised 
in terms of the form factor as 

(V(p, V ) lL{R) ( q ,e)\Q^\B} 

= ~ mtTf^iO) [e^e; V tp p q a ± i{{e* ■ rj*)(p ■ q) - (e* ■ p)(rj* ■ q)}] 

= -^m b T^ v {Q)S L{R) , (6.10) 

where Sl.r are the helicity amplitudes corresponding to left- and right-handed photons, 
respectively, and e M (?7 M ) is the polarisation four- vector of the photon (vector meson). 
The leading-order diagram is given in Fig. 16.11 which is also the leading diagram for 
the form factor Tf^ v . The factorisation formula (I6.10p is therefore trivial to leading 



Figure 6.1: The leading contribution to B — > Vj due to the electromagnetic dipole operator 

order in a s and the heavy-quark limit; the matrix element given by the standard form 
factor, the scattering kernel by a purely kinematical function and Tj 1 does not feature. 
The electroweak penguin operators Q7 v ..,io appear at higher-order and safely neglected 
in the analysis. All other operators begin to contribute at 0(a s ). The hard- vertex 
corrections contribute to T/ yielding functions of Tn^ c /ml and originate from penguin 
contractions of the operators Qi,...,6 an d the chromomagnetic operator Q 8g as shown in 
Fig. 16.21 The hard-spectator scattering diagrams of Fig. 16.31 in which the spectator quark 
of the B meson participates, contribute to T/ 7 and involve the same operators as the 
hard-vertex corrections. The hard-gluon exchange probes the momentum distribution of 
the B and vector mesons and so requires the introduction of the mesons' light-cone DAs, 
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Figure 6.2: Penguin contractions of Qi,...^ (top line) and the chromomagnetic dipole operator 
Qs (bottom line) contributing to the hard- vertex corrections of T* at 0(a s ). Crosses denote 
possible photon emission vertices. 

as suggested by the factorisation formula; it is in these contributions that the B meson 
DA parameter \b and decay constants /# and fy appear. Also, the dominant power- 




Figure 6.3: Penguin contractions of Qi 6 (left) and the chromomagnetic dipole operator Qg g 
(right) contributing to the hard-scattering kernel Tf 1 at 0(a s ). Crosses denote possible photon 
emission vertices at leading order. Photon emission from the other quark lines power-suppressed. 
Photon emission from the final state meson for Q$ g breaks factorisation. 

suppressed weak annihilation (WA) contributions, shown in Fig. 16.41 are calculable in the 
QCDF approach, and involve the operators (5i,...,6- WA contributions are 0(l/mb); photon 
emission from the b quark and the quarks in the vector meson is further suppressed and 
0(l/ml) - unless the weak interaction operator is Qn } e, which can be Fierz transformed 
into (.0(1 + 75)9) (<?(1 — 75) and picks up an additional factor m B from the projection 
onto the B meson DA thus resulting in this contribution being 0(l/ra&). Consequently, 
due to the large Wilson coefficients these contributions are sizeable and important 
phenomenologically, see Chapter [7J 
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Figure 6.4: Weak annihilation contributions, which are suppressed by one power of 1/m;,. 
Crosses denote possible photon emission vertices at leading order. The dominant mechanism for 
Ql,...,4 is the emission of the photon from the light quark in the B meson and for Qs^q it is the 
emission from the final state vector meson quarks. Other possible emissions are either vanishing 
or more strongly suppressed. 



The decay amplitude is then given by 

A(B -> V lL[R) ) = ^ {\y u 7L(R) (V) + A? a c 7L(R) (V)) (V lL{R) \Q% R) \B) , (6.11] 



where the left-handed coefficients are given, to leading order in QCDF, by 

a £QCDF (y) =Cr + Q ^ 1/mb) ^ (g 12) 

and the right-handed parameters, for a b — ► D transition, by [170] 

a u,QCBF = Cj ^R + (l/m b , a s /m b ) . (6.13) 
rrib 

Explicit expressions for the 0(a s ) corrections to the left-handed coefficients can be found 
in Refs. [16, 17] and will be considered in Chapter alongside the dominant power- 
suppressed corrections. 
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Chapter 7 

B — > I/7 Beyond QCD Factorisation 



In this chapter we perform a phenomenological analysis of the exclusive radiative B de- 
cays to vector mesons. We make use of the QCDF framework outlined in Chapter [6] 
and investigate the impact of the leading power-corrections on the branching ratios, CP 
asymmetries and isospin asymmetries for all b — * D transitions; B U ( i — > (p,u, K*)'-/ and 
B s — > ((f), K*)^. Weak annihilation effects, although power-suppressed, are calculable in 
QCDF, and are included for all decay modes in this analysis. The other power-suppressed 
contributions "beyond QCDF" considered are; soft photon emission from the soft B spec- 
tator quark [72]; and long-distance contributions from heavy quark loops [170] and light 
quark loops [70] which have been estimated from LCSR. The estimation of the light quark 
loop contribution is new to the present analysis. Whereas the branching ratios are gen- 
erally dominated by the leading contributions, and power-suppressed contributions play 
a minor role, the same cannot be said for the CP and isospin asymmetries for which the 
impact of power-corrections is in fact crucial. 

The motivation to study radiative B decays stems from a variety of sources: 

• as loop-induced, penguin mediated decays, they allow the extraction of the CKM 
matrix element |V^(d s )| complimentarily to the determination from B mixing and 
also that from the SM UT analysis based on the tree- level observables IK^/V^I and 
the angle 7. 

• They are sensitive to new physics contributions, which may occur within the pen- 
guin loops, with the time-dependent CP asymmetry a very promising avenue of 
investigation. They are also subject to large short-distance QCD corrections, which 
now approach next-to-next-to-leading-order accuracy, see Refs. [186,187]. 
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• The decay rates are of order G^ojqed and are enhanced with respect to other loop- 
induced non-radiative rare decays which are of order G 2 F aQ ED . Also, the b — > s 
modes are CKM-favoured. Consequently there exist good experimental results for 
the exclusive branching ratios; B — > K*^ is known to 5%, but the b —>■ d transitions 
are not so well known. 

As discussed in Chapter El the QCDF framework for B — ► V'j relies on the leading- 
twist vector meson DA 0^y. Moreover, the LCSR calculations of the form factors 
and the parameters entering expressions for the soft-quark contributions rely also on the 
higher-twist DAs of the vector mesons and thus we find immediate use for the results 
of the twist-2 and twist-3 DA parameters of Chapter HJ as presented in Tab. 14.11 and 
Tab. Si 

We begin with an introduction, and then go on to discuss the power-suppressed contri- 
butions and investigate their impact on the decay observables. We extract the CKM 
parameter [V^/V^l from the branching ratio results, assuming no new physics contribu- 
tions, and discuss possible new physics contributions to the CP and isospin asymmetries. 
The material covered in this chapter follows that of Ref. [70]. 



7.1 Introduction 

B — > V'j decays are a very rich and promising probe of flavour physics. Both the inclusive 
decay B — > X s 7 and the exclusive decays B — > (K*,p)f have been under scrutiny for many 
years, see for example Refs. [171,172]. The experimental results for B — > (p, u,K*)j are 
shown in Tab. 17.11 For B s — > 07 only an upper bound B(B S — > 07) < 120 x 10~ 6 exists 
and no experimental information is available for B s — > K*-f [27]. 

In the SM the decays are flavour-changing-neutral-current (FCNC) b — > D^j transitions, 
mediated by penguin diagrams; they are therefore loop-suppressed and potentially very 
sensitive to new physics. To determine the relative sizes of contributions to the decays 
one must consider the following points: 

• the leading term is loop-suppressed ~ l/(47r) 2 and proportional to Cj ~ —0.3. 

1 The analysis presented in Ref. [70] used preliminary input for the DA parameters, values for which 
were later finalised in Ref. [40]. The conclusions and numerics of the analysis are unaffected, due somewhat 
to the large errors attributed to the soft quark loop calculations in which the twist-3 DA parameters 
feature. 
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B x 10 6 


BaBar [189] 


Belle [190] 


B x 10 6 


HFAG [37] 


B -> (p,u;j7 

-> p+7 
£?° -> p° 7 
5° ^^7 


1 OC+0.25 _i_ pi pin 

i.zo_o 24 ^ u.uy 

1.1018:33 ±0.09 

0.79l8:| ± 0.06 
< 0.78 


i OO+0.34+0.10 
L - oz -0. 31-0.09 

n C-C-+0.42+0.09 
U - JO -0. 36-0.08 

i 9 c-+0.37+0.07 
i - ZlJ -0.33-0.06 

n Q«+0.34+0.05 
u - JU -0.27-0.10 


D — > is. 7 


A n Q -1- O P. 
4U.O ± Z.O 

40.1 ± 2.0 



Table 7.1: Experimental branching ratios of exclusive b — > (d,s)j transitions. All entries are 
CP averaged. The first error is statistical, the second systematic. B — > (p, w)7 is the CP average 
of the isospin average over p and cj channels: 

B(B - (p,u,) 7 ) = \ {B(B± -+ p± 7 ) + ^ [8(1?° - p° 7 ) + £(5° - w 7 )] }. 

• Evidently from Eq. (16. lip for each mode there are two amplitudes proportional to 
different CKM factors A^ . For b — > d transitions both A« and Ac are ~ A 3 , 
however, for b — > s transitions A«^ ~ A 4 and A^ ~ A 2 ; there is a relative CKM 
suppression of the up-quark contribution. 

• Power suppressed corrections from WA are formally ~ l/m& although come with 
large Wilson coefficients C\ ~ —0.3 and C 2 ~ 1 and are not loop suppressed. The 
WA contributions drive the isospin asymmetries. 

• The production of "wrong" helicity photons is suppressed by m D /m b 06.131) . The 
interplay of both helicity amplitudes generates the time-dependent CP asymmetries, 
which are small in the SM due to this suppression. 

7.2 Wilson Coefficients 

Considerable effort has gone into calculating the Wilson coefficients to NLO accuracy. 
Using the expressions for the NLO anomalous dimension matrices available in the litera- 
ture we employ the renormalisation techniques of Eqs. f ll.24H1.3"li) to calculate the Wilson 
coefficients at the required scales. Numerical values of all the NLO Wilson coefficients C, 
used in the analysis are given in Tab. 17.21 The situation is complicated by the fact that 
the QCDF results of Ref. [17] are given in terms of two bases. The first, the so-called 
BBL basis named after the authors of Ref. [173], is that of Eqs. (11. 22ft and HI . 231) except 
with Qi and Q 2 exchanged with respect to the basis of Ref. [16]. The second is the 
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so-called CMM basis of Ref. [174,175]. The two bases differ except for QBBl = Q^ M . 
Following Ref. [17], the CMM set is used for calculating hard- vertex corrections to the 
QCDF formulas and the BBL set at the lower scale fih ~ V^-h H (with ~ 0.5 GeV 
and \x = 0(m b )) is used to calculate hard-spectator corrections. Power corrections are 
calculated from the BBL set at scale mj. 

NLO accuracy is mandatory only for C7, as it is for this term only that the hadronic matrix 
element is also known to NLO accuracy. We evaluate all 0(a s ) and power-suppressed 
corrections using both LO and NLO scaling for Wilson coefficients and hadronic matrix 
elements and include the resulting scale dependence in the theoretical uncertainty. 



C? MM (m b ) 


C 2 CMM (m b ) 


C 3 CMM (m 6 ) 


C 4 CMM (m b ) 


C 5 CMM (m 6 ) 




C 7 CMM (m b ) 


-0.322 


1.009 


-0.005 


-0.087 


0.0004 


-0.001 


-0.309 


Cf BL (m b ) 


C BBL (m b ) 


C BBL (m b ) 


C BBL (m b ) 


C BBL (m b ) 


C BBL (m b ) 


c 8 CMM K) 


-0.189 


1.081 


0.014 


-0.036 


0.009 


-0.042 


-0.170 




c BBL M 




c? BL M 








-0.288 


1.133 


0.021 


-0.051 


0.010 


-0.065 


-0.191 



Table 7.2: NLO Wilson coefficients to be used in the analysis, at the scales m b = 4.2 GeV and 
fj,h = 2.2 GeV. The coefficients labelled BBL correspond to the operator basis of Ref. [173] and 
given in Eq. (fL23|) . whereas CMM denotes the basis of Ref. [174]. We use a s (m z ) = 0.1176 [27] 
and mt(m t ) = 163.6 GeV [201]. Note that C BBL and C BBL are exchanged with respect to 
the basis of Ref. [16] and that Cf ( f } L = C° ( M M . Following Ref. [17], the CMM set is used for 
calculating hard-vertex corrections to the QCDF formulas and the BBL set at the lower scale 
fj-h is used to calculate hard-spectator corrections. The BBL set at scale m b is used for the 
calculation of power-corrections. 



7.3 Leading and Power Suppressed Contributions 

It proves convenient to split to the coefficients in Eq. (16. lip into three contributions which 
we will investigate separately: 

o?l(V) = a^ CBF (V) + aT n (V) + a^ i \V) + ..., 

*7r(V) = a^ CDF (V) + a^ nn (^) + * S ° ft (^) + .-., (7.1) 
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where the leading term in the l/m& expansion is given by Eq. (16.121) and all other terms 
are suppressed by at least one power of m&. The dots denote terms of higher order in a s 
and further l/m& corrections to QCDF, most of which are incalculable. We only include 
those power-suppressed terms that are either numerically large or relevant for isospin and 
CP asymmetries. 



7.3.1 Leading Contributions 

The diagrams giving the leading QCDF contributions are given in Chapter [HJ It turns 
out that, at the level of two decimal places, all equal and so are a^ QCDF @ For 

central values of the input parameters of Tab. 17.81 we obtain 

Vertex Hard-Spectator 
Corrections Corrections 



a c 7 § CDF {V) = - (0.41 + 0.03z) - (0.01 + O.Olz) , 

a^ CDF (V) = -(0.45 + 0.07z) + (0.02 - Oz) . (7.2) 



The size of the hard-spectator corrections is set by the factor 

27T 2 Ml 

nv 9 m B T^ v (0)X B - 



(7.3) 



For B s decays one has to set f B — » f Bs and correspondingly for the other B meson 
parameters. We estimate the value of X Bs , the first inverse moment of the twist-2 B- 
meson light-cone DA, from \ Bd by a simple scaling argument: 

^ ( Aqcd + ms ) = ^ A QCD , (7.4) 

A B S * Bq 

which follows from the assumption that the B q DA peaks at the spectator momentum 
k + = Aqcd, whereas that of B s peaks at Aqcd + m s - Its numerical value is given, along 
with all the other input parameters, in Tab. 17.81 



7.3.2 Weak Annihilation 

a^ a ™ encodes the 0{\jm^) contribution of the WA diagram of Fig. 17.1( a) which drives 
the isospin asymmetries and has been calculated in QCDF in Ref. [17] with a s corrections 

2 Explicit formulas for a^£^ CDF , complete to 0(a s ), can be found in Ref. [17]. 
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given in Ref. [176] for p and K* and in Ref. [148] for u. WA receives contributions from the 
current-current operator Q2, which for b — > s transitions is doubly CKM suppressed, and 
QCD penguin operators Q3,...,6, which are not CKM suppressed. Formulas for aj^ nn (p, K*) 
in QCDF can be found in Refs. [17,148]; in this approximation, there is no contribution 




(a) (b) 



Figure 7.1: (a) Weak annihilation diagram where photon emission from the B meson spectator 
quark is power-suppressed. The crosses denote possible photon emission vertices for Q§fi only, 
(b) soft-gluon emission from a quark loop, where there is also a second diagram in which the 
gluon is picked up by the B meson. 



Preliminary results for the 0(a s ) corrections to WA in B — > p7 were presented in 
Ref. [177]. In QCDF, the a^ nn are expressed in terms of the hadronic quantities 



Tf~ v (0) m B m„\ B ' " ^-"(0) m B m t J 

and d% , obtained by replacing 1/v — > 1/v in the integrand; (p^-v i s the twist-2 DA of a 
transversely polarised vector meson, (12.281) . Numerically, one finds, for instance for the 
p, b p = 0.22 and d p = -0.59, at the scale p = 4.2 GeV. As T x ~ l/m 3 b /2 and f B ~ m~ 1/2 
in the heavy-quark limit, these terms are 0(l/rrib), but not numerically small because of 
the tree-enhancement factors of tt 2 . 

For u, K* and we obtain 

a^ ann (^)| QCDF = Q d b w {a\ + 2(a 3 + a 5 ) + a 4 ) + Qd« + d%)a Q , 
a7'™ n (cj)| QCDF = Q d b UJ (2(a- i + a 5 ) + a 4 ) + Qd(d% + d%)a§ , 
a^ ann (0) = Q s b*(a 3 + a 5 ) + Q s (dt + dl)a 6 , 

QCDF 

a^ n (r) = Q s b R *a A + Q s (df* Q d /Q s + df)a 6 , (7.6) 

QCDF 
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WA 


B~ -> K*~ 


B° -> iT" 


B^(p,u) 


£? s ^0 


B S ^K* 


induced by 


C (and P) 


P 


C and P 


P 


P 


CKM 


A 2 (and 1) 


1 


1 


1 


1 



Table 7.3: Parametric size of WA contributions to B — > V7. C denotes the charged-current 
operators <2i,2> P the penguin operators Qz,...fi'-, their Wilson coefficients are small - see Tab. 17.21 
CKM denotes the order in the Wolfenstein parameter A with respect to the dominant amplitude 
induced by Q7. 

with ai = Ci + C 2 /3, 03 = ^3 + C4/3, a 4 = C 4 + C3/3, a 5 = C 5 + C 6 /3, a 6 = C 6 + C 5 /3@ 
The expressions for and ^* are new; for u, we do not agree with [148]. Apart from for 
p and u, all the WA coefficients are numerically small and do not change the branching 
ratio significantly; the terms in a^, however, are relevant for the isospin asymmetries. 

In Tab. 17.31 we show the relative weights of these diagrams in terms of CKM factors 
and Wilson coefficients. The numerically largest contribution occurs for B ± — > p ±r j: it 
comes with the large combination of Wilson coefficients Ci + C\/3 = 1.02 and is not 
CKM suppressed. For B° — > (p ,u)^ it comes with the factor C\ + C2/3 = 0.17 instead 
and an additional suppression factor 1/2 from the electric charge of the spectator quark 
[d instead of u). For all other decays, WA is suppressed by small (penguin) Wilson 
coefficients. Apart from B — > (p, cjj)j, WA is not relevant so much for the total values of 
ajLi but rather for isospin breaking, which is set by photon emission from the spectator 
quark. WA is the only mechanism to contribute to isospin asymmetries at tree-level; see 
Ref. [176] for 0(a s ) contributions. 

In view of the large size of a"£ mn (p) it is appropriate to have a look at further corrections. 
The most obvious ones are 0(a s ) corrections to the QCDF expressions, shown in Fig. 17.21 
As it turns out, the corrections to the B vertex in Fig. 17.2( a) are known: they also enter 
the decay B — > •yiu and were calculated in Ref. [71, 178]. Numerically, they are at 
the level of 10%. Fig. 17.2( b) shows the vertex corrections to the V vertex, which are 
actually included in the decay constant fy. For the non-factorisable corrections shown in 
Fig. 17.2( c) preliminary results have been reported in Ref. [177] according to which these 
corrections are of a size similar to the B vertex corrections. 

In Ref. [176] also another class of l/m& corrections to B — > K*^ was calculated, namely 

3 Note that a\ <-> as compared to [17] as in our operator basis (i.e. the BBL basis) Q\ and Q2 are 
exchanged. 
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(a) (b) (c) 

Figure 7.2: Example radiative corrections to weak annihilation. The corrections to the B 
vertex in (a) are known [71, 178] and those to the V vertex in (b) are included in fy. For the 
non-factorisable corrections in (c) only preliminary results are available [177]. 

0(a s ) corrections to the isospin asymmetry in this decay. As these corrections break 
factorisation (require an infra-red cut-off in the momentum distribution of the valence 
quarks in the K* meson) and are numerically small, we do not include them in our 
analysis. 

7.3.3 Long-Distance Photon Emission 

Another class of corrections is suppressed by one power of m& with respect to the QCDF 
contributions and is due to long-distance photon emission from the soft B spectator 
quark. A first calculation of this effect was attempted in Ref. [179] and was corrected and 
extended in Ref. [72]. The long-distance photon emission from a soft-quark line requires 
the inclusion of higher-twist terms in the expansion of the quark propagator in a photon 
background field, beyond the leading-twist (perturbative) contribution; a comprehensive 
discussion of this topic can be found in Ref. [69]. The quantity calculated in Ref. [72] is 

(p~ (ph(q) | (du) v _ A (ub) v _ A \B- (p + q)) = 

= e {F v e^el Pp q a - iF A [e*»(p ■ q) - q»{e* ■ p)}} 

TUB 

= -e^i^F Y (S L + S^ + ^F A (S L -S R )\ (7.7) 

in terms of the photon-helicity amplitudes Sl^B In QCDF, F A Y are given by QuJb/^b 
and induce a term Q u d2b p in a"£ inn (p~). The long-distance photon contribution to Fy A 

4 Eq. (|7.7p differs from the one given in [72] by an overall sign, which is due to the different convention 
used in [72] (and in [69]) for the covariant derivative: = — ieQfA^ instead of = + ieQfA^ 
as in this analysis. 
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was found to be [72] 



TJISoft 

b A 



with G A + G\ 



0.07 ±0.02 = Q U G A , 
G A 



= -0.09 ± 0.02 = Q U G V . (7.8) 



-0.24 ±0.06 and G v 

7L(R) 



>fl 



0.030 + 0.015H In order to obtain concise 



expressions for oHff^, it proves convenient to define one more hadronic quantity: 



9 P l,r 



7l 



mpfp 



T{ m b m B 



(G v ± G A ) 



(7.9) 



and correspondingly for other mesons, is 0(1/ ml) as Gy + G A has the same power 
scaling in m& as T\, i.e. ~ m b 3 ^ 2 , as one can read off from the explicit expressions in [179]. 
The difference Gy — Ga, on the other hand, is a twist-3 effect due to three-particle light- 
cone DAs of the photon and is suppressed by one more power of mj, i.e. gn ~ 1/wif- This 



quantity will enter the CP asymmetry. Our final expressions for a 7 j^ then read: 



f/,ann 
l 7L 

£/,ann 
l 7R 



(V) 
(V) 



(7,ann / 
l 7L 



(V) 



QCDF 



(7,ann / 
l 7L 



(V) (6 



QCDF 



(7.10) 



Numerically, one has g p L /b p = —0.3, so these corrections, despite being suppressed by 
one more power in l/m&, are not small numerically and larger than the known 0(a s ) 
corrections to QCDF from B ^lv. Based on this, we feel justified in including these 
long-distance corrections in our analysis, while dropping the radiative ones of Figs. 17.2( a) 
and (c). For central values of the input parameters we find the following numerical values 
for the various WA and long-distance photon contributions, including in particular those 
to which Qi^ contribute (with no Cabibbo suppression): 



a 7L 



l 7L 



(K*°) 

n V) 



b 7L 



-0.013- 0.001 LD, 
-0.001 - 0.004 LD, 
-0.024 + 0.003 LD. 



CWr (K 



b 7L 

it, ami 
l 7L 



(P 



0.004 + 0.001 LD, 
0.149 - 0.043 LD, 



(7.11; 



5 Again, there is a relative sign with respect to the results in [72]. This comes from the fact that the 
product eF A °y is independent of the sign convention for e, and as we have changed the overall sign of 
(|7.7p with respect to [72], we also have to change the sign of F A v . Stated differently: the relative sign 
between F A ° l v and F A y in [72] is wrong because of a mismatch in sign conventions for e in the covariant 
derivative. 
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The contribution from the long-distance photon emission is labelled "LD" (LD— > 1 at the 
end). The unexpectedly small dj^ nn (p ) is due to a numerical cancellation between the 
charged- current and penguin-operator contributions. Comparing these results with those 
from QCDF, Eq. (17.21) . it is evident that WA is, as expected, largely irrelevant for the 
branching ratios, except for B ± — > p ±r ~f. 



7.3.4 Soft Quark Loops 

a 7L(R) enc °des soft-gluon emission from a (light or heavy quark) loop as shown in Fig. l7.lf b). 
Soft-gluon emission from a charm loop was first considered in Ref. [182] as a poten- 
tially relevant long-distance contribution to the branching ratio of B — > K*~f, however, 
the same diagram also contributes dominantly to the time-dependent CP asymmetry in 
B° — > K* 0r y [181]. As for atf R , the dominant contributions to a^ R {K*) were calculated in 
Ref. [170] and new to this analysis is their generalisation to the other vector mesons and 
the inclusion of contributions from light-quark loops. Motivation to include light quark 
loops stems from the fact that they are doubly CKM-suppressed for b S7 transitions, 
but not for b —>■ dj, for which they are on an equal footing as the heavy quark loops. 
The quark loop contributions are suppressed by one power of m b with respect to a^ CDF , 
but they also induce a right-handed photon amplitude which is of the same order in 
l/m,b as a^jj? CDF (16.131) . and this amplitude induces the time- dependent CP asymmetry. 
The asymmetry is expected to be very small in the SM and cx ma/mi, due to the chiral 
suppression of the leading transition (16 . 1 3j) . but could be drastically enhanced by new 
physics contributions - thus constituting an excellent "null test" of the SM [7,170]. It 
was noticed in Refs. [180, 181] that the chiral suppression is relaxed by emission of a gluon 
from the quark loop and contributes dominantly to the time-dependent CP asymmetry 
in B° — > K* 0/ ~f, which motivates the inclusion of these contributions. The task of the 
present analysis, however, is not so much to calculate these contributions to high accu- 
racy, but to exclude the possibility of large contributions to the CP asymmetry. With 
this in mind, the theoretical uncertainties of the results are very generously estimated - 
which is somewhat unavoidable due to the current uncertainties of the relevant hadronic 
input parameters. 

Potentially the most important contribution to the soft-gluon emission diagram in Fig. 17. 1( b) 
comes from the charged- current operator Q\ with the large Wilson coefficient C2 ~ 1; 
it vanishes for QY by gauge invariance. In addition, the penguin operators Qs^q give 
a non-zero contribution. Details of the derivation of a^' soft can be found in Ref. [70] in 
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which the following expression is obtained: 



a 




(7.12) 



Here the sum over q runs over all five active quarks it, d, s, c, b. The contribution from Q 5 
is proportional to mo and hence helicity suppressed and neglected. Assuming SU(3)f- 
flavour symmetry for the light quark loops, one has l u — Id — h, and ditto for l Uy d, s , 
which causes a cancellation of these contributions in the last term of Eq. 07.12)) . SU(3)p- 
breaking effects are estimated to be around 10% [70]. The parameters l c (K*) and l c (K*) 
were first calculated from three-point sum rules in Ref. [182] and were re-calculated in 
the more suitable method of LCSR via a local OPE in Ref. [170]. The analysis therein 
as been updated and extended to lb, If, and the other particles p, u, K*, <fr for the present 
analysis [70]. The results for l c and l c are given in the upper table of Tab. 17.41 Those for 
lb and lb are obtained as 



For light-quark loops the photon is almost at threshold and the local OPE does not 
apply. In Ref. [70] a method was developed for calculating these contributions via LCSRs. 
Similar to the method of Ref. [160] used for the calculation of soft-gluon contributions to 
B — > 7T7T, a dispersion relation approach is used to connect the off-shell matrix element to 
the physical regime q 2 = 0. The results are presented in the lower table of Tab. 17.41 

7.4 Phenomenological Results 

In this section we combine the different contributions to the factorisation coefficients o,yl(r) 
and give results for the observables, namely the branching ratios, the isospin asymmetries 
and the time-dependent CP asymmetries. 




(7.13) 
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lc + lc 


B - 


» K* 


-355 ± 280 


-596 ± 520 


242 ± 370 


-952 ± 800 


B - 


» (p,w) 


-382 ± 300 


-502 ± 430 


120 ± 390 


-884 ± 660 


B s - 


-> K* 


-347 ± 260 


-342 ± 400 


-4 ± 300 


-689 ± 600 


B s - 


-0 


-312 ±240 


-618 ±500 


306 ± 320 


-930 ± 750 





lu 




lu lu 


f>u H~~ 


B - 


> K* 


536 ± 70% 


635 ± 70% 


-99 ± 300 


1172 ±70% 


B - 


> (p, to) 


827 ± 70% 


828 ± 70% 


-1 ±300 


1655 ± 70% 


B s - 


-> K* 


454 ± 70% 


572 ± 70% 


-118 ±300 


1025 ± 70% 


B s - 


-0 


156 ± 70% 


737 ± 70% 


-581 ±300 


893 ± 70% 



Table 7.4: Soft-gluon contributions from c-quark (upper table) and w-quark (lower table) loops 
in units KeV. The quantities l CiU and l C)U are denned in Ref. [70]. We assume equal parameters 
for p and to. lb is obtained as Z& = l c m 2 /ml and correspondingly for It,. The uncertainty for l u — l u 
is given in absolute numbers because of cancellations. In the SU(3)F-flavour limit assumed in 
this calculation one has l u = Id = l s = l„ 



7.4.1 Branching Ratios 



The (non-CP-averaged) branching ratio of the b — > Dj decay B — > is given by 



r B G 2 F a QET) m 3 B m 2 



B(B — > Vj) = 



■in 



2 \ 3 



X 



327T 4 

U =u,c 



111 



B 



U =u,c 



(7.14) 



with the isospin factors c p ±^* 4 



1 and c p o jU> 



a/2. The branching ratio for the CP- 



conjugated channel B — > V"f (b — > D7 at parton level) is obtained by replacing — > 
(A[f )*. With the input parameters from Tab. 17.81 and the lifetimes given in Tab. 17.71 we 
find the following CP-averaged branching ratios for B — > if* 7, making explicit various 
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sources of uncertainty: 



Ti n v ch iu,c other 
B{B~ ^ K*~-f) = (53.3 ± 132 ±^8 ±^8 ±^9^ ± 1.3) x 10~ 6 

= (53.3 ± 13£±5.8) x 10" 6 , 

Ti 

Ti m (u, c other 

B(B° -> AT*°7) = (54.2±T3^±^0 V ±l5 V ±^8 V ± l^T) x 10~ 6 

= (54.2 ± ^g±6.7) x 10" 6 . (7.15) 

We have added all individual uncertainties in quadrature, except for that induced by the 
form factor. The uncertainty in fi is that induced by the renormalisation-scale dependence, 
with /i = mb(mb) ± 1 GeV. "Other" sources of uncertainty include the dependence on the 
parameters in Tab. 17.61 on the size of LD WA contributions and the replacement of NLO 
by LO Wilson coefficients. The above results agree, within errors, with the experimental 



ones given in Tab. 17.11 within the large theoretical uncertainty induced by the form factor. 

As the uncertainties of all form factors in Tab. 17.81 are of roughly the same size, one might 
conclude that the predictions for all branching ratios will carry uncertainties similar to 
those in (17.151) . This is, however, not the case: the accuracy of the theoretical predictions 
can be improved by making use of the fact that the ratio of form factors is known much 
better than the individual form factors themselves. The reason is that the values given 
in Tab. 17.81 which were calculated using the same method, LCSRs, and with a common 
set of input parameters, include common systematic uncertainties (dependence on fg, 
etc.) which partially cancel in the ratio. In Ref. [183] the ratio of the K* and p form 
factors was found to be 

T B ~* K *(Q) 

Co = ^-b — = 1.17 ±0.09. (7.16) 

The uncertainty is by a factor 2 smaller than if we had calculated £ p from the entries in 



Tab. 17.81 analogously for to one finds 

_ 7f-^(0) 



1.30 ±0.10. (7.17) 



The difference between £ p and £ w is mainly due to the difference between and /±. For 
the B s form factors, we also need the ratio of decay constants fB s /fB d - The status of 
fs from Lattice QCD was reviewed in Ref. [184]; the present state-of-the-art calculations 
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are unquenched with Nf = 2 + 1 active flavours [185], whose average is /e s //s d = 1.23 ± 
0.07. Again, this ratio is fully consistent with that quoted in Tab. 17.81 but has a smaller 
uncertainty. One then finds the following ratios for B s form factors: 



T x (0) 
lf w (0) 



1.01 ±0.13 



n 



B^K* 



(0) 



(0) 



1.09 ±0.09. 



(7.18) 



The uncertainty of is smaller than that of ^ because the input parameters for K* 
and K* are the same (except for G-odd parameters like a±) and cancel in the ratio; the 
uncertainty is dominated by that of /s s //s d . To benefit from this reduced theoretical 
uncertainty in predicting branching ratios, one has to calculate ratios of branching ratios, 
which mainly depend on £v and only mildly on T\ itself: in addition to the overall nor- 
malisation, Ti also enters hard-spectator interactions and power-suppressed corrections, 
whose size is set by hadronic quantities oc 1/71. As these corrections are subleading (in 
a s or 1/rrib), however, a small shift in 7\ has only very minor impact on the branching 
ratios. The absolute scale for the branching ratios is set by the CP- and isospin-averaged 
branching ratio with the smallest experimental uncertainty, i.e. B — > K*^\ from Tab. 17.1} 
one finds: 



B{B -> K*j) = - \b{B ± -> K*±i) + — B(B° -» K* ^) \ = (41.6 ± 1.7) x 10~ e 

2 I T B o 



That is, we obtain a theoretical prediction for B[B — > Vj) as 



B(B -> Vi) 



th 



B{B -> Vj) 

B{B -»■ K*i) 



B[B K*i) 



th 



cxp 



(7.19) 



(7.20) 



where [. . .] th depends mainly on £y and only in subleading terms on the individual form 
factors Ti~* and T^ v . It is obvious that, except for these subleading terms, this 
procedure is equivalent to extracting an effective form factor Tf^* K *(0) \ ff from B — > 
K*j and using Tf^ v (0)\ s = Tf~* K *(0)\ /£y for calculating the branching ratios for 
B-^Vj. From (03]) we find 



off 



th cxp 

0.267 ±(101? ±(1006 



0.267 ±0.018 



(7.2i; 



where the theoretical uncertainty follows from the second uncertainty given in 07.15p . 
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Eqs. (17TB . dZHD and fT7~T8l then yield 



Tf~* p (0) ^ = 0.228±0.023, T?^ w (0) \ cS = 0.205 ± 0.021 



cff 



t b s ^k*^ = 0.245 ±0.024, Tf s ^(0) = 0.260 ± 0.036 . (7.22) 



cff 



elf 



Note that all effective form factors agree, within errors, with the results from LCSRs 
given in Tab. 17.81 which confirms the results obtained from this method; the crucial 
point, however, is that the uncertainties are reduced by a factor of 2 (except for Tf B ~*^). 
We would like to stress that the motivation for this procedure is to achieve a reduction 
of the theoretical uncertainty of the predicted branching fractions in B — > (p, u)'~f and B s 
decays. The effective form factors do not constitute a new and independent theoretical 
determination, but are derived from the experimental results for B — > K*^ under the 
following assumptions: 

• there is no new physics in B — > if*7^| 

• QCDF is valid with no systematic uncertainties; 

• LCSRs can reliably predict the ratio of form factors at zero momentum transfer. 
From (17.141) and (IT. 22j) . we then predict the following CP-averaged branching ratios: 







Ti Other 




B(B~ - 


* P 7) 


= (1.16 i'oSi '013 s ) X 


10~ 6 . 


B(B°- 


-A) 


= (0.55 ±0.11 ±0.07) x 


10" 6 , 


B(B° 




= (0.44 ±0.09 ±0.05) x 10" 6 


B(B S - 


>K*l) 


= (1.26 ± 0.25 ± 0.18) x 


10~ 6 , 


B(B S 




= (39.4 ± 10.7 ±5.3) x 


10" 6 , 



(7.23) 

where the first uncertainty is induced by the effective form factors and the second includes 
the variation of all inputs from Tab. 17.81 except for the angle 7 of the UT, which is 
fixed at 7 = 53° ll The total uncertainty in each channel is ~ 20%, except for B s — > 
07, where it is 30%. The results for p and ui agree very well with those of BaBar, 



6 Which is motivated by the results from inclusive B —* X s j decays [186]. 

7 The value of the UT angle 7 in Tab. l7.8l comes from Belle's Dalitz-plot analysis of the CP asymmetry 
in B~ — > (Kgir + ir~)K~ , with Kgir + ir~ [188] being a three-body final state common to both D° and D° . 
Other determinations all come with theoretical uncertainties and/or possible contamination by unresolved 
new physics, so we take this result as a reference point. 
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Tab. 17.11 but less so with the Belle results, although present experimental and theoretical 
uncertainties preclude a firm conclusion. Our prediction for B s — > 07 is well below the 
current experimental bound 120 x 10~ 6 [27]. A branching ratio of the size given in ( 17.231) 
implies that (9(10 3 ) B s — > 07 events will be seen within the first few years of the LHC. 

In Tab. 17.51 we detail the contributions of individual terms to the branching ratios. In all 
cases B is dominated by the QCDF contribution, with WA most relevant for B~ — > p~^y. 
This is expected as WA enters with the large Wilson coefficient Ci ~ 1. The effect is 
extenuated by long-distance (LD) photon emission, which itself is compensated by soft- 
gluon emission. The other channels follow a similar pattern, although the size of the 
effects is smaller. 



B x 10 6 


QCDF 


+ WA (no LD) 


+ WA (inc. LD) 


+ soft gluons 


B~ -> p~7 


1.05 


1.17 


1.11 


1.16 


B° p° 7 


0.49 


0.53 


0.53 


0.55 


B° ^uo-i 


0.40 


0.42 


0.42 


0.44 


B~ -> K*-*y 


39.7 


38.4 


38.3 


39.4 


BO K*°-f 


37.1 


39.7 


39.9 


41.0 


B° s -> i?*°7 


1.12 


1.22 


1.23 


1.26 


B° s - 07 


34.6 


38.2 


38.3 


39.4 



Table 7.5: Contributions to CP-averaged branching ratios, using effective form factors and 
central values of all other input parameters, Tab. 17.81 (in particular 7 = 53°). LD stands for long- 
distance photon-emission contributions. Each column labelled "+X" includes the contributions 
listed in the previous column plus the contribution induced by X. The entries in the last column 
are our total central values. 

We would like to close this section by making explicit the dependence of the three B — > 
(p, uo) r j branching ratios on 7. In Fig. 17. 31 we plot these branching ratios, for central values 
of the input parameters, as functions of 7. We also indicate the present experimental 
results from BaBar [189], Tab. 17.11 within their la uncertainty. 



Ill 



7 



7 




Figure 7.3: CP-averaged branching ratios of B — > (p, iv)-f as function of UT angle 7, using 
the effective form factors and central values of other input parameters, (a): — > /0 =t 7, (b): 
B° — ► p°7, (c): B° — ► wy. The boxes indicate the la experimental results from BaBar 
[189], Tab. 17.11 Note that the resulting value of 7 from the average of all three channels is 
7 = (61.0+^;^(exp)+^)° - see Section [73J 



7.4.2 Isospin Asymmetries 

The asymmetries Ai(p), Aj(K*), and A(p,u) are given by 

Y(B° -> K* ^) - T(B ± — > K* ±r A 
A^K*) = 1) 1} ^— ^, (7.26) 

where the partial decay rates are CP-averaged. Let us first discuss A(p,u) and Ai(p) 
which are relevant for the experimental determination of BiB — > (p,(jj)j), which in turn 
is used for the determination of |Kd/Ks| (or 7), see Section I7T51 The present experimental 
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7 


40° 


50° 


60° 


70° 


Mp) 


-(5.3 ±6.9)% 


(0.4 ±5.3)% 


(5.7 ±3.9)% 


(10.5 ±2.7)% 



Table 7.6: Isospin asymmetry Aj(p) for different values of 7. 



statistics for b — > dj transitions is rather low, so the experimental value of BiB — > (p, w)7) 
is obtained under the explicit assumption of perfect symmetry, i.e. r(S ± — > p ±r y) = 
2T(B° — > p°7) = 2T(B° — > wy). In reality, the symmetry between p° and uj is broken by 
different values of the form factors, and isospin symmetry between neutral and charged p is 
broken by photon emission from the spectator quark, the dominant mechanism of which 
is WA. From the formulas for individual branching ratios, Eq. f)7.14p . and the various 
contributions to the factorisation coefficients a^ L f R \, we find 

th. 

A(p, oj) = -0.20 ± ffo^ . (7.27) 

The uncertainty is dominated by that of the form factor ratio T 1 B ~* w (0)/T 1 B ~ >p (0) = 0.90± 
0.05§ The dependence on all other input parameters is marginal. The result (I7.27P is not 
compatible with the strict isospin limit A(p,u) = 0. Aj(p), on the other hand, is very 
sensitive to 7, whereas the form factors drop out. It is driven by the WA contribution 
and, in the QCDF framework, vanishes if WA is set to zero. In Fig. 17.4( a) we plot Aj(p) 
as function of 7, including the theoretical uncertainties. As suggested by the findings 
of Ref. [177], these results are not expected to change considerably upon inclusion of the 
non-factorisable radiative corrections of Fig. 17.2( c). In Tab. 17.61 we give the corresponding 
results for several values of 7, together with the theoretical uncertainty. Our result agrees 
very well with that obtained by the BaBar collaboration: Ar(p)B a Bar = 0.56 ±0.66 [189]. 

Aj(K*) was first discussed in Ref. [176], including power-suppressed 0(a s ) corrections 
which unfortunately violate QCDF, i.e. are divergent. It is for this reason that we decide 
to drop these corrections and include only leading-order terms in a s . We then find 

(J. NLO^LO /s other 

Ai(k*) = (5.4±^j?± 'ofT ± / olT±'a6 N )% 

= (5.4 ±1.4)%, (7.28) 



8 Note that this result is dominated by the ratio of decay constants given in Tab. 17.81 and discussed in 
Ref. [70] . The experimental results entering these averages have a large spread which may cast a shadow 
of doubt on the averaged hnal branching ratios for (p°,uj) — > e + e _ quoted by PDG [27]. 



113 



C 12. 




Figure 7.4: Left panel: isospin asymmetry Aj(p), Eq. (I7.25p . as function of the UT angle 
7. Solid line: central values of input parameters; dashed lines: theoretical uncertainty. Right 
panel: Aj(K*), Eq. (|7.26p . in percent, as function of the ratio r = ag/cig M of the combination 
of penguin Wilson coefficients a§ = Cq + C5/3. Solid line: central value of input parameters, 
dashed lines: theoretical uncertainty. The box indicates the present experimental uncertainty 
and the straight black lines the theory uncertainty in r. 

where NLO LO denotes the uncertainty induced by switching from NLO to LO accu- 
racy in the Wilson coefficients and "other" summarises all other sources of theoretical un- 
certainty. As can be inferred from the entries in Tab. 17.11 the present experimental result 
is A/(i\"*) cxp = (3.2±4.1)%. In Ref. [176] it was pointed out that A^K*) is very sensitive 
to the values of the Wilson coefficients Cff L in the combination = Cf BL + C<? BL /3. In 
the SM, varying the renormalisation scale as n = m^m^) ± 1 GeV and switching between 
LO and NLO accuracy for the Wilson coefficients, one has = —0.039 ± 0.008, which 
actually induces the bulk of the uncertainty in Eq. (17.281) . In Fig. 17.4( b) we plot Ai(K*) 
as function of a 6 /<2g M , with a| M = —0.039. The figure clearly indicates that, although 
there is presently no discrepancy between theoretical prediction and experimental result, 
a reduction of the experimental uncertainty of Aj(K*) may well reveal some footprints of 
new physics in this observable. 

7.4.3 CP Asymmetries 

The time-dependent CP asymmetry in B° — > V ^ is given analogously to Eq. (11.161) as 

Acp(t) = S{Vj) sin(Am D t) - C{Vj) cos(Am D t) . (7.29) 

The above equation is technically only valid for Ar = and while this is a good as- 
sumption for B® decays, it is not so for B® decays. Although Eq. (17.291) can easily be 
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adapted to non-zero Ar s we refrain from doing so: the whole point in calculating the 
CP asymmetry is not so much to give precise predictions for S and C, but rather to 
exclude the possibility of large corrections to the naive expectation S ~ mo/mi,. With 
this is mind, small corrections from a non-zero Ar s are irrelevant. The time-dependent 
CP asymmetries are given in terms of the left- and right-handed photon amplitudes 
by 

21m U(AIA L + A R A R )) 

(7.30) 

With Al,r and Al,r as given in (16. lip . The indirect CP asymmetry S(Vj) relies on 
the interference of both left- and right-helicity amplitudes and vanishes if one of them is 
absent; it thus probes indirectly the photon helicity. The direct CP asymmetry C(V^y) 
is less sensitive to Ar, but very sensitive to the strong phase of Al and vanishes if the 
radiative corrections to a^ CDF , Eq. (17.21) . are neglected. As the accuracy of the prediction 
of strong phases in QCDF is subject to discussion, and in any case C(Vj) is less sensitive 
to new physics than S(V , -f), we shall not consider direct CP asymmetries in this analysis. 

Let us briefly discuss the reason for the expected smallness of 5". In the process b — > Z>y, 
in the SM, the emitted photon is predominantly left-handed in b, and right-handed in b 
decays. This is due to the fact that the dominant contribution to the amplitude comes 
from the chiral-odd dipole operator Qj. As only left-handed quarks participate in the 
weak interaction, an effective operator of this type necessitates, in the SM, a helicity 
flip on one of the external quark lines, which results in a factor m;, (and a left-handed 
photon) in b R — > Dl'Jl and a factor mo (and a right-handed photon) in b^ — > Dr'Jr. 
Hence, the emission of right-handed photons is suppressed by a factor mrj/mj, which 
leads to the QCDF prediction (I6.13P for atf R . The interesting point is not the smallness 
of the CP asymmetry per se, but the fact that the helicity suppression can easily be 
alleviated in a large number of new physics scenarios where the spin flip occurs on an 
internal line, resulting in a factor rrii/mb instead of mo/mf,. A prime example is left- 
right symmetric models [191], whose impact on the photon polarisation was discussed in 
Refs. [180,181,192]. These models also come in a supersymmetric version whose effect 
on b — > S7 was investigated in Ref. [193]. Supersymmetry with no left-right symmetry 
can also provide large contributions to b — > D~/r, see Ref. [194] for recent studies. Other 
potential sources of large effects are warped extra dimensions [195] or anomalous right- 
handed top couplings [196]. Unless the amplitude for b — > D^ R is of the same order as 
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the SM prediction for b — > D'yi, or the enhancement of b — > -D7_r goes along with a 
suppression of b — > -D7L, the impact on the branching ratio is small, as the two helicity 
amplitudes add incoherently. This implies there can be a substantial contribution of new 
physics to b — > D7 escaping detection when only branching ratios are measured. 

We can calculate S directly from (I7.30p and obtain, making explicit the contributions 
from different sources: 





m D /m b 


soft 

LD WA gluons 




SQry) = 


('ooT 


+ ^02^ + ^20^ ± 1.6)% 


= (0.2 + 1.6)%, 


S(ujj) = 


(0.01 - 


- 0.08 + 0.22 ± 1.7)% 


= (0.1 + 1.7)%, 


S(K^) = 


-(2.9- 


+ 0.6 + 1.6)% 


= -(2.3 + 1.6)%, 


S(K^) = 


(0.12 4 


-0.03 + 0.11 + 1.3)% 


= (0.3 + 1.3)%, 


S(<h) = 


(0 + 


+ 5.3 + 8.2) x 10^ 2 % 


= (0.1 + 0.1)%. 



(7.31) 

Including only the helicity-suppressed contribution, one expects, for B — > K*j, neglecting 
the doubly Cabibbo suppressed amplitude in A« 

SiK^) | no soft gluons = -2 -i sin <p d « -2.7% . (7.32) 

'ii'b 

For B s — > 07, one expects the CP asymmetry to vanish if the decay amplitude is pro- 
portional to \[ s \ which, at tree-level, precludes any contributions of type sm(<f) s )m s /mt, 
and also any contribution from WA. This is because the mixing angle <p s is given by 
arg[(A^) 2 ], Eq. (jl.l5p . and the interference of amplitudes in (17.301) also yields a factor 
(A^) 2 , if the individual amplitudes are proportional to A^ or (A$ )*, respectively; this is 
indeed the case for the helicity-suppressed term m s /mb induced by the operator Q-j and 
the WA contributions to a^ R (<f>), Eqs. f!7.6p and f)7.10p . so that the phases cancel in (17.3 



The actual results in f)7.3ip disagree with the above expectations because of the contribu- 
tions from soft-gluon emission, which enter Sf R . Moreover, for S(({rf) this is because the 
soft-gluon emission from quark loops is different for u and c loops so that a? R 7^ a% R and 
hence Ar (Al) is not proportional to A^ ((A^)*). Note that a substantial enhancement 
of S(4ry) by new physics requires not only an enhancement of |Ar| (and | Al|), but also the 
presence of a large phase in (17.301) : this could be either a large B s mixing phase which will 
also manifest itself in a sizable CP violation in, for instance, B s — > J/ip<ft, see Ref. [28,197]; 
or it could be a new weak phase in Ar (and Al)] or it could be a non-zero strong phase in 
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one of the coefficients. Based on the light quark loop results there is not much scope 
for a large phase in dj R (whose contribution is, in addition, doubly Cabibbo suppressed), 
but the situation could be different for , where only the leading-order term in a l/m c 
expansion is included, which does not carry a complex phase [70] . It is not excluded that a 
resummation of higher-order terms in this expansion will generate a non-negligible strong 
phase — which is not really relevant for our results in Eq. ( I7.3ip . but could be relevant 
for the interpretation of any new physics to be found in that observable. For S(K* / -/), on 
the other hand, no new phases are required, and any enhancement of \Ar\ (and by 
new physics will result in a larger value of S(K*'y). 

For all S except S(K*^), the uncertainty is entirely dominated by that of the soft-gluon 
emission terms l U)C — l UyC , whose uncertainties have been doubled with respect to those 
given in Tab. 17.41 The smallness of S((p, oj)^f) is due to the fact that the helicity factor 
is given by rrid/mi, (we use m u ^/m s = 1/24.4 from ChPT). For K*, the suppression 
from the small mixing angle is relieved by the fact that both weak amplitudes in Xjj 
contribute so that the CP asymmetry is comparable with that of p and uj. Despite the 
generous uncertainties, it is obvious that none of these CP symmetries is larger than 4% 
in the SM, which makes these observables very interesting for new physics searches. The 
present experimental result from the B factories, S(K*j) = —0.28 ± 0.26 [37], certainly 
encourages the hope that new physics may manifest itself in that observable. While a 
measurement of the b — > d CP asymmetries is probably very difficult even at a super- 
flavour factory, S(K*-~f) is a promising observable for B factories [6], but not for the LHcE 
B s — > (j>{— > K + K~) / -f, on the other hand, will be studied in detail at the LHC, and in 
particular at LHCb, and any largely enhanced value of 5"(07) will be measured within the 
first years of running. 



7.5 Extraction Of CKM Parameters 

Let us now turn to the determination of CKM parameters from the branching ratios 
determined in Section [7.4.11 In this context, two particularly interesting observables are 

_ B(B^(p,uj) 1 ) BjB -> pry) 



9 K* has to be traced via its decay into a CP eigenstate, i.e. Ksn°. Neutrals in the final state are not 
really LHC's favourites. 
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given in terms of the CP- and isospin-averaged branching ratios of B — > (p, uj) r y and 
B — > p7, respectively, and B — > i^*7 decays. -R p / W has been measured by both BaBar 
and Belle [189,190], a first value of _R P has been given by BaBar [189]. The experimental 
determinations actually assume exact isospin symmetry, i.e. T(B ± — > p 7) = 2T(B° —> 
p°7), and also T(B° — > p°7) = r(5° — > ^7); and as we have seen in Section 17.4.21 
these relations are not in fact exact. Hence, the present experimental results for R p / U are 
theory-contaminated. As the isospin asymmetry between the charged and neutral p decay 
rates turns out to be smaller than the asymmetry between p° and u, it would actually be 
preferable, from an experimental point of view, to drop the uo channel and measure R p 
instead of R p / U , as done in the most recent BaBar analysis on that topic [189]. We will 
give numerical results and theory uncertainties for both R p / U and R p . 

One parametrisation of R p / U often quoted, in particular in experimental papers, is 



R 



Vtd 

V ta 



1 — m 2 /m% \ 3 1 
l-mj c ./m 2 B J Q 



(7.34) 



with AR = 0.1 ± 0.1 [198] and again assuming isospin symmetry for p and u. This 
parametrisation creates the impression that AR is a quantity completely unrelated to 
and with a fixed value independent of |Vtd/Ks|- We would like to point out here that this 
impression is wrong: AR contains both QCD (factorisable and non-factorisable) effects 
and such from weak interactions. In Ref. [183] AR is expressed in terms of the factorisation 
coefficients a% L , assuming isospin symmetry for p° and u, as 



1 + AR 



l 7L 



(P) 



'7L 



Re (5a± + 5a c 



R\ 



Rb cos 7 



2Rb cos 7 



{\6a ± \ 2 + \5a \ 2 ) 



Rl 



2R b cos 7 



R 2 



(7.35) 



with 5ao,± = «7i,(p ' ± )/a7L(p ' ± ) _ 1- Eq. (17.351) shows explicitly that AR depends both 
on QCD (5a± t o) and CKM parameters (^,7). The point we would like to make is that 
the calculation of AR requires input values for R b and 7. Once these parameters (and 
the Wolfenstein parameter A) are fixed, however, |Vfd/Vf S | is also fixed and given by 



Vtd 
V ts 



AW1 - 2R b cos 7 + R 



1 + -(1 

2 v 



2R b cos7)A 2 + C(A 4 



(7.36) 



Hence, as iV^/Vy and (Rb, 7) are not independent of each other, it is impossible to extract 
from (17.341) with a fixed value of AR. Of course R p j w and R p of (17. 33ft can be 
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used to extract information about CKM parameters, but in order to do so one has to 
settle for a set of truly independent parameters. Based on (I7.36p . one can exchange, say, 
7 for |Vt(j/T4s|@ So we can either consider Ry as a function of the CKM parameters Rb 
and 7 (let us call this the 7 set of parameters) or as a function of Rb and (V^/Vy (to be 
called the \V tx \ set). Using the 7 set, a measurement of Rv{l, Rb) allows a determination 
of 7, whereas Ry(\ Vtd/Vt s \, Rb) allows the determination of |Kd/Ks|- In either case, the 
simple quadratic relation (17.341) between Ry and |Kd/Ks| becomes more complicated. 

In Figs. 17.51 and 17.61 we plot the resulting values of |V^/V^| 2 and 7, respectively, as a 
function of Ry. Although the curve in Fig. 17.5( a) looks like a straight line, as naively 
expected from ( 17. 34ft . this is not exactly the case, because of the dependence of AR on 
|Vtd/T4s|. In Fig. 17.5( b) we plot AR for the \V tx \ set of parameters. The dependence of 
AR on I Vtd/Vts| is rather strong. Apparently indeed AR = 0.1 ±0.1 in the expected range 
0.16 < I Vtrf/Vts| < 0.24, but this estimate does not reflect the true theoretical uncertainty 
which is indicated by the dashed lines in the figure. 

It is now basically a matter of choice whether to use R p / W to determine |Vtd/Vt s | or 7. 
Once one of these parameters is known, the other one follows from Eq. (17.361) . In Fig. 17.61 
we plot 7 as a function of R p / U , together with the theoretical uncertainties. In Fig. 17.71 
we also compare the central values of R p / W with those of R p , as a function of |Vtd/Vt s |. 
Although the difference is small, R p is expected to be larger than R p / U . R p / U and R p are 
dominated by the uncertainties of £ p and as discussed in Ref. [183], a reduction of this 
uncertainty would require a reduction of the uncertainty of the transverse decay constants 
f$ of p and K*. With the most recent results from BaBar, R p/u = 0.030 ± 0.006 [189], 
and from Belle, R p/u) = 0.032 ± 0.008 [190], we then find 



BABAR: 
Belle: 



V t 



td 



V ts 

v td 



exp t jj exp th 

0.199^±o"oi4 <-> 7 = (61.0±ii;g 



V 



ts 



0.207 ±8:SH ±8:8it - i = (65.7i£?iS3)° • (7.37) 



These numbers compare well with the Belle result [188] from tree-level processes, 7 = 
(53 ± 20)°, quoted in Tab. 17.81 an d results from global fits [31]. We also would like to 
point out that the above determination of 7 is actually a determination of cos 7, via 
Eq. (17.36j) . and implies, in principle, a twofold degeneracy 7 «-> 2tt — 7. This is in contrast 



10 Strictly speaking, ()7.36p only fixes COS7 as function of |Vtd/Vt s |, leaving a twofold degeneracy of 7. 
Eq. ()7.35j) . however, only depends on cos 7, so that indeed one can unambiguously replace 7 by |Vtd/Vt s |. 
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0.025 0.03 0.035 



0.16 



0.2 0.24 
\V td /V t3 \ 



.28 



Figure 7.5: Left panel: |^d/Vf S | 2 as function of R p / W , Eq. (|7.33p . in the \Vt x \ basis - see text. 
Solid line: central values. Dash-dotted lines: theoretical uncertainty induced by ^ p = 1.17±0.09, 
(|7.16p . Dashed lines: other theoretical uncertainties, including those induced by \ V u b\, \ V c b\ and 
the hadronic parameters of Tab. 17.81 Right panel: AR from Eq. (|7.35p as function of IV^/V^I 
in the \ Vt x \ basis. Solid line: central values. Dashed lines: theoretical uncertainty. 



7 80. 
70 . 
60. 
50 . 
40. 
30. 

0.025 0.03 .0 35 0.04 

Figure 7.6: The UTangle 7 as function of R p / W in the 7 set of CKM parameters. Solid 
lines: central values of input parameters. Dash-dotted lines: theoretical uncertainty induced by 
£p = 1.17 ± 0.09. Dashed lines: other theoretical uncertainties. 

R 
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0.04 
0.03 
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\V id /V ts \ 

Figure 7.7: Central values of R p / U (solid line) and R p (dash-dotted line) as functions of 
\Vtd/V u \. 
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to the determination from B — > D^*'K^*' in [188], which carries a twofold degeneracy 
7 «-> 7i + 7. Obviously these two determinations taken together remove the degeneracy 
and select 7 « 55° < 180°. If 7 « 55° + 180° instead, one would have \V td /V ts \ « 0.29 
from (I7.36p . which is definitely ruled out by data. Hence, the result (17.371) confirms the 
SM interpretation of 7 from the tree-level CP asymmetries in B — > D^K^*\ 
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tb±/t~b» 


tb°/tb° 


1.530(9) ps 


1.071(9) 


0.958(39) 



Table 7.7: B lifetimes from HFAG [37]. 



CKM parameters and couplings 



A [27] 


\V cb \ [202] 


\v ub \ 


7 [188; 


a s (m z ) [27] 


«QED 


0.227(1) 


42.0(7) x 10" 3 


4.0(7) x 10" 3 


(53 ±20)° 


0.1176(20) 


1/137 


B parameters 


f Bq [184] 


f Bs [184] 


Ab,W [183] 








200(25) MeV 


240(30) MeV 


0.51(12) GeV 


0.6(2) GeV 


2.2 GeV 




p parameters 


f P 


J P 




4(p) 






216(3) MeV 


165(9) MeV 





0.14(6) 


0.27(4) 




uj parameters 


u 






4M 






187(5) MeV 


151(9) MeV 





0.15(7) 


0.25(4) 




K* parameters 




fx* 


af(K*) [61 


4(#*) 


rf^*(o) 


lf^*(0) 


220(5) MeV 


185(10) MeV 


0.04(3) 


0.15(10) 


0.31(4) 


0.29(4) 


<fi parameters 


f<t> 




4(0) 


4(0) 


2f*^(0) 




215(5) MeV 


186(9) MeV 





0.2(2) 


0.31(4) 




quark masses 


m s (2GeV) [199] 


m b (m b ) [202] 


m c (m c ) [200] 


m t (m t ) [201] 


100(20) MeV 


4.20(4) GeV 


1.30(2) GeV 


163.6(2.0) GeV 



Table 7.8: Summary of input parameters. The value of \V ub \ is an average over inclusive and 
exclusive determinations and the result from UTangles Refs. [31,37,203]. None of our results 
is very sensitive to \V ub \. For an explanation of our choice of the value of the UT angle 7, see 
text. Xb s is obtained from A# 9 , see Eq. (|7.4|) . The vector meson decay constants fy, fy are 
discussed in Ref. [70]; the values of the Gegenbauer moments af~ are compiled from various 
sources [54, 55, 59, 183] and include only small SU(3)F-breaking, in line with the findings for 
pseudoscalar mesons [57]. The form factors T\ are updates of previous LCSR results [112], 
including the updated values of the decay constants /p, w ,0 and of a\{K*) [61,62]. All scale- 
dependent quantities are given at the scale (i = 1 GeV unless stated otherwise. 
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Chapter 8 

Summary and Conclusions 



This thesis has consisted of three main analyses centred on the investigations and determi- 
nations of meson light-cone distribution amplitudes. We have seen how the determinations 
of decay observables in B decays are reliant on the sound understanding of both theoreti- 
cal and experimental uncertainties with the work presented in this thesis striving towards 
the former. To summarise: 

We began, in Chapter with a brief introduction defining the QCD Lagrangian, dis- 
cussing CP violation and the AB = 1 effective Hamiltonian. 

In Chapter [2] we investigated the structure of vector mesons distribution amplitudes to 
twist-3 accuracy. We included all SU(3)p-breaking and G-parity violating effects. The 
QCD equations of motion were implemented to unpick the interwoven relations between 
the distribution amplitudes ultimately expressing the two-particle twist-3 distribution 
amplitudes in terms of the three-particle twist-3 and two-particle twist-2 distribution am- 
plitudes. The equations of motion result in integral equations which are readily solved 
order-by-order in conformal spin and to the order considered all the distribution ampli- 
tudes are then expressed by a small number of non-perturbative parameters. Finite quark 
mass effects appear in the equation of motion and therefore impact the two-particle twist- 
3 distribution amplitudes (12.48112301) . Such effects also cause mixing between the twist-3 
hadronic parameters under renormalisation scale evolution, see Eq. ( 12.631 ). 

In Chapter [3] we discussed the methods of QCD sum rules (the SVZ method) and QCD 
sum rules on the light-cone. We outlined the procedures with example correlation func- 
tions and ended the chapter with an example calculation of the a s corrections to the gluon 
condensate contribution to a K meson sum rule. The calculation made use of the back- 
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ground field technique and served to illustrate the calculation of radiative corrections to - 
and extraction of - vacuum condensates in the SVZ method. The result of the calculation 
is in conflict with that in the literature, see Eqs. ( 13.41D and ( 13.42H . 

In Chapter H] we determined the leading hadronic parameters defined in Chapter [2] via 
SVZ sum rules. We calculated the three-particle twist-3 parameters to NLO in conformal 
spin, also including all G-parity violating terms and finite strange quark mass effects. 
The determination of the twist-3 parameters is new for K* and <fi. The results for the p 
agree within uncertainties with previous determinations and are presented in Tabs. 14.11 
and 14.21 We also calculate 0(a s ) and 0{m 2 s ) corrections to the quark condensate for the 
sum rules for a)n (V), which for n = 2 is the first non-trivial Gegenbauer coefficient of 
the G-even particles p and <fi. We add this contribution to the existing sum rules taken 
from the literature and update the value of a| ,_L (0) which we find to be consistent with 
that found for K* and p; a,i 2 (V) — a i2(^0 within uncertainties. The results find direct 
application in QCD factorisation descriptions of B — > V decays, and the light-cone sum 
rule analyses of B — > V transition form factors. 

In Chapter [3] we calculated the form factors of B — > rj' semileptonic transitions from light- 
cone sum rules, including the gluonic singlet contributions. We built upon the previous 
light-cone sum rule determination of the B — > rj form factor by casting the calculation 
consistently within the phenomenologically motivated 77-77' mixing scheme of Refs. [117, 
118]. We found that, as expected, these contributions are more relevant for than for 
/+ and can amount up to 20% in the former, depending on the only poorly constrained 
leading Gegenbauer moment B 2 of the gluonic twist- 2 distribution amplitude of 77'. The 
numerical results, with each contribution listed separately, are given by Eqs. (15.451) and 
(15.461) . Consequently, it seems unlikely that the large exclusive B — > rj K and inclusive 
B — > rj'X branching ratios can be explained by a large B 2 , as it would have to assume a 
very extreme value. We also found that the form factors are sensitive to the values of the 
twist-2 two-quark Gegenbauer moments a 2 ' n which, given the uncertainty of independent 
determinations, we have set equal to a 2 , see Fig i5.6l 

The ratio of branching ratios BiB — > rjev)jB(B — > rjeu) is sensitive to both a 2 and B 2 and 
may be used to constrain these parameters, once it is measured with sufficient accuracy, 
see Fig. 15.81 The extraction of \ V u b\ from these semileptonic decays, in particular B — > rjeu, 
with negligible singlet contribution, although possible in principle, at the moment is 
obscured by the lack of knowledge of a^- We would also like to stress that, in the framework 
of the quark-flavour mixing scheme for the 77-77' system as used in this analysis, B — > 77' 
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transitions probe only the r] q component of these particles. The rj s component could be 
probed directly for instance in the b — > s penguin transition B s r]'£ + £~, although such 
a measurement would also be sensitive to new physics in the penguin diagrams. 

In Chapter [U] we discussed the QCD factorisation (QCDF) approach of Refs. [10,11] and 
its application to the radiative B decays B — > V'j of Refs. [16,17]. We discussed the 
appearance of distribution amplitudes in the factorisation formulas and focused on the 
leading contributions to the B — > V'j decays. 

In Chapter [7] we performed a phenomenological analysis of the radiative B decays to 
vector mesons B — > Vj, using the framework discussed in Chapter El We investigated 
the most relevant power-suppressed corrections to the QCDF predictions for the radiative 
decays B u ^ — > (p, u, K*)j and B s — > (<f), K*)^. We use the QCDF framework presented in 
Refs. [16, 17] in which we find use for the twist-2 DA parameters determined in Chapter HI 
Besides the leading QCDF contributions we included long-distance photon emission and 
soft-gluon mission from quark loops. These effects, although formally ~ l/m,b with respect 
to the leading contributions, augment the QCDF predictions for the branching ratios, CP 
and isospin asymmetries. 

The impact of the power-suppressed corrections on the branching ratios is found to be 
very small, with the exception of the weak annihilation contributions to B — > p ±r y which 
are large due to a large combination of Wilson coefficients C-x + Ci/3 = 1.02 and no 
CKM-suppression. Moreover, long-distance photon emission also impacts most here, see 
Eq. (17. lip . An explicit break down of the results are given in Tab. 17.51 

The isospin asymmetries A(p,u), Ar(p) and Ai(K*) are driven by weak annihilation and 
long-distance photon emission contributions. We found a non-zero asymmetry A(p,u) = 
—0.20 ± 0.09 which suggests the explicit assumption of perfect symmetry, i.e. r(B ± — ► 
p ±, ~f) = 2T(B° — > p°7) = 2T(B° — > ury) used to obtain the experimental value of B(B — > 
(p, uj)^) is not so well justified. We found Ai(p) to depend strongly on the UT angle 7, 
as shown in Tab. 17.61 With our central value of 7 = 53° (see Tab 17.81) our result agrees 
very well with the BaBar result A I (p) BsLBar = 0.56 ± 0.66 [189]. For Ai(K*) we found 
a result consistent with the experimental result Ai(K*) cxp = (3.2 ± 4.1)% and, via its 
sensitivity to the Wilson coefficient combination C5 + Cq/3 conclude that a reduction 
in the experimental uncertainty may uncover signs of new physics contributing to these 
Wilson coefficients, see Fig. 17.41 

The indirect CP asymmetries S(Vf) are caused by the interference between the ampli- 
tudes describing the production of left and right-handed photons, see Eqs. (16.81) and (17.30j) . 
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The right-handed amplitude is suppressed by mj)/m(, with respect to the left-handed one 
for B = bq decays (and vice versa for B decays). Due to this natural suppression in 
the SM we expect the CP asymmetries to be small, and this suppression can be relieved 
by many new physics senarios. We investigated the soft-gluon effects arising from soft 
heavy and soft quark loops. The calculation of these contributions makes use of the three- 
particle twist-3 DA parameters determined in Chapter HI They contribute to both the left 
and right-handed amplitudes, and so may also relieve to SM suppression. We found that 
although they do divert the results from the values naively expected, there is no scope for 
a large enhancement due to these power-suppressed contributions. The results are given 
in Eq. (E3I]). 

Finally, using the most recent results from BaBar and Belle, we extracted the CKM 
parameter ratio |Vtd/^L| and equivalently the UT angle 7 from the ratio of branching 
ratios R p / U . The results are 



BABAR: 
Belle: 



Vtd 
V u 
Vtd 



cxp 



th 



cxp 



th 



0.199 ±8:825 =•= 0.014 



n 9D7 +0-028 +0.014 
U.ZUf _ .033 -0.015 



7 = (61.0 



+13.5 +8.9 no 
-16.0 -9.3 I 



7 



(at; 7+17.3 +8.9\o 
\VO. I _ 20 .7 -9.2) 



(8.1) 



and agree well with the Belle result 7 = (53 ± 20)° obtained from tree-level processes, 
and results from global fits [31]. The result confirms the SM interpretation of 7 from the 
tree-level CP asymmetries in B — > D^*'K^*' . 
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Appendix A 
Light-cone Co-ordinates 



To perform the light-cone expansion one relate the meson's 4-momentum P^, polarisation 
vector and the coordinate x M to two light-like vectors p^ and z^. We have the usual 
relations 



P =0, 



0. 



(a.i; 



and 



m 



so that the limit m 2 K * 
that 



P ■ e (A) = 0, 



(A.2) 



gives p — > P and x 2 — > gives 2 — > x. From this it follows 



A* * 



X • P 



x-py 



2 /rn 2 



arm 



A'* 



2^ 



K* 



p ■ Z 



X, 



2^2 



x m 



A* 



4(z ■ pf 



x 



2 p ■ z 



0(x 4 
(A. 



The meson's polarization vector can be decomposed into projections onto the two 
light-like vectors and the orthogonal plane 



p ■ z 



P» 



p ■ z 



, (A) 



m 



A* 



2p ■ z 



2 U +e 



(A) 



(e (A) • x) ^ X ' ^ x v m K* 



[x ■ py 



M""A* , (A) 
2™2 C ±M ' 



arm 



(A.4) 



A* 
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We also need the projector g^ u onto the directions orthogonal to p and z 



div = 9^u — (p^Zv + PuZp) . (A.5) 

p • z 



Some useful scalar products are 



z ■ P = z 



■ p = ■ P) 2 — x 2 m 2 K * 



p . e W = _<i z . e (A) 

y 2pz 

z . e W = x -e (A) . (A.6) 

Will use the notations 

a z = a^, b p = b^, </=l^, d^ = g^ v d v , (A.7) 

for arbitrary Lorentz vectors a M , b^, c M and and 

x» = x.n^ + x + n^ + , (A. 8) 

for null unit vectors n 2 = n 2 = and n ■ n = 1. The following notation is also used: 

a_ = , a = a M - a + Zn . (A. 9) 

p ■ z p p ■ z 
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Appendix B 

Useful formulas for sum rule 
determinations 



B.l Loop Integrals 

Here we summarise the loop integrals needed for calculating the twist-3 correlation func- 
tions in Chapter HI At one loop, one has (z 2 = 0) [59] 

X f dw W D/2 ~ l -\l - w )DH+n-l-a ? 

Jo 

(B.l) 

where the integration measure is defined as dPk = i/(4ii) 2 \d L k\ and is an arbitrary 
numerical factor, which in the cases considered in Chapter H] is either v or v. One also 
needs the integral 



J 11 (i 2 ) c ((i - k) 2 ) d 

2=4 w, vT{c+d-D/2) f 1 



( _i)^ {k 2 ) UIM ( k -p)(k ■ z y 1 r(c)rW J duel{1 ' u)hk ' ZuDIM ^- u ) DI2+3 ~ c 

x[ due i{l ~ u)flk - z u D/2 ~ d (l-u) D/2 ~ 1+ ^ c (j + if l (l-u){k- z)) . (B.2) 
Jo 
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Two-loop integrals are obtained by combining the above one-loop integrals. 



B.2 Borel Subtraction 



To derive the sum rules from tt^. v , ir^.y and ir^: v we use the relation 



hm s [-q 2 - zOf = *° e(s) 

7r L J F{— a)r(l + a) 



(B.3) 



where s = —q 2 , to find the imaginary part. Using the following notation for the Borelisa- 
tion and continuum subtraction procedure 



B S ub [AT] 



dse- s/M2 -lm s X , 

7T 



(B.4) 



and the definition of the Borel transform (13.141) allows one to write the required results 

as 



B 



sub 



B sub [q 2 H-q 2 )] 
ln(—q 2 



(a - 1)!(M 2 )«~ 1 ' 

roo 

-M 4 + / dse- s/M2 s, 



B sub [\n(-q 2 )} = -M 2 + / dse~ s ' M 



Bsub 
Bsub 



ln(-g 2 ) 



so 



lE -\nM 2 + / dse~ s/M 

Jso s 



/•oo 

B sub [ln(-g 2 ) 2 ] = 2M 2 ( 7e - In M 2 ) + 2 / e" s/M2 In s 



(B.5) 



where 7# is Euler's constant. 

B.3 Input Parameters 

For the twist-2 and twist-3 DA parameter sum rule determinations of Chapter [4] we use 
the following input parameters: 



To evaluate the sum rules for the three-particle twist-3 DA parameters we use the following 
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\QQ) — I — u.z4 it u.ui ) vje v 

(qag s Gq) = m\ (qq) 




\ ss ) - I 1 ~ d 3j W<7/ 

larrn He\ — (1 A Mnrrn Hn\ 

\scrg s us) — [i — o^)\qag s (jq) 


(~ = (°- 012 ± °- 003 ) GeV4 






ml = (0.8±0.1)GeV 2 , 5 3 


= 0.2 ±0.2, <5 5 = 0.2±0.2 


m s (2 GeV) = (100 ± 20) MeV <- 
m 9 (/i) = m a (n)/R, 


-> m s (l GeV) = (133 ± 27) MeV 
i2 = 24.6 ± 1.2 


a s (M z ) = 0.1176 ±0.002 *- 




-> a s (lGeV) = 0.497 ±0.005 



Table B.l: Input parameters for sum rules at the renormalisation scale fi = 1 GeV. The value 
of m s is obtained from unquenched lattice calculations with Nf = 2 flavours as summarised 
in [204], which agrees with the results from QCD sum rule calculations [205]. fn q is taken from 
chiral perturbation theory [206]. a s (Mz) is the PDG average [27]. 



values of the continuum threshold so 

s J(p) = (1.3 ± 0.3) GeV 2 , sl(K*) — (1.3 ± 0.3) GeV 2 , 4(0) = (1 A ±0.3) GeV 2 , 
4(p) = (1.5 ±0.3) GeV 2 , 4(iT) = (1.6 ± 0.3) GeV 2 , s^((j)) = (1.7 ± 0.3) GeV 2 . 

(B.6) 

The threshold for the p channel is from [13]. 
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